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Abstract 

Recent  developments  in  communication  systems  technology  promise  to  greatly  improve 
the  performance  of  point-to-point  communications  for  both  commercial  and  tactical 
networks.  These  developments  include  the  use  of  electronically  steerable  antenna  arrays, 
space-time  multiple  input  multiple  output  (MIMO)  signal  processing  techniques,  and 
improved  techniques  for  error  correction.  In  this  project  we  will  address  the  challenging 
question  of  how  these  technological  developments  can  best  be  exploited  in  a  tactical 
networking  context,  where  signal  interference  and  channel  uncertainty  issues  have  a 
tremendous  impact  on  end-to-end  system  performance. 

Tactical  applications  pose  unique  requirements  for  the  network,  including  decentralized 
control  to  eliminate  single  points-of- failure,  vulnerability  to  jamming  and  electronic 
warfare,  and  mission  critical  latency  bounds  for  end-to-end  data  delivery.  Moreover,  a 
tactical  network  is  generally  composed  of  mobile  nodes  and  the  routing  protocols  must 
deal  with  a  range  of  node  mobilities  and  time  varying  channel  conditions.  Consequently 
this  project  is  focused  on  the  design  of  ad  hoc  networking  architectures  that  utilize 
MIMO  transmitters  and  receivers  at  each  node.  The  goal  of  this  program  is  to  define  the 
best  way  to  utilize  multiple  transmit  and  receive  antennas  at  each  node  to  improve  the 
robustness,  capacity,  and  quality  of  service  of  the  network. 
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1.  Introduction 

The  project  team  consists  of  fourteen  faculty  members  from  four  campuses  of  the 
University  of  California  (San  Diego,  Irvine,  Santa  Cruz,  and  Riverside);  Brigham  Young 
University,  Provo,  Utah;  and  McMaster  University,  Hamilton,  Ontario,  Canada.  The 
individual  faculty  members  were  selected  for  specific  areas  of  expertise  that  together 
spanned  the  wide  ranging  research  concentration  areas  defined  in  the  initial  BAA  for  the 
topic  area.  The  specified  concentration  areas  spanned  “the  physics  of  RF  propagation 
and  signal  processing;  the  electrical  engineering  of  antenna  array  design  and  electronics; 
computer  science  of  networking;  and  the  mathematics  of  information  and  control  theory”. 

The  specified  objective  of  this  MURI  topic  was  to  “create  network  protocols  and  signal 
processing  algorithms  necessary  to  implement  adaptive  beam  steering  and  spatial  channel 
reuse  in  mobile  wireless  communication  networks  with  the  specific  objective  of  enabling 
reuse  of  radio  channels  to  double  network  capacity  and  improve  protection  for  military 
communications.  The  research  should  also  result  in  the  science  that  will  allow  for  the 
decision  of  which  spatial  reuse  technique  to  use  (space-time  coding  (STC)  or  transmit 
beam  forming),  if  any,  based  on  topology  and  network  load. 

A  series  of  project  meetings  have  been  held  at  UCSD  and  UC,  Irvine  over  the  course  of 
this  project.  The  goal  of  these  meetings  was  to  first  define  a  hierarchy  of  problems  that 
require  resolution  and  then  outline  the  approaches  to  solve  these  problems  so  that  robust 
and  reliable  mobile  ad  hoc  networks  can  be  developed  for  tactical  applications.  These 
group  meetings  have  led  to  joint  work  between  Pis,  joint  publications,  co-advising  of 
graduate  students,  and  numerous  visits  between  subgroups  of  Pis  to  identify  viable 
approaches  for  solving  cross-layer  networking  problems. 

The  problem  definition,  research  issues  and  the  overall  project  research  goals  were 
defined  in  the  previous  interim  annual  reports  and  the  work  during  the  past  year  has 
focused  on  satisfying  these  established  goals. 

This  report  describes  the  specific  accomplishments  for  the  current  year  from  August  1 , 
2008  through  July  31,  2009.  During  this  period  we  have  made  significant  progress  in 
meeting  the  objectives  of  “providing  cross-layer,  energy  efficient  MIMO  signal 
processing  algorithms  for  mobile  multiuser  ad  hoc  networks  employing  directional 
antenna  arrays  and  space  time  coding  (STC)  for  tactical  applications.”  We  have 
successfully  developed  techniques  for  exploiting  beamforming  (BF)  and  STC  for 
improved  collision  avoidance.  Our  results  clearly  indicate  however,  that  the  full  potential 
for  the  integration  of  the  physical  layer  information  from  MIMO  antennas  into  the 
routing  and  scheduling  protocols  for  a  tactical  ad  hoc  network  will  only  be  achieved  by 
developing  a  fundamentally  new  approach  to  the  protocol  development  that  utilizes  the 
antennas  to  enable  multiple  simultaneous  transmissions  from  multiple  nodes  to  provide 
multi-packet  reception  capablilities  at  each  node.  The  results  show  that  the  network 
performance  will  be  optimized  by  utilizing  the  spatial  processing  gains  afforded  by  the 
antennas  to  manage  the  interference  in  a  controlled  fashion  rather  than  simply  trying  to 
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improve  existing  network  routing  and  seheduling  protoeols  to  reduee  the  number  of 
eollisions. 

In  addition,  the  results  elearly  show  the  advantages  of  eooperative  eommunieations 
between  nodes  using  either  relay  nodes  and/or  virtual  MIMO  nodes.  These  approaehes 
provide  energy  effieient  utilization  of  the  network  resourees  and  enable  the  development 
of  protoeols  with  signifieant  spatial  spreading  gain  even  for  nodes  with  limited  antenna 
resourees.  Sueh  eooperation  has  also  been  shown  to  provide  improved  spaee-time 
seheduling  in  terms  of  fairness  and  quality  of  serviee. 

The  utilization  of  MIMO  assets  forees  the  entire  question  of  channel  capacity  in  MIMO 
ad  hoc  networks  to  be  reexamined  as  a  function  of  how  the  information  is  disseminated. 
The  use  of  unicast,  multicast,  broadcast  and  various  forms  on  one-to-one,  one-to-many 
and  many-to-many  routing  can  be  done  more  efficiently  with  multiple  antennas  at  the 
nodes.  The  optimization  of  the  routing  protocols  depends  fundamentally  on  the  quality  of 
the  channel  state  information  (CSI)  that  is  available  at  the  nodes  of  the  network  however, 
and  significant  effort  has  been  made  to  define  the  best  transmission  strategies  as  a 
function  of  the  available  CSI.  A  testbed  has  been  developed  and  tested  to  evaluate  end- 
to-end  performance  of  different  topologies  in  multi-hop  ad  hoc  networks. 

The  relations  between  MIMO  signal  processing  and  network  coding  have  been  explored 
and  novel  rate  adaptation  techniques  that  can  trade  off  rate  for  coding  and  diversity  gain 
have  been  developed.  In  addition  network  metrics  such  as  information  efficiency  and 
throughput  have  been  defined  in  terms  of  physical  layer  parameters  such  as  the 
modulation  and  coding  schemes,  the  channel  statistics,  the  Doppler  rate  and  the  spatial 
diversity  order.  The  selection  of  the  appropriate  number  of  training  symbols  to  assure 
given  levels  of  throughput,  delay  and  efficiency  were  also  evaluated.  The  dependence  of 
the  information  efficiency  and  network  throughput  on  the  underlying  physical  parameters 
of  the  network  was  also  extended  to  random  multi-hop  networks. 

One  of  the  initial  barriers  to  crosslayer  integration  of  the  physical  layer  information  into 
the  routing  and  scheduling  protocols  was  the  significant  differences  between  the  time 
scales  associated  with  the  relative  time  scales  of  channel  variations  at  the  physical  layer 
that  occur  at  the  symbol  rate  and  the  temporal  stability  (nominally  a  packet  interval)  that 
is  deemed  necessary  for  reliable  routing  and  scheduling.  Stable  subspaces  in  the  channel 
impulse  response  have  been  identified  for  which  stable  transmission  can  be  sustained 
without  instantaneous  channel  state  information  at  the  transmitter.  This  was 
accomplished  using  channel  distribution  information  instead  of  CSI.  The  reduction  in  the 
required  feedback  requirements  for  CDI  vs  CSI  have  been  quantified.  The  effect  of 
quantized  noisy  feedback  in  the  CSI  estimates  was  also  investigated  and  the  concepts  of 
network  beamforming  and  distributed  beamforming  have  been  examined. 

A  number  of  methods  to  minimize  the  impact  of  mutual  interference  in  a  multi-user 
network  have  been  examined.  In  addition  to  the  use  of  beamsteering  to  avoid  nodes  that 
are  generating  interference  in  the  network,  the  use  of  interference  cancellation  techniques 
based  on  multi-user  detection  has  been  investigated. 
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The  overall  performanee  attainable  by  a  mobile  ad  hoe  network  depends  fundamentally 
on  the  MIMO  channel  characteristics.  During  the  past  year  the  previous  models  were 
refined  using  several  different  modeling  approaches  for  time  varying  MIMO  channels 

During  this  period  the  project  team  has  written  a  total  of  40  journal  papers  that  have  been 
published  or  accepted  for  publication.  We  have  also  presented  35  conference  papers  and 
completed  17  manuscripts  that  have  been  submitted  to  peer  reviewed  journals  and  are 
currently  under  review.  In  addition,  a  book  Modeling  MIMO  Channel  Models  was  co¬ 
authored  by  one  of  the  Pis  and  his  graduate  student  and  accepted  for  publication  by 
Wiley  Press  in  2009.  Also,  a  chapter  “Interference  Rejection  and  Management”  was  co¬ 
authored  by  three  of  the  Pis  on  the  project  and  published  in  the  book  New  Directions  in 
Wireless  Communications  Research  by  Springer- Verlag  in  2009.  The  lists  of  the  papers 
published  and  submitted  for  publication  under  this  project  is  provided  in  Appendix  3. 

A  summary  of  the  major  research  activities  conducted  and  results  that  were  obtained  is 
provided  below.  Detailed  descriptions  of  the  overall  results  are  provided  in  the  papers 
listed  in  the  Appendicees.  In  addition  these  papers  are  posted  on  the  ARO  website  and 
also  on  the  project  website  at  http://zeidler.ucsd.edu/~muri 
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2.0  Scientific  Progress  and  Accomplishments 


2.1  Antenna  Design  and  Multi-User  Feedback  Reduction  for  MIMO 
Communications  in  Mobile  Ad  Hoc  Networks 

PI:  Michael  A.  Jensen 
GSR:  Y an  Shi,  Daniel  Evans 

Prior  work  in  optimal  MIMO  antenna  design  had  not  constrained  the  solution  to  result  in 
physically  realizable  antennas.  We  extended  the  prior  work  to  incorporate  a  constraint  of 
physical  realizability  in  order  to  design  optimal  MIMO  antennas  that  are  practically 
implementable. 

Feedback  Complexity  Reduction  for  CDl-based  Signaling:  Our  joint  work  with  UCSD  has  focused  on  the 
use  of  covariance-based  transmission  for  time -variant  channels  which  reduces  the  feedback  of  the 
frequency  but  increases  the  complexity  of  each  feedback  event.  This  work  demonstrates  the  application  of 
matrix  compression  techniques  for  dramatically  reducing  the  feedback  complexity  for  CDl-based 
signaling. 


In  the  following  sections,  we  summarize  our  progress  in  these  areas. 


Optimal  MIMO  Antenna  Design 


Recent  work  has  shown  how  to  specify  antenna  radiation  patterns  which  maximize  the  diversity  gain 
given  the  stochastic  power  angular  spectmm  (PAS)  of  a  multipath  propagation  environment  and  an 
aperture  within  which  the  antennas  must  reside  [1].  However,  achieving  these  characteristics  is  difficult  if 
not  impossible,  since  they  require  that  the  current  distributions  for  the  different  elements  overlap  in  the 
aperture.  This  work  presents  several  variants  of  the  approach  which  can  determine  near-optimal  radiation 
patterns  under  the  constraint  that  the  apertures  (and  therefore  currents)  for  the  elements  do  not  overlap. 
Computational  examples  show  that  some  of  the  methods  can  consistently  provide  performance  which 
rivals  that  achievable  with  the  optimal  results  for  overlapping  currents. 


Theoretical  Framework 

We  assume  that  the  vector  incident  field  impinging  on  the  antenna  array  is  a  zero-mean 

complex  Gaussian  stochastic  process  with  the  field  arriving  at  one  angle  uncorrelated  with  that  arriving  at 
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another  angle.  The  dyadic  PAS,  or  the  average  impinging  power  per  unit  angle,  is 
P(6',  (f)  =  E{pi^^  {e,  ^  {e,  ^)} ,  where  {•}*  is  a  conjugate  and  E{-}  denotes  an  expectation.  If 

{O,  (f)  represents  the  open-circuit  radiation  pattern  of  the  mth  antenna,  the  covariance  matrix  for  the 
open-circuit  antenna  voltages  has  elements  [1] 

Cl 

The  diversity  gain  may  be  computed  by  creating  an  equivalent  system  of  uncorrelated  antennas  with  the 
branch  gains  given  by  the  eigenvalues  of  the  covariance  matrix  [2].  Therefore,  antennas  which  maximize 
these  eigenvalues  will  maximize  the  system  diversity  gain. 


The  synthesis  of  the  optimal  antenna  patterns  begins  by  using  standard  electromagnetic  analysis  to  relate 
the  mth  radiation  pattern  to  a  vector  current  function  (r  )  confined  to  an  aperture.  If  the  currents  for 

different  elements  (values  of  m)  are  allowed  to  overlap  within  the  same  aperture,  we  expand  the  current 
functions  using  orthogonal  vector  basis  functions  using 

n 

where  represents  an  unknown  weighting  coefficient.  With  this  expansion,  the  covariance  becomes 
R  =  B^CB* ,  where  is  a  transpose  and  C  can  be  computed  from  known  quantities.  Additionally, 
constraints  are  placed  on  the  coefficients  such  that  1)  each  pattern  has  unit  radiated  power  and  2)  antenna 
superdirectivity  is  limited.  These  constraints  produce  a  modified  version  of  the  matrix  C,  and  the 
dominant  eigenvectors  of  this  modified  matrix  are  used  to  construct  the  current  coefficients  which 
maximize  the  diversity  gain  [1].  Unfortunately,  the  resulting  antennas  are  difficult  if  not  impossible  to 
synthesize  since  the  specified  currents  overlap  in  the  aperture.  Our  goal  is  therefore  to  determine  non¬ 
overlapping  currents  (antennas)  whose  performance  comes  close  to  that  of  the  optimal  elements. 

Covariance  Method:  Consider  first  applying  the  eigen-analysis  to  each  aperture  individually  and  selecting 
the  dominant  eigenvector  as  the  coefficients  for  each  antenna.  This  maximizes  the  power  received  by  each 
aperture  but  neglects  cross-correlation,  and  therefore  the  approach  should  work  well  only  for  reasonably- 
spaced  apertures. 


Modified  Covariance  Method:  In  this  case,  we  determine  the  current  distribution  of  Aperture  1  using  the 
Covariance  Method  and  compute  the  diversity  gain  obtained  using  the  radiation  characteristics  of 
Aperture  1  coupled  with  the  current  distribution  associated  with  each  of  the  possible  eigenvectors  for 
Aperture  2.  The  eigenvector  producing  the  maximum  diversity  gain  is  then  selected.  For  each  subsequent 
aperture,  this  procedure  is  repeated  based  upon  the  already-selected  currents.  This  helps  to  incorporate  the 
impact  of  correlation  in  the  computation. 
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Current  Method:  Let  represent  the  mth  column  of  B  representing  the  optimal  coefficients  for  the  mth 
current  distribution,  with  corresponding  to  the  elements  of  b^  representing  the  currents  in  Aperture  i 

(this  assumes  each  basis  function  has  support  in  only  a  single  aperture).  Furthermore,  let  w'  represent  the 
unknown  coefficients  describing  the  currents  in  Aperture  i.  To  find  the  vectors  w'  that  form  a  good  basis 
for  b^,  we  solve 

w‘  =  argmax  V  cos^  6'^ 

w 

m 

where  is  the  angle  between  b^  and  w  .  The  weight  allows  placing  emphasis  on  dominant  modes 

II  i  iF 

and  is  selected  according  to  =  A^  bM  ,  where  A^  is  the  eigenvalue  of  the  covariance  matrix 
corresponding  to  the  current  weight  b^  .  Since  cos^  6^^  =  w^b^b^w/|^||b^||  w^wj,  we  can  write 

,  w^B'w 

w  =  arg  max - - - 

w  w  w 

where  B'  =  ^  optimal  value  of  w'  is  the  dominant  eigenvector  of  B' . 


Numerical  Optimization:  Starting  with  the  currents  from  the  Modified  Covariance  Method,  we  apply  a 
simplex  optimization  to  refine  the  currents  to  maximize  the  system  diversity  gain.  While  this  does  not 
guarantee  a  global  maximum,  it  does  demonstrate  the  scale  of  improvement  that  can  be  expected  from 
numerical  optimization. 

Results 

For  the  computations  used  in  evaluating  these  techniques,  we  restrict  ourselves  to  a  two-dimensional 
scenario,  where  the  basis  functions  are  two-dimensional  square  pulse  functions.  These  two-dimensional 
representations  are  suitable  for  environments  where  the  propagation  has  a  single  polarization  and  is 
confined  to  the  horizontal  plane.  The  PAS  consists  of  multiple  clusters  each  characterized  by  a  Laplacian 
function  defined  by  a  mean  angle  in  the  azimuthal  plane  computed  as  a  random  variable  uniformly 

distributed  on  [0,  360°)  ,  a  relative  strength  computed  as  a  random  variable  uniformly  distributed  on  [1, 
10],  and  an  angle  spread  (width  of  the  Laplacian).  For  each  computation,  1000  different  statistical 
realizations  of  the  PAS  are  constructed,  and  the  diversity  gain  results  represent  average  values.  The 
diversity  gains  are  computed  at  the  1%  probability  level  assuming  maximal  ratio  combining  with  an 
isolated  dipole  used  as  the  reference  [2]. 

Figure  1(a)  shows  the  geometry  of  the  aperture  arrays  used  in  the  computations.  Each  array  consists  of  a 
set  of  square  apertures,  0/2  on  a  side,  and  configured  either  in  a  linear  array  of  four  apertures  or  a  circular 
array  of  six  apertures.  In  each  aperture,  25  pulse  functions  arranged  in  a  5  by  5  grid  are  used  to  represent 
the  currents.  The  diversity  gain  obtained  using  the  radiation  characteristics  computed  using  the  methods 
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outlined  above  is  compared  with  that  obtained  from  Hertzian  dipole  antennas  placed  at  the  centers  of  the 
apertures. 

Figure  1(b)  plots  the  average  diversity  gain  for  the  linear  array  of  apertures  computed  with  the  different 
methods  as  a  function  of  the  spacing  5  between  the  aperture  edges.  The  number  of  PAS  clusters  is 
uniformly  drawn  from  a  discrete  set  of  integers  between  one  and  four,  while  the  angle  spread  of  each 

cluster  is  a  uniformly  distributed  random  variable  on  [l5° ,  50°  ]  .  These  results  show  that  as  the  spacing 
increases,  the  diversity  gain  for  the  three  approximate  methods  converges  to  the  same  value.  This  is 
because  much  of  the  diversity  is  created  by  the  wide  array  spacing,  and  therefore  the  techniques  must  only 
create  radiation  patterns  which  maximize  the  gain.  Furthermore,  the  performance  of  the  approximate 
techniques  approaches  that  of  the  optimal  radiation  characteristics  (overlapping  currents)  for  the  same 
reason.  All  of  the  techniques  produce  radiation  characteristics  which  significantly  outperform  those  of  the 
dipole  array,  emphasizing  the  fact  that  shaping  the  element  response  can  provide  significant  benefit. 
Finally,  it  is  interesting  to  recognize  that  the  performance  of  the  optimal  elements  and  the  dipoles  remain 
relatively  constant  for  spacings  larger  than  a  wavelength. 
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Figure  1 :  (a)  Geometry  of  the  apertures  used  for  the  computational  examples,  (b)  Average  diversity  gain  versus  s  for 
the  four-element  linear  array. 


Next,  we  consider  the  case  of  the  circular  array  of  apertures  and  the  corresponding  array  of  dipoles. 

Figure  2(a)  plots  the  average  diversity  gain  as  a  function  of  the  number  of  clusters  in  the  PAS  when  the 
angle  spread  is  chosen  randomly  as  outlined  above.  These  results  show  that  increasing  multipath  creates  a 
reduction  in  diversity  gain  for  all  methods,  since  the  apertures  are  no  longer  able  to  effectively  increase 
performance  by  focusing  the  radiation  pattern  in  specific  directions.  In  contrast,  the  diversity  gain  of  the 
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dipoles  increases  as  the  increased  multipath  richness  creates  a  better  environment  for  diversity  processing. 
These  curves  also  demonstrate  that  the  Modified  Covariance  Method  is  on  average  best  able  to  determine 
the  radiation  characteristics  appropriate  for  the  environment.  Figure  2(b)  shows  the  diversity  gain  versus 
the  cluster  angle  spread  when  four  clusters  are  used  in  the  PAS.  These  results  reinforce  the  observations 
concerning  multipath  richness  and  performance,  and  further  show  that  the  performance  degrades  when  the 
angle  spread  becomes  too  narrow. 


Table  I:  Diversity  Gain  Performance  Relative  to  those  of  the  Optimal  Radiation  Characteristics 


Single  PAS  (dB) 

Average  over  1000  PAS  (dB) 

Numerically  Optimized 

-0.59 

- 

Covariance 

-3.69 

-3.50 

Modified  Covariance 

-1.35 

-1.39 

Current 

-3.02 

-2.89 

Table  1  shows  the  values  of  diversity  gain  in  dB  relative  to  the  performance  of  the  optimal  overlapping 
antennas  for  a  single  representative  PAS  as  well  as  the  average  taken  over  1000  random  PAS  shapes.  The 
average  numerically  optimized  results  are  not  shown  due  to  the  large  computation  time  for  this  method. 
These  results  show  that  the  numerically-optimized  results  can  nearly  achieve  the  performance  of  the  ideal 
radiation  characteristics  and  confirm  that  the  Modified  Covariance  Method  is  superior  to  the  other 
methods. 


Number  of  Clusters 


Cluster  Angular  Spread  (degrees) 
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Figure  2:  (a)  Average  diversity  gain  versus  the  number  of  clusters  in  the  PAS  for  the  six-element  circular  array,  (b) 
Average  diversity  gain  versus  the  angle  spread  of  the  four  clusters  in  the  PAS  for  the  six -element  circular  array. 


Feedback  Reduction  for  the  MIMO  Broadcast  Channel 

The  use  of  MIMO  technology  in  the  multi-user  broadcast  channel  is  an  intriguing  possibility,  as  it  enables 
dramatic  increase  in  network  throughput  over  single-antenna  implementations.  However,  when  the 
transmit  precoding  is  performed  based  on  the  channel  state  information  (CSI)  at  the  transmitter  [3],  [4], 
the  performance  is  highly  sensitive  to  CSI  errors  due  to  inaccurate  channel  estimation,  stale  CSI  due  to 
channel  time  variations,  or  inadequate  knowledge  of  the  channel  frequency  response  in  frequency- 
selective  channels  [5]-[7].  This  has  motivated  the  development  of  a  precoding  strategy  based  on  channel 
distribution  information  (CDI)  at  the  transmitter,  in  the  form  of  spatial  correlation  matrices,  that  offers 
stable  performance  in  the  time -varying  and  frequency-selective  MIMO  broadcast  channel  [6],  [7].  While 
the  multi-user  sum  rate  of  such  an  approach  is  suboptimal  compared  to  the  case  of  non-linear  dirty  paper 
coding  (DPC)  or  linear  capacity  optimal  beamforming  techniques  with  perfect  CSI,  the  fact  that  the  CDI 
remains  more  constant  in  time  and  frequency  enables  the  approach  to  offer  superior  average  throughput 
with  infrequent  feedback  events. 


The  cost  of  using  CDl-based  techniques,  however,  is  that  the  complexity  of  each  feedback  event  can  be 
very  high.  For  example,  given  a  MIMO  system  with  transmit  and  receive  antennas,  CSI-  and 

CDl-based  schemes  respectively  require  feedback  of  an  N^.  x  channel  matrix  and  x  A^^A^, 

covariance  matrix  for  each  user.  This  dramatic  difference  in  feedback  complexity  motivates  the 
development  of  new  methods  for  reducing  the  feedback  communication  requirements.  Recently,  it  has 
been  shown  that  by  using  simple  models  -  such  as  the  well-known  Kronecker  model  -  for  the  covariance 
structure,  the  volume  of  required  feedback  data  can  be  reduced  to  a  level  similar  to  that  required  for  CSl- 
based  precoding  [6].  However,  further  research  focused  on  additional  reductions  in  the  required  feedback 
and  detailed  comparisons  showing  how  reduced-feedback  CDl-based  precoding  performs  relative  to  CSl- 
based  precoding  for  the  MIMO  broadcast  channel  has  yet  to  be  completed. 


The  goal  of  this  work  is  to  demonstrate  use  of  the  Karhunen-Loeve  (KL)  Transform  to  compress  the 
amount  of  information  that  must  be  fed  back  to  the  transmitter  when  using  CDI  signaling.  The 
performance  of  the  approach  is  assessed  using  long  channel  time  series  created  with  a  clustered  physical 
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channel  model  based  on  experimental  measurements.  The  results  show  that  excellent  results  can  be 
obtained  with  significant  reduction  in  feedback  complexity. 


System  Description  and  Channel  Model 

The  MIMO  broadcast  channel  consists  of  a  single  transmitting  node  with  N^  antennas  and  receiving 
nodes  each  with  antennas.  We  represent  the  x  channel  matrix  from  the  transmitter  to  the  jth 
user  and  at  the  kth  sample  time  as  .  Assuming  that  transmit  precoding  (non-linear  or  linear)  has  been 
applied  such  that  represents  the  x  1  precoded  vector  destined  for  the yth  user  at  the  kth  sample 
time,  the  x  1  received  vector  for  the  jth  user  at  the  kth  sample  time  can  be  expressed  as 

yf  =  +r\P 

i=U*j 

where  i\P^  is  an  A^^.  x  1  vector  of  zero-mean  unit-variance  additive  white  Gaussian  noise.  Note  that  this 

model  can  be  extended  to  frequency-selective  channels  by  adding  a  frequency  index,  although  this  is  not 
included  here  as  our  focus  is  on  narrowband  time-varying  channels. 


Transmit  Precoding 

Given  this  system  description,  we  turn  our  attention  to  the  precoding  methods  used  for  the  MIMO 
broadcast  channel.  While  dirty  paper  coding  represents  the  optimal  communication  strategy  [5],  prior 
research  has  clearly  demonstrated  that  its  performance  is  highly  sensitive  to  accuracy  of  the  transmit  CSl, 
resulting  in  rapid  throughput  degradation  as  the  CSl  becomes  stale  [6].  Since  the  behavior  of  dirty  paper 
coding  with  error  in  the  CSl  is  similar  to  that  encountered  with  linear  beamforming,  such  as  the  capacity 
optimal  regularized  channel  inversion  (CO-RCl)  [6],  we  keep  the  comparison  simple  and  use  CO-RCl  as 
the  benchmark  technique  for  CSl-based  beamforming  (CSl-BF).  This  approach  derives  the  multi-user 
beamforming  vectors  at  the  transmitter  directly  from  the  CSl  which  must  be  fed  back  from  each  receiver 
to  the  transmitter.  From  a  practical  perspective,  since  feedback  cannot  occur  for  data  transmission  at  all 
sample  times,  the  beamforming  vectors  are  derived  from  the  channel  and  are  used  for  some  finite 

time  window  characterized  by  a  finite  set  of  samples  Uq  <  k  <  +  K  . 


For  CDl-based  beamforming  (CDl-BF),  we  estimate  the  covariance  of  the  channel  as 
Ry  =  },  where  is  a  vector  formed  from  column-wise  stacking  of  ,  and  the 

expectation  is  performed  over  time  (with  index  k).  We  use  the  algorithm  developed  in  [6]  to  construct 
beamforming  vectors  from  the  matrices  R  ^ ,  1  <  j  <  .  As  demonstrated  by  the  results  obtained  with 
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measured  channel  data  assuming  the  receiving  nodes  are  mobile  [6],  we  can  assume  that  this  covariance 
matrix  remains  relatively  constant  over  some  reasonable  time  window. 


As  we  consider  the  construction  of  the  beamforming  vectors  for  these  two  signaling  methods,  we 
immediately  see  that  CSl-BF  and  CDl-BF  respectively  require  feedback  of  matrices  of  size  x 

and  x  N^N^.  While  CDl-BF  requires  less  frequent  feedback,  to  be  competetive  either  the  feedback 

frequency  must  be  reduced  by  a  factor  of  or  the  amount  of  data  must  be  dramatically  reduced.  For 

systems  where  ~  A^t  >  ^^is  can  be  effectively  accomplished  by  parameterizing  the  covariance  matrices 

using  a  physical  model  such  as  the  Kronecker  model.  However,  for  situations  where  is  small,  use  of 
CDl-BF  with  these  parameterizations  still  requires  feedback  of  larger  matrices  than  would  be  required  for 
CSl-BF  in  the  same  situation.  This  observation  motivates  the  work  presented  in  this  paper. 


In  comparing  these  beamforming  methods,  we  quantify  their  performance  using  the  sum  rate  of  the  multi¬ 
user  channel.  Specifically,  let  represent  the  sum  mutual  information  achieved  by  a  specific 

precoder  x  >  where  x  selected  from  the  set  {CSl-BF,  CDl-BF},  at  the  time  index  k  when  the 
precoding  is  constructed  from  CSl  or  CDl  at  time  index  ko  <  k  .  The  sample  expected  throughput  (SET) 

is  then  defined  by  averaging  this  mutual  information  over  different  starting  times  for  the  same  time  offset 
using 


K-A 


"  X 


k  m=\ 


where,  as  a  reminder,  K  is  the  total  number  of  time  samples.  We  also  use  the  average  SET  (ASET) 
defined  by 


1  M-l 


A=0 


where  M  represents  the  time  displacement. 


The  CO-RCl  algorithm  considered  in  this  paper  assumes  a  single  data  stream  per  user,  such  that  the 
symbol  x^P  destined  for  the  jth  user  is  transmitted  using  the  beamforming  weights  constructed 
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from  the  CSl  Beamforming  weights  constructed  from  the  CSl  are  applied  at  the 

receiver.  The  signal-to-interference-plus-noise  ratio  (SINK)  for  the  jth  user  under  this  scenario  and 
assuming  phase  synchronization  at  the  receiver  becomes  [6] 
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It  is  possible  that  the  beamforming  weight  determination  excludes  a  user  from  the  channel  [8],  in  which 
case  such  that  p ^  =  0  .  The  sum  mutual  information  obtained  by  summing  the  mutual  information 

values  for  the  users  is  given  by 


/^(ko,k)  _ 

'^CSI-BF  ~ 


Jl0g2(l  +  P,). 


The  goal  of  the  CDl-based  beamforming  approach  considered  in  this  work  is  to  maximize  the  throughput 
averaged  over  a  specified  time  window  [6],  Once  the  beamforming  weights  are  constructed  from  the 
algorithm,  the  sum  mutual  information  constructed  by  using  the  appropriate  weights  in 

the  formulation  above.  The  results  in  [6]  demonstrate  that  these  beamforming  weights  provide  a  stable 
throughput  in  a  time -varying  channel. 


Channel  Model 

Evaluations  of  the  CSl-BF  and  CDl-BF  algorithms  considered  here  using  time-varying  channels 
measured  in  several  different  environments  reveal  that  while  the  CSl-BF  technique  has  high  initial  SET 
(when  the  CSl  is  of  high  quality),  its  performance  degrades  rapidly  with  time  offset  .  In  contrast,  the 

CDl-BF  technique  provides  stable  sum  rate  performance  over  a  relatively  long  time  period.  The  problem 
with  this  analysis,  however,  is  that  the  temporal  extent  of  each  measured  channel  time  series  available  for 
this  study  is  inadequate  for  determining  the  temporal  extent  over  which  the  CDl-BF  vectors  remain  valid. 


As  a  result  of  this  limitation  in  the  measured  data,  in  this  paper  we  resort  to  the  use  of  a  physics-based 
channel  model  [9].  This  recently  developed  model  extracts  the  dominant  multipath  clusters  from 
measured  data  and  uses  an  auto-regressive  model  to  simulate  the  time  variation  in  the  key  cluster 
parameters  (angle  of  departure  or  arrival  and  channel  gain).  The  model  includes  a  mechanism  for 
allowing  new  clusters  to  appear  and  existing  clusters  to  disappear  based  on  the  probability  distributions 
observed  in  real  measurements.  Analysis  of  the  channel  matrices  obtained  with  this  model  show  that  they 
exhibit  behaviors  very  consistent  with  the  behaviors  observed  in  the  measured  data.  We  therefore  use  this 
model  in  our  analyses.  In  the  implementation,  each  multipath  cluster  consists  of  50  components  or  rays. 
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with  the  angle  of  departure  (or  arrival)  for  each  ray  relative  to  the  cluster  mean  angle  specified  as  a 
random  variable  satisfying  a  truncated  Laplacian  distribution  with  an  angle  spread  of  30° .  The  model 
implementation  assumes  that  each  time  sample  represents  a  physical  spacing  of  0.0086  wavelengths, 
which  matches  the  sampling  of  the  underlying  measured  data. 


It  is  important  to  also  emphasize  that,  as  in  the  analyses  performed  previously,  channel  variation  is 
created  in  this  study  as  the  receiving  nodes  move.  Therefore,  while  this  is  manifest  as  a  time -variation  in 
the  channel,  the  key  physical  phenomenon  that  creates  the  variation  is  node  spatial  displacement. 
Therefore,  the  results  shown  in  this  paper  are  all  given  as  a  function  of  node  displacement  in  wavelengths 
at  the  center  frequency  of  the  simulation.  Given  a  velocity  for  the  nodes,  it  is  a  simple  matter  to  convert 
the  displacement  independent  variable  to  time. 


KL  Transform  for  Feedback  Reduction 


Methodology 

We  emphasize  at  this  point  that  CDl-BF  is  designed  to  maximize  the  average  sum  rate  of  the  broadcast 
channel  communication,  which  implies  that  each  beamforming  vector  provides  good  average  signal 
strength  to  the  intended  user  while  controlling  the  average  interference  to  the  other  users  in  the  network. 

If  we  define  the  effective  null  space  of  as  the  subspace  characterized  by  eigenvalues  smaller  than 

some  threshold  (since  in  most  multipath  channels  a  true  null  space  of  the  covariance  will  not  exist),  then 
good  CDl-BF  performance  implies  we  can  find  a  beamforming  vector  for  user  j  which  lies  in  the  effective 
range  space  of  and  the  effective  null  space  of  R; ,  i  ^  j .  Practically  speaking,  if  this  is  to  work  for  a 

reasonable  number  of  users,  it  implies  that  each  covariance  matrix  is  characterized  by  a  few  dominant 
eigenvectors. 


It  is  therefore  intuitive  to  try  to  use  this  property  to  reduce  the  feedback  complexity,  which  naturally 

motivates  an  evaluation  the  KL  Transform  for  feedback  data  compression.  Specifically,  let 

R^  =11 J  A  represent  the  eigenvector  decomposition  of  R^ ,  where  because  R^  is  Hermitian  we 

know  that  H  j  is  the  unitary  matrix  of  eigenvectors  and  A  j  is  the  diagonal  matrix  of  real,  non-negative 

eigenvalues  constructed  so  that  the  eigenvalues  fall  in  order  of  decreasing  value.  The  feedback  reduction 
comes  by  feeding  back  only  the  dominant  eigenvectors  and  the  corresponding  eigenvalues  for  each  user. 


SET  V5.  Node  Displacement 
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Consider  the  case  where  =  5  and  the  total  allocated  transmit  power  for  the  beamforming 

methods  is  set  at  A*  =  10  [6].  Figure  3(a)  plots  the  SET  as  a  function  of  the  receiver  displacement 
assuming  that  only  one  or  all  25  eigenvectors  and  eigenvalues  are  fed  back  to  the  transmitter  at  the  initial 
location  (displacement  =  0).  The  values  represent  averages  taken  over  50  random  channel  realizations 
from  the  stochastic  model.  The  vertical  line  in  the  plot  shows  the  time  window  over  which  the  covariance 
is  estimated,  which  means  that  the  portion  of  the  curve  representing  true  communication  throughput  for 
CDl-BF  is  that  to  the  right  of  this  solid  line.  It  is  noteworthy  that  feedback  of  a  single  eigenvector 
provides  performance  close  to  that  obtained  using  full  covariance  feedback  (25  eigenvectors).  The  SET 
for  the  case  of  CSl-BF  assuming  the  CSl  at  the  receiver  is  known  perfectly  at  all  times  but  the  CSl  at  the 
transmitter  is  fed  back  only  at  the  initial  location  is  also  shown  for  comparison.  In  terms  of  feedback 
complexity,  CSl-BF  requires  feedback  of  25  complex  numbers,  which  is  the  same  as  that  required  for 
CDl-BF  with  a  single  eigenvector.  The  difference,  however,  is  that  the  CDl  remains  useful  for  a  longer 
duration,  demonstrating  the  effectiveness  of  the  technique. 


Figure  3(b)  repeats  the  plot  of  Fig.  3(a)  for  the  case  where  the  covariance  matrix  fed  back  is  modeled 
using  a  Kronecker  product  approximation  with  the  one-sided  covariance  matrices  computed  using  the 
Rank-1  approximation  [6].  In  this  case,  each  one-sided  covariance  is  5x5,  and  feedback  of  a  single 
eigenvector  for  each  of  the  two  one-sided  covariance  matrices  results  in  feedback  of  two  vectors  each 
with  five  complex  numbers  (along  with  two  real  eigenvalues).  The  performance  in  this  case  is  similar  to 
that  in  Fig.  3(a)  despite  the  additional  reduction  in  feedback  complexity. 


Achievable  System  Throughput 

We  can  further  manipulate  these  metrics  to  represent  achievable  performance  in  a  practical 
implementation.  Conceptually,  if  a  receiver  possesses  either  CSl  or  CDl  at  a  point  in  time  and  must  feed 
that  information  back  to  the  transmitter,  and  if  the  receiver  (or  transmitter)  is  mobile,  then  by  the  time  the 
receiver  has  fed  back  the  information  to  the  transmitter,  that  information  is  stale,  resulting  in  performance 
degradation  relative  to  the  case  of  perfect  CSl  or  CDl.  Based  on  the  fact  that  degradation  with 
displacement  is  less  for  CDl-BF  than  for  CSl-BF  (as  shown  in  Fig.  3(a)),  we  expect  that  a  system  using 
CDl  will  suffer  less  relative  degradation,  although  the  actual  throughput  performance  depends  on  the 
capabilities  of  the  feedback  channel. 
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(a)  (b) 

Figure  3 :  SET  Curves  for  a  5  x  5  MIMO  network  with  =  5  and  P  =  10  with  different  numbers  of  covariance 

matrix  eigenvectors  fed  back  to  the  transmitter.  The  KL  Transform  of  the  (a)  full  covariance  matrix  or  (b)  rank-1 
approximation  to  the  Kronecker  product  of  the  covariance  matrix  is  used. 


To  get  a  more  accurate  understanding  of  the  throughput  behavior,  we  assume  that  at  zero  displacement, 
the  receivers  have  estimates  of  either  their  CSI  or  their  CDl  as  appropriate  and  that  the  nodes  are  capable 
of  full  duplex  communication.  Each  receiving  node  quantizes  the  required  feedback  information  into  B 
bits  assuming  64  bits  are  used  for  each  real  number  and  subsequently  feeds  the  information  back  over  the 
channel  with  a  bit  rate  (which  incorporates  feedback  channel  bandwidth,  modulation,  and  coding). 

With  each  receiver  node  moving  at  a  velocity  v,  the  maximum  Doppler  frequency  is  =  2vl  A  , 

where  □  is  the  wavelength  at  the  communication  center  frequency.  We  refer  to  the  distance  moved  (in 
wavelengths)  by  the  nodes  over  the  time  during  which  the  feedback  is  taking  place  as  the  latency 
(consistent  with  the  horizontal  axis  in  Fig.  3(a)),  which  can  be  computed  using 


fD.,n 


2R. 


B. 


We  then  compute  the  ASET  for  the  communication  assuming  feedback  occurs  at  an  interval  (in 
wavelengths)  of  A  ^  ,  where  the  window  over  which  the  covariance  information  is  estimated  is  that  shown 
in  Fig.  3(a)  (which  corresponds  to  200  samples  at  0.0086  wavelengths  per  sample). 


In  the  following  computations,  we  assume  the  receiving  nodes  move  at  vehicular  speeds  of  72  km/hour 
and  communication  occurs  at  a  center  frequency  of  2.4  GHz,  resulting  in  =  3.2  Hz  .  ASET  is 

computed  over  a  window  of  1000  samples  (node  displacement  of  8.6  wavelengths),  where  the  number  of 
feedback  events  within  that  window  is  determined  by  the  feedback  interval  A^  .  For  CDl-BF,  feedback 
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uses  a  single  eigenvector  for  each  of  the  one-sided  covariance  matrices  computed  using  the  Rank-1 
approximation  with  the  Kronecker  model.  The  system  again  uses  =  5  . 


Figure  4(a)  plots  the  ASET  for  this  system  as  a  function  of  the  feedback  rate  for  two  different  values 

of  the  feedback  interval  for  both  CSl-BF  and  CDl-BF.  These  results  show  that  the  performance  of 

both  beamforming  strategies  increases  with  the  feedback  rate  because  of  the  reduced  latency  between  the 
time  of  estimation  and  time  of  use  of  the  CSI  or  CDI.  The  impact  of  rate  on  CSI-BF  is  particularly 
significant,  which  can  be  explained  by  the  rapid  degradation  in  performance  with  receiver  motion  shown 
in  Fig.  3(a)  for  CSl-BF.  In  other  words,  for  CSl-BF,  increased  latency  in  the  feedback  results  in  less 
accurate  CSI  at  the  outset  of  its  use,  and  therefore  the  system  cannot  achieve  the  high  performance 
obtainable  for  small  displacements  in  Fig.  3(a).  In  contrast,  CDl-BF  is  more  forgiving  of  stale  CDI, 
resulting  in  reduced  sensitivity  to  feedback  rate  in  Fig.  4(a).  More  importantly,  these  results  show  that 
unless  the  feedback  rate  is  high  and  the  feedback  interval  is  small,  CDI-BF  outperforms  CSI-BF  in  terms 
of  average  system  throughput. 


(a)  (b) 

Figure  4:  Effective  system  throughput  (ASET)  of  CSI-BF  and  CDI-BF  (a)  versus  the  feedback  rate  for  two 
different  values  of  the  feedback  interval  and  (b)  versus  the  feedback  interval  for  two  different  values  of  the 
feedback  rate  r^. 


These  observations  are  reinforced  by  the  results  in  Fig.  4(b)  which  plots  the  ASET  for  the  system  as  a 
function  of  the  feedback  interval  A  ^  for  two  different  values  of  the  feedback  rate  .  Once  again,  we  see 

that  CSI-BF  outperforms  CDI-BF  only  when  the  feedback  interval  is  small  and  the  feedback  rate  is  high. 
These  results  further  highlight  that  the  improvement  offered  by  CDl-BF  is  particularly  enhanced  when  the 
feedback  rate  is  low,  since  the  increased  latency  significantly  degrades  the  performance  of  the  CSI-BF 
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approach.  Finally,  the  curves  emphasize  that  CSl-BF  experiences  increased  sensitivity  to  the  feedback 
interval. 
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2.2  Cross-Layer  Design  and  Analysis  of  MAC  and  Routing  Protocols 
for  Ad  Hoc  Networks  with  Multiple  Antennas 


PI:  Prof.  Michele  Zorzi 


Introduction 


During  the  past  year,  Dr.  Zorzi's  research  efforts  have  been  focused  on  different  issues  concerning  the 
PHY.  MAC  and  routing  protocol  layers.  Novel  results,  as  well  as  wrap-up  research  completing  the 
activities  carried  out  during  the  last  year,  include: 


1)  The  proposal  of  a  novel  practical  system  for  decode-and-forward  physical  layer  network  coding,  a 
topic  to  which  increasing  effort  has  been  devoted  by  the  research  community  in  the  past  four 
years.  While  remarkable  improvements  have  been  carried  out  from  a  theoretical  point  of  view,  the 
definition  of  practical  architectures  have  not  yet  received  as  much  attention.  In  this  field,  we  have 
virtually  solved  all  issues  that  would  hinder  implementation  by  mixing  some  well-known 
techniques  (such  as  OFDM)  and  new  ideas  (network  demodulation)  in  the  right  context 

2)  The  routing  problem  in  MIMO  ad  hoc  networks  has  been  addressed  again  in  relation  to  MAC 
issues  and  PHY  issues  such  as  imperfect  channel  estimation.  The  main  objective  here  is  to  wrap 
up  the  results  of  the  previous  years  concerning  packet  forwarding  strategies  for  MIMO  ad  hoc 
networks  yet  provide  new  theoretical  results  explaining  some  behaviors  and  projecting  the 
evolution  of  these  behaviors  as  such  usual  assumptions  as  perfect  channel  estimation  are 
removed;  checking  whether  previous  results  hold  or  not  is  key  to  a  correct  protocol  design. 

3)  The  cooperative  paradigm  for  ad  hoc  networking,  which  exploits  diversity.  Hybrid  ARQ,  and 
MIMO  techniques,  has  been  explored  further.  Routing  protocols  have  been  designed  that  take 
explicit  advantage  of  these  advanced  PHY  techniques  in  a  multi-hop  network  setting.  We  have 
shown  that  this  approach  has  the  potential  to  significantly  improve  the  networking  performance, 
and  is  deeply  cross-layer  as  it  couples  routing  choices  and  opportunistic  networking  behaviors 
with  PHY  and  MAC  mechanisms. 


In  this  document,  we  provide  an  overview  of  the  main  technical  issues  and  contributions,  as  well  as  some 
sample  performance  results,  in  these  three  areas,  referring  to  the  specific  papers  for  a  detailed  description 
and  for  an  extensive  set  of  results. 
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MIMO,  Network  Coding,  and  PHY  integration 


During  the  year  2008,  Prof.  Zorzi’s  efforts  on  network  coding  have  focused  on  designing  new  techniques 
that  better  integrate  this  technique  with  the  physical  layer.  In  the  last  months  new  work  has  been  carried 
out  in  this  area,  especially  around  physical  layer  network  coding  (PNC).  This  concept  has  lately  attracted 
a  certain  interest  in  the  network  coding  community  [l]-[7].  The  basic  idea  can  be  illustrated  in  the  two- 
way  relay  channel  (TWRC,  Fig.  1).  Two  nodes  (A  and  B)  must  exchange  two  packets  (A's  X  and  B's  Y) 
through  an  intermediate  relay  R.  In  classic  Network  Coding  (NC),  A  would  send  X  in  time  slot  0,  B 
would  transmit  Y  in  slot  1  and  R  would  broadcast  X©Y  in  slot  2  (Fig.  1  a).  In  physical  layer  network 
coding,  A  and  B  would  simultaneously  send  X  and  Y,  while  R  would  relay  a  function  of  X  and  Y,  which 
is  invertible  in  X  or  Y  as  soon  as  the  other  variable  is  known  (Fig.  lb).  Given  that  A  and  B  know  their 
own  packets,  they  can  each  potentially  decode  the  other  node's  frame.  This  method  reduces  the  number  of 
required  slots  from  3  to  2.  Note  that  R  need  not  decode  X  and  Y  separately,  but  it  is  enough  to  directly 
decode  a  linear  combination  Z.  This  can  potentially  reduce  the  error  rate  because  less  information  needs 
to  be  extracted  from  the  received  signal  [1],  [3],  [4]. 


At  least  two  main  ideas  have  been  developed  in  this  context.  The  first  one  is  called  amplify-and-forward 
physical  layer  network  coding  (AF-PNC,  also  called  analog  network  coding  [3],  [6]).  The  relay  amplifies 
the  analog  superposition  of  X  and  Y  and  broadcasts  this  signal,  that  will  be  called  S.  The  end-user 
terminals  A  and  B  can  recover  the  signal  sent  by  the  other  node  by  subtracting  their  own  waveform  out  of 
the  received  signal  sent  by  R.  The  resulting  signal  depends  only  on  the  intended  packet,  under  the 
hypothesis  of  perfect  cancellation.  On  the  other  hand,  in  decode -and-forward  physical  layer  network 
coding  (DF-PNC  [3], [4]),  R  decodes  a  linear  combination  L  of  X  and  Y  directly  from  the  analog 
superposition.  This  packet  L  is  broadcast  and  A  and  B  can  remove  their  own  frame  from  L  to  recover  the 
desired  data  unit. 


Such  systems  are  of  interest  in  the  context  of  ad  hoc  MIMO  wireless  networks  for  a  variety  of  reasons. 
First  of  all,  the  presence  of  multiple  antennas  at  the  terminals  is  not  precluded,  although  most  of  the 
current  studies  focus  on  single  antenna  terminals.  We  note  that  even  in  the  presence  of  single-antenna 
terminals,  the  processing  at  the  receiver  is  inherently  MIMO,  where  the  multiple  input/output 
relationships  are  generated  by  the  network  coding  technique  rather  than  by  multi-element  antennas. 
Actually,  the  introduction  of  multiple  antennas  at  R  can  dramatically  improve  the  performance  of  DF- 
PNC:  DF-PNC  suppresses  the  noise  at  every  decoding  step.  Hence,  it  is  able  to  work  at  higher  SNR  than 
AF-PNC.  Clearly,  the  higher  the  diversity,  the  smaller  the  bit  error  rate.  Hence,  the  introduction  of 
multiple  antennas  at  least  at  the  relay  can  significantly  enhance  the  performance  of  DF-PNC  over  AF- 
PNC.  It  is  thus  of  great  interest  to  explore  practical  architectures  for  DF-PNC,  since  they  would  benefit 
more  than  other  competing  systems  from  the  properties  of  MIMO  techniques.  Prof  Zorzi’s  efforts  have 
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focused  on  SISO  nodes  as  a  necessary  and  preliminary  step  towards  MIMO  nodes,  but  studies  of  MIMO 
centric  DF-PNC  are  a  natural  extension. 


Furthermore,  physical  layer  network  coding  is  especially  suited  for  bidirectional  traffic.  Such  type  of 
traffic  is  of  special  importance  in  military  or  emergency  ad  hoc  networks,  where  VoIP  or  video  sessions 
are  necessary  among  the  participants.  Moreover,  NC  is  well  known  to  be  beneficial  also  for  multicast 
packets,  and  these  frames  are  again  very  commonplace  for  wireless  ad  hoc  networks. 


Figure  1:  Comparison  of  Digital  and  Physical  Layer  Network  Coding 
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While  AF-PNC  has  been  implemented  in  a  real-world  testbed  [6],  no  actual  implementation  of  DF-PNC 
has  yet  been  built  (or  at  least  disclosed).  Some  of  the  biggest  problems  in  this  field  are  the  necessity  of 
having  symbol  synchronous  signals  at  R  and  flat  fading  channels.  Such  a  symbol  synchronization  induces 
a  certain  structure  in  the  received  signals  that  eases  the  decoding  process,  for  instance  enabling  the  usage 
of  sophisticated  lattice-based  schemes  [4].  In  addition,  the  demodulation  is  rather  different  from  the 
conventional  one  and  nobody  has  yet  defined  how  this  process  would  yield  soft  values  for  the  channel 
decoder.  So  far,  only  hard  detection  has  been  designed  and  this  can  be  a  relevant  hindrance  for  modern 
channel  codes. 


Prof.  Zorzi's  effort  in  this  area  has  aimed  to  solve  the  aforementioned  practical  issues  for  DF-PNC 
(namely,  design  a  DF-PNC  system  that  makes  no  assumptions  on  symbol  synchronization  and  multipath 
fading  and  can  work  with  modern  channel  codes).  The  main  result  has  been  the  proposal  of  a  transceiver 
architecture  (depicted  in  Fig.  2).  The  main  elements  of  this  system  are  the  following,  while  the  details  can 
be  found  in  [7].  First  of  all,  data  is  channel  encoded  by  a  potentially  non  binary  channel  code.  The 
resulting  Galois  symbols  are  mapped  into  a  set  of  OFDM-modulated  QAM  symbols  which  are  pulse 
shaped  and  sent  into  the  wireless  medium.  The  receiver  will  collect  the  superposition  of  potentially 
different  quasi-simultaneous  waveforms  (i.e.,  roughly  frame  synchronized  but  not  necessarily  symbol 
synchronized)  from  multiple  transmitters.  The  receiver  will  first  estimate  the  channel,  perform  OFDM 
demodulation,  map  the  resulting  samples  into  soft  values  by  means  of  so  called  “network  demodulation” 
and  finally  decode  these  values  into  information  bits. 
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Figure  2:  System  block  architecture 


The  network  demodulation  block  is  a  signal  processing  algorithm  that  takes  as  an  input  symbols  that  are 
symbol  synchronous  (thanks  to  OFDM)  and  yields  the  soft  values  of  the  Galois  symbols  for  an  arbitrary 
linear  combination  of  the  collided  packets.  The  network  demodulator  works  on  a  symbol-by-symbol  basis 
and  lists  all  most  likely  sets  of  modulated  symbols  sent  by  the  transmitters.  The  key  property  is  that 
different  sets  of  transmitted  symbols  lead  to  the  same  linear  combination.  Hence,  the  probability  that  one 
element  of  the  resulting  linear  combination  assumed  a  certain  value  can  be  computed  by  marginalization, 
i.e.,  by  summing  up  the  probabilities  of  all  pairs  of  transmitted  symbols  that  yield  the  same  output.  This 
enumeration  and  summation  can  be  performed  efficiently  by  sphere  decoding. 


Some  brief  remarks  are  in  order  so  as  to  fully  appreciate  the  benefits  of  this  architecture: 


Usage  of  nonbinary  codes:  Network  Coding  requires  a  potentially  large  Galois  field  to 
properly  work  in  a  generic  network.  In  order  to  do  it,  it  is  necessary  to  work  with  nonbinary 
channel  codes. 

OFDM:  OFDM  is  a  necessary  block  in  our  architecture,  since  it  allows  to  compensate  for  a 
relative  delay  between  the  users  up  to  the  cyclic  prefix.  This  is  in  stark  contrast  with  single 
carrier  DF-PNC,  where  the  relative  delay  can  be  up  to  a  small  fraction  of  the  symbol  duration; 
hence  the  system  robustness  has  increased  by  2-3  orders  of  magnitude.  In  addition,  frequency 
selective  channels  are  naturally  handled  by  OFDM  and  no  special  provision  must  be  taken  for 
it. 

Network  demodulation:  network  demodulation  is  the  core  and  novel  element  of  our 
architecture.  It  yields  an  arbitrary  linear  combination  of  the  colliding  packets  by  listing  all 
most  likely  sets  of  received  symbols.  In  addition,  it  is  flexible  enough  to  compute  the  soft 
values  of  the  decoded  symbols,  hence  enabling  the  usage  of  modern  channel  decoders  (such 
as  LDPC  or  turbo  codes). 


The  performance  of  AF-PNC  and  our  proposed  scheme  have  been  compared  in  the  TWRC  for  Rayleigh 
fading  channel.  As  shown  in  Fig.  3,  the  proposed  scheme  (blue  line  called  “Heuristic  Digital”)  performs 
very  close  to  the  synchronous  case  (i.e.,  the  signals  arrive  simultaneously  to  the  relay).  In  addition,  a  gain 
of  about  2.5  dB  is  possible  with  respect  to  AF-PNC. 

The  simulation  results  have  shown  that  even  a  suboptimal  implementation  of  our  DF-PNC  system 
concept  is  able  to  outperform  AF-PNC.  In  addition  to  the  performance  improvements,  other  advantages 
are  brought  by  this  technique.  For  example,  we  expect  DF-PNC  to  be  more  robust  to  channel  estimation 
errors,  since  in  AF-PNC  the  interference  cancellation  process  is  subject  to  these  errors;  instead  DF-PNC  is 
not  affected,  because  A  and  B  subtract  their  own  digital  packet  (X  and  Y,  respectively),  rather  than  the 
estimate  of  the  amplified  analog  signal  at  R.  Moreover,  AF-PNC  needs  oversampling  to  yield  satisfactory 
performance,  whereas  our  architecture  still  achieves  good  results  without  it.  Furthermore,  AF-PNC  can 
transmit  only  a  signal  proportional  to  the  sum  of  the  colliding  packets.  Instead,  DF-PNC  can  potentially 
yield  any  linear  combination,  and  this  flexibility  is  beneficial  for  NC  in  realistic  topologies.  Finally,  AF- 
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PNC  does  not  work  in  frequency  selective  channels,  while  our  proposal  naturally  accommodates  these 
environments  because  of  OFDM.  All  these  points  encourage  additional  investigation  on  this  technique. 


0 
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MIMO  MANETs  with  Decision-Feedback  Multiuser  Detection 


1 )  Integrated  Multiuser  detection  /  MAC  protocol  /  routing 


Dr.  Zorzi’s  main  line  of  efforts  has  been  directed  toward  integrating  PHY,  MAC  and,  ultimately,  routing 
design.  Designing  a  MAC  layer  helps  understand  the  interaction  between  network  protocols  and  PHY, 
giving  insight  on  how  the  spatial  de -multiplexing  algorithm  should  be  “driven”  in  order  to  achieve 
transmission  effectiveness.  The  main  problem  to  face  here  is  to  answer  the  following  questions:  how 
should  the  nodes  decide  what  to  transmit  and  to  which  receiver?  How  much  should  transmission 
parallelism  be  encouraged?  What  is  the  impact  of  the  underlying  policies  and  protocols  on  routing?  The 
studies  carried  out  during  the  previous  years  led  to  the  design  of  a  MAC  protocol  that  tries  to  address 
these  problems.  This  year,  this  protocol  has  been  the  subject  of  further  investigations  aimed  at  assessing 
its  actual  capabilities  and  at  substantiating  the  assumptions  behind  it.  This  has  led  to  a  journal  paper  that 
has  been  submitted  to  the  IEEE  Transactions  on  Wireless  Communications  and  is  currently  being 
reviewed  for  publication. 


In  order  to  fix  ideas,  let  us  describe  shortly  the  basics  of  the  MAC  protocol  and  the  routing  protocol 
working  on  top  of  it.  Our  MAC  exploits  multiuser  detection  as  offered  by  advanced  MIMO -BLAST 
signal  processing  to  support  multiple  access  in  a  different  way  than  in  standard  CSMA/CA  networks. 
Namely,  the  concept  of  collision  is  softer,  and  nodes  are  encouraged  to  access  the  channel  simultaneously, 
leaving  it  to  signal  processing  to  decouple  superimposed  signals.  This  access  paradigm  bears  the  potential 
to  boost  network  performance;  however,  it  requires  some  distributed  control  on  the  nodes’  part  so  that 
simultaneous  transmissions  do  not  exceed  the  decoupling  capabilities  of  PHY -level  signal  processing.  To 
reach  this  objective,  the  MAC  protocol  works  in  frames.  Each  frame  is  composed  of  4  slots,  which  are 
used  to  transmit  Request-To-Send  (RTS),  Clear-To-Send  (CTS),  data  and  ACKnowledgment  packets, 
respectively.  RTSs  and  CTSs  do  not  carry  the  same  meaning  as  in  standard  CSMA/CA:  they  are  not 
employed  to  stop  neighbors  from  transmitting,  but  only  to  advertise  transmit  intentions,  inform  receivers, 
and  get  back  instructions  as  to  controllable  transmission  parameters  (namely,  the  number  of  antennas  that 
should  be  used  to  transmit,  which  directly  affects  the  achievable  bit  rate  as  well  as  the  receiver  load). 
However,  traffic  control  does  not  take  place  only  at  the  receivers:  transmitters  control  the  injected  traffic 
by  limiting  both  the  number  of  receivers  they  call  upon  in  the  CTSs  and  the  overall  number  of  antennas 
used.  In  particular,  nodes  are  divided  into  classes;  the  class  of  a  node  conveys  information  about  the 
maximum  number  of  antennas  to  be  used  when  transmitting  to  a  set  of  neighbors  which  includes  that 
node.  RTSs  are  formed  so  as  never  to  exceed  a  neighbor’s  class  [8].  We  usually  indicate  classes  within 
curly  brackets:  for  example,  {2,  4,  6,  8}  means  that  a  node  can  only  take  a  class  value  equal  to  one  of  the 
numbers  in  the  brackets.  Upon  receiving  (possibly  many)  RTSs  the  receivers  decide  whether  to  grant 
transmissions  or  not  by  applying  a  so-called  CTS  policy.  The  policy  takes  into  account  how  many  packets 
represent  wanted  traffic  and  how  many  represent  interference,  decides  whether  some  of  this  interference 
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should  be  detected  and  canceled  to  improve  the  reception  of  wanted  packets,  and  correspondingly 
compiles  instructions  in  CTS  packets  according  to  a  balance  between  these  two  factors.  This  is  the 
concept  behind  the  Follow  Traffic  (FT)  policy  [8]. 

After  CTS  reception,  data  transmissions  follow,  using  as  many  antennas  as  indicated  in  the  CTS,  and 
sending  one  Packet  Data  Unit  (PDU)  of  fixed  length,  through  each  antenna.  Finally,  ACKs  are  sent  back 
by  receivers  to  provide  feedback.  Transmitters  receiving  no  CTS  in  response  to  a  sent  RTS  back  off  for  a 
random  amount  of  time,  chosen  within  an  exponentially-increasing  window.  Two  backoff  policies  are 
considered,  Dest-Lock  (that  blocks  attempts  toward  the  node  that  denied  the  transmission)  and  Node- 
Lock  (that  blocks  transmissions  toward  all  nodes,  and  thus  conservatively  blocks  excess  traffic).  Nodes 
need  not  be  strictly  symbol-synchronous  to  work  in  accordance  with  this  MAC  protocol:  only  frame- 
synchronization  is  required,  as  the  actual  performance  decrease  due  to  the  absence  of  symbol 
synchronization  can  be  obviated  by  oversampling  the  receive  antenna  inputs  and  running  the  de¬ 
multiplexing  algorithm  on  a  larger  set  of  samples.  This  turns  the  synchronization  problem  into  a  receiver 
hardware  complexity  problem. 


MAC  and  routing  are  tightly  related  in  MIMO  ad  hoc  networks.  Take  for  example  our  protocol  above: 
shortest  path  routing  is  not  necessarily  the  best  choice.  In  fact,  short  paths  tend  to  rely  on  longer  links;  in 
a  MIMO  network,  long-reach  links  may  require  transmissions  with  few  antennas  in  order  to  keep  the 
transmit  power  high  and  transmit  signals  effectively.  Therefore  a  shortest  path  is  also  a  lower  bit  rate  end- 
to-end  connection.  On  the  other  hand,  adaptive  routing  policies  could  show  improved  communication 
performance  by  choosing  locally  optimal  relays  over  a  multihop  path  in  a  dynamic  fashion.  The 
comparison  below  shows  the  apparent  benefits  yielded  by  dynamic  relay  selection.  On  the  left  (Figure  4), 
we  considered  a  typical  Poisson  traffic  model,  whereby  packets  are  generated  randomly  at  given  rate  by 
all  nodes:  in  this  very  simple  case  dynamic  selection  outperforms  even  a  pre-computed  set  of  shortest 
paths,  obtained  through  the  Djikstra  algorithm  and  optimized  for  the  minimum  number  of  hops.  On  the 
right  (Figure  5),  we  have  set  up  a  more  complex  scenario  involving  a  more  realistic  type  of  traffic,  more 
bursty  in  nature  and  whose  generation  is  correlated  in  time.  According  to  the  model  in  [10],  traffic  is 
decomposed  in  sessions,  each  originating  flows;  each  flow  is  composed  of  a  random  number  of  packets. 
The  number  of  flows  within  a  session  is  also  random.  In  time,  packets  and  flows  are  generated  according 
to  a  Poisson  process  of  given  rate.  In  the  correlated  traffic  case,  the  dynamic  relay  selection  policy 
outperforms  shortest  path  routing  more  than  with  uniform  traffic,  and  increases  throughput  by  about  25%. 
With  either  traffic  type,  dynamic  relay  selection  accounts  for  many  node  parameters,  such  as  current 
queue  backlog  level,  advancement  provided,  class  of  the  neighbor  and  backoff  timer.  A  specifically 
designed  weighing  function  yields  a  synthetic  utility  number  between  0  and  1 :  the  available  node  bearing 
the  largest  utility  is  chosen. 
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Figure  4:  End-to-end  throughput  as  a  function  of 
traffic  generation  rate,  Poisson  model. 


Figure  5:  End-to-end  throughput  as  a  function  of 
traffic  generation  rate,  correlated  model. 


We  recall  that  routing  protocols  require  a  means  to  measure  not  only  the  overall  quality  of  a  node,  but 
also  of  the  link  over  which  the  communication  is  to  take  place.  In  our  protocol,  we  are  inherently 
provided  with  such  a  measure,  which  is  the  class  of  the  node.  We  have  previously  noted  [10]  that  a  fixed 
class  is  however  not  a  good  choice,  as  it  does  not  keep  track  of  a  node’s  actual  capability  to  exploit 
concurrent  access  as  allowed  by  the  MAC.  To  this  end,  we  deployed  a  class  update  policy  that  allows  to 
convey  more  precise  information  about  the  amount  of  spatial  multiplexing  the  node  can  support.  This  is 
key  in  MAC  operations:  only  by  knowing  the  exact  limits  of  its  receiver  does  a  node  make  correct 
decisions  as  to  the  traffic  to  inject.  Correct  MAC  operations  then  reflect  on  routing  (through  correct 
transmissions  and  through  more  correct  information  to  add  up  to  the  relay  utility).  To  this  end,  we 
restricted  the  class  of  a  node  to  within  the  set  {2,  4,  6,  8},  chose  the  success  ratio  as  a  measure  of  the 
reliability  of  a  link,  and  considered  the  following  link  class  updating  algorithm: 


N 


r  =  a  — 


ACK 


N 


+  (l-a)r 


old 


SEM 


where  ^oid  and  ^new  are  the  old  and  new  reliability  measures,  respectively,  and  ^  is  a  weight  used  to 
tune  the  reaction  speed  of  the  algorithm  to  link  reliability  changes.  The  new  reliability  is  then  compared  to 
different  increasing  and  a  decreasing  thresholds.  If  either  is  exceeded,  the  class  of  the  link  is  respectively 
increased  or  decreased,  though  always  allowed  to  have  values  in  the  set  {2,  4,  6,  8}.  Every  time  the  class 
is  changed,  the  link  reliability  variable  is  reinitialized.  In  terms  of  routing,  the  choice  of  updating  the  class 
dynamically  yields  apparent  benefits.  We  investigated  more  deeply  in  this  direction  to  yield  such  results 
as  those  in  Figures  7  and  6,  which  depict  end-to-end  throughput  and  success  ratio  for  different  class 
update  and  MAC -level  backoff  policies.  The  first  observation  in  order  here  is  that  the  more  aggressive 
Dest-Lock  policy  achieves  a  better  throughput  too.  In  a  multihop  network,  this  result  is  easily  explained: 
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as  a  node  fails  to  contact  a  node  for  routing,  it  is  allowed  to  try  another  one  in  the  following  frame.  Node 
lock  does  not  allow  so,  and  leaves  the  node  in  backoff  for  a  while  before  retrying.  However,  regardless  of 
the  specific  backoff  policy  (which  in  the  end  only  affects  the  general  throughput  trend),  we  see  from  the 
figures  that  updating  the  class  indeed  yields  a  benefit,  and  a  more  apparent  one  in  terms  of  success  ratio. 
Being  able  to  maintain  higher  levels  of  throughput  at  higher  success  ratio  indicates  that  fewer  resources 
are  wasted  and  that  only  transmissions  with  a  chance  to  succeed  actually  take  place.  This  is  very 
important  in  an  ad  hoc  network,  where  resources  may  be  limited  and  thus  must  be  carefully  exploited. 


lambda  [pkt/s] 
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Figure  7:  End-to-end  throughput  as  a  function  of  Figure  6:  End-to-end  success  ratio  as  a  function 

traffic  generation  rate,  for  the  class  update  and  of  traffic  generation  rate,  for  the  class  update  and 

fixed  class  policies  .  fixed  class  policies  . 


While  the  observations  above  make  very  important  points  towards  the  development  of  increasingly 
effective  routing  protocols,  special  attention  should  be  posed  to  PHY -originated  impairments  such  as 
those  due  to  imperfect  channel  estimation.  A  significant  fraction  of  Prof  Zorzi  efforts,  last  year,  was 
devoted  to  the  characterization  of  imperfect  channel  knowledge  through  the  statistics  of  the  matrix 
describing  the  deviation  of  channel  estimates  from  the  actual  realization.  In  ad  hoc  networks,  this  was 
shown  to  be  different  than  in  studies  of  point-to-point  links,  which  form  the  largest  part  of  the  available 
literature.  In  particular,  it  was  shown  that  the  channel  matrix  can  be  written  as 


H  =  H  +  JH  +  Z''  =  H  +  AH 


(1) 
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where  the  last  term  highlights  the  source  of  uncertainty  in  the  channel  estimation  procedure:  this  term 
contains  both  noise,  through  the  matrix  Z”,  and  a  component  related  to  the  instantaneous  channel 
realization,  JH.  The  matrix  J,  in  this  term,  describes  the  effect  of  simultaneous  transmissions  between 
other  node  pairs,  which  cause  a  “footprint”  of  the  signal  to  appear  in  the  matched  fdter  output  of  the 
current  receiver.  This  output  depends  on  the  channel  between  the  different  transmitters  and  the  receiver, 
which  explains  the  term  JH.  The  wrong  channel  estimates  in  (1)  can  easily  be  simulated  as  it  can  be 
shown  that  the  variance  of  the  elements  of  J  is  1/6A^,  where  N  is  the  length  of  the  training  sequence 
employed  for  channel  estimation.  The  elements  of  Z”  also  have  a  close-form  variance,  which  is  expressed 
as  No/NT,  where  No/2  is  the  double-sided  noise  power  spectral  density,  and  T  is  the  symbol  time.  Last 
year,  it  was  highlighted  that  acting  on  N  allows  to  drive  different  kinds  of  tradeoffs.  More  specifically,  by 
calling  Nsig  the  length  of  the  training  sequence  associated  to  signaling  messages,  and  Ndata  the  length  of 
the  sequence  used  for  data  packets,  then  one  can  choose  whether  to  move  the  system  working  point  on  a 
better  throughput,  better  efficiency,  better  success  ratio  mode,  or  to  a  combination  of  those.  This  year, 
these  tradeoffs  have  been  studied  to  a  larger  extent,  yielding  such  results  as  those  in  Figure  9,  which 
shows  the  contour  curves  of  average  transmission  efficiency  (the  ratio  of  purely  information  bits  sent  per 
frame  over  all  bits  sent  in  a  frame),  success  ratio  (the  ratio  of  correctly  transmitted  packet  data  units,  i.e., 
data  chunks  suitable  to  be  sent  over  one  antenna,  to  all  packet  data  units  sent)  and  delivery  delay  (time  in 
ms  required  to  deliver  one  packet  to  its  destination). 

The  superposition  of  these  curves  allows  to  understand  which  metrics  are  competing  and  which  can  be 
optimized  almost  simultaneously.  For  example,  from  Figure  9,  we  infer  that  success  ratio  and  efficiency 
are  competing  metrics:  to  opt  for  a  high  success  ratio  (e.g.,  to  reduce  retransmissions)  means  to  accept  a 
limited  transmission  efficiency,  because  the  large  Ndata  required  to  support  the  wanted  success  ratio  will 
exceedingly  increase  the  frame  time  duration.  With  the  help  of  tradeoff  curves  such  as  those  in  Figure  8, 
one  can  also  determine  the  sensitivity  of  network  metrics  to  local  variations  on  the  choice  of  Nsig  and 
Ndata.  In  this  figure,  which  compares  the  performance  in  terms  of  throughput  and  transmission  efficiency 
for  varying  Nsig  and  Ndata,  the  most  desirable  operating  point  is  as  close  as  possible  to  the  top-right 
corner.  In  this  case  we  see  that  increasing  or  decreasing  Nsig  for  fixed  Ndata  around  the  point  of 
optimality,  corresponding  to  traveling  the  bold  black  lines)  yields  worse  results  than  varying  Ndata  for 
fixed  Nsig  (which  corresponds  to  traveling  dashed  grey  lines).  In  fact  [12],  decreasing  Nsig  below  the 
optimality  point  makes  the  detection  of  RTS  messages  more  prone  to  errors:  if  this  happens,  links  cannot 
be  set  up.  On  the  other  hand,  if  Nsig  is  too  high,  RTSs  are  given  too  much  protection,  and  thus  their 
detection  will  likely  end  successfully.  Therefore,  more  links  will  be  set  up,  more  data  packets  transmitted, 
and  the  channel  estimation  accuracy  offered  by  Ndata  turns  out  to  be  too  low  to  support  the  increased 
amount  of  data  transmissions. 
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Figure  9:  Contour  curves  of  transmission 
efficiency,  PDU  success  ratio,  and  Delivery  delay 
(in  ms)  as  a  function  of  Nsig  and  Ndata.  Each 
curve  is  obtained  by  cutting  the  surface  plot  of 
each  metric  with  a  plane  corresponding  to  the 
metric  value  indicated  on  the  curve. 


Figure  8:  Tradeoff  curves  between  of 
transmission  efficiency  and  throughput.  Black 
arrows  correspond  to  variations  in  Ndata  for 
fixed  Nsig,  vice-versa  for  grey  arrows. 


Another  extension  to  last  year’s  work  includes  the  investigation  of  the  following  question,  which  directly 
arises  from  the  previous  evaluation:  since  the  number  of  packet  data  units  sent  in  a  frame  is  the  most 
important  variable,  and  perhaps  the  most  difficult  to  support  through  a  proper  choice  of  the  PHY 
parameters,  does  any  artificial  bound  on  the  number  of  PDUs  sent  on  a  link  help  at  all?  To  see  this,  we 
plot  throughput  against  traffic  in  Figure  10,  and  limit  the  number  of  antennas  to  be  used  in  a  transmission 
by  limiting  the  maximum  class  of  a  node.  So  while  in  the  previous  years  we  almost  always  used  a  {2,  4, 
8}  configuration  (meaning  that  a  node  can  be  transmitted  to  using  at  most  either  2,  4,  or  8  antennas),  now 
we  reduce  that  to  {2,  4,  7},  {2,  4,  6}  and  {2,  4,  5}.  These  configurations  limit  to  at  most  7,  6,  or  5 
antennas  the  maximum  spatial  multiplexing  that  a  transmit  node  is  ever  allowed.  As  we  see  from  the 
Figure,  though,  this  hardly  yields  any  advantage.  Although  there  is  a  slight  throughput  improvement,  this 
does  not  justify  the  operation  of  limiting  the  number  of  transmit  antennas.  For  comparison,  we  also  put  in 
the  picture  the  {1,  2,  4}  and  {3,  6,  12}  configuration,  which  assumes  that  the  number  of  physically 
available  antennas  equipped  at  each  node  changes  from  8  to  4  and  12,  respectively.  The  results  expectedly 
show  worse  performance  for  the  {1,  2,  4}  configuration  (4  antennas  equipped  at  each  node)  and  better 
performance  for  the  {3,  6,  12}  configuration  (12  antennas  at  each  node).  From  this  comparison  and 
previous  results,  we  infer  that  playing  on  the  tradeoff  arising  between  physical-level  parameters  and 
higher-level  performance  indicators  (while  allowing  the  upper  layers  to  automatically  manage 
communications)  is  a  much  better  way  to  increase  performance  than  to  empirical  limit  transmission 
parallelism.  While  these  are  just  examples  of  the  performance  investigations  carried  out  by  Prof.  Zorzi, 
the  reader  is  referred  to  [11],  [12]  and  [13]  for  a  more  in-depth  discussion. 
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Figure  10:  Throughput  as  a  fuuction  of  traffic  for 
differeut  uumher  of  auteouas  aud  class 
coufiguratioos. 


Routing  in  Cooperative  MIMO  Ad  Hoc  Networks 


In  the  recent  literature,  cooperation  has  been  envisioned  to  increase  the  capacity  of  wireless  networks. 
The  cooperative  rationale  provides  the  transmission  of  information  content  by  nodes  that  overheard  the 
transmission  by  the  original  source.  The  goal  may  be  either  to  increase  the  capacity  of  the  original 
communication  or  to  decrease  its  outage  probability. 


Most  of  the  previous  work  looks  at  cooperation  from  a  single-link  perspective.  In  general,  simple 
topologies  with  few  nodes  are  analyzed,  often  disregarding  practical  issues  such  as  nodes  availability  for 
cooperation,  packet  arrivals  and  interference. 


The  single-link  perspective  not  only  biases  the  performance  evaluation,  but  also  fails  to  enlighten  some 
aspects  of  cooperation  that  may  be  extremely  relevant  in  networking,  thus  potentially  limiting  the 
effectiveness  of  the  design.  We  therefore  look  at  cooperation  from  a  network-wide  perspective. 


In  particular,  we  consider  a  MIMO  ad  hoc  network  where  the  potential  of  the  physical  layer  is  exploited 
to  improve  the  parallelism  of  the  network,  rather  than  to  increase  the  capacity  of  each  single  link  and  the 
interference  rejection  capabilities  granted  by  the  MIMO  technology  is  used  to  let  multiple  transmitters  to 
concurrently  access  the  channel.  Thus,  the  access  protocol  does  not  follow  a  collision  avoidance 
rationale,  and  the  outage  probability  of  a  transmission  strongly  depends  on  the  interference  coming  from 
the  other  transmissions  taking  place  at  the  same  time. 
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In  our  previous  work,  we  addressed  the  design  of  cooperation  in  such  a  network,  considering  single-hop 
communications  only.  We  proposed  a  cooperative  protocol  that  integrates  error  control  and  cooperation  in 
order  to  increase  the  efficiency  of  the  communications. 


Due  to  the  start  and  the  end  of  interfering  transmissions,  the  conditions  of  the  channel  during  a 
communication  may  significantly  vary.  Therefore,  the  capabilities  of  the  communication  protocol  to  adapt 
the  transmission  rate  to  the  instantaneous  conditions  of  the  channel  and  the  evolution  of  the  interference  is 
a  key  design  goal.  Hybrid  Automatic  Retransmission  reQuest  (HARQ),  by  providing  packet  encoding, 
retransmission  and  combining,  has  been  shown  to  grant  adaptiveness  and  efficiency  to  wireless 
communications  in  time  varying  channels.  We  focus  in  particular  on  Incremental  Redundancy  HARQ, 
where  in  each  retransmission  of  a  packet,  the  source  transmits  a  different  fragment  of  a  long  codeword 
obtained  by  encoding  the  original  packet.  The  receiver  collects  the  different  fragments  and  combines 
them  in  order  to  get  a  code  with  lower  rate,  thereby  increasing  the  probability  of  successful  decoding. 


The  proposed  cooperative  protocol  allows  the  neighbors  of  a  transmitting  source  that  decoded  the  packet 
sent  by  the  latter  to  transmit  redundancy  in  its  place.  This  increases  the  channel  diversity  of  the 
communication  and  allows  an  improved  efficiency. 


In  the  considered  multiple  simultaneous  access  scenario,  cooperation  has  also  a  beneficial  effect  on  other 
ongoing  communications.  In  fact,  shorter  and  more  effective  communications  result  in  a  lower  per  packet 
interference.  Moreover,  the  possibility  for  simultaneous  transmissions  considerably  reduces  the  overhead 
needed  to  set  up  cooperation. 


The  cooperative  HARQ  protocol  described  above  has  been  shown  to  significantly  improve  the 
performance  of  the  network.  However,  by  shifting  the  design  from  single  link  to  the  entire  network 
perspective  (Fig.  11),  it  is  possible  to  achieve  a  deeper  integration  of  cooperation  in  the  system. 


In  particular,  cooperation  may  be  intended  as  a  means  to  improve  packet  forwarding  in  the  network, 
rather  than  to  improve  the  efficiency  of  a  single  link.  Herein,  we  propose  an  effective  protocol  that 
integrates  error  control  and  cooperation  with  multi-hop  routing  in  order  to  more  effectively  face  channel 
impairments  due  to  interference. 
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Figure  11:  Cooperative  communications  need  to  be  considered  in  a  network-wide  perspective. 


It  can  be  observed  that  traditional  routing  protocols  may  not  be  effective  in  a  network  where  channel 
conditions  quickly  and  unpredictably  vary.  In  fact,  route  formation  cannot  deal  with  the  quick  oscillations 
of  the  channel  quality,  and  may  fail  to  select  an  efficient  sequence  of  relays  to  reach  the  destination. 


In  order  to  cope  with  these  issues,  we  propose  a  cooperative  protocol,  referred  to  as  cooperative  HARQ 
and  opportunistic  routing  (CHOR)  protocol  in  the  following,  that  manages  packet  forwarding  by 
adaptively  and  opportunistically  selecting  an  efficient  route  to  the  destination.  We  underscore  again  that 
an  inefficient  route,  where  each  hop  requires  long  delivery  attempts,  i.e.,  attempts  with  many  transmission 
phases,  may  increase  the  interference  load  of  the  network  and  damage  other  communications  by 
increasing  both  their  failure  rate  and  their  waiting  time  before  packet  transmission. 


The  CHOR  protocol  adds  two  main  features  to  the  basic  communication  protocol: 

it  enables  other  nodes  that  overheard  the  transmission  by  the  source  to  send  fragments  of  the  codeword 
obtained  by  re-encoding  the  received  packet 

it  opportunistically  switches  to  a  different  next  hop  relay  in  order  to  minimize  the  amount  of  redundancy 
sent  and  to  guarantee  a  geographical  advancement  to  each  packet. 


In  the  cooperative  error  control  we  previously  proposed  the  cooperators  only  transmit  incremental 
redundancy  to  strengthen  the  chosen  route  to  the  destination.  In  the  cooperative  integrated  error  control 
and  routing  protocol  those  nodes  may  opportunistically  take  charge  of  packet  forwarding  in  order  to  grant 
an  efficient  advancement  of  the  packet  toward  the  destination  (see  Fig.  12). 
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Figure  12:  Cooperative  rationale  plainly  applied  to  multi-hop  communications. 


As  in  the  cooperative  error  control  protocol,  the  possibility  for  simultaneous  channel  access  granted  by 
the  MIMO  physical  layer  allows  an  efficient  deployment  of  the  cooperative  relay  selection. 


Figure  13:  State  diagram  of  the  proposed  cooperative  packet  forwarding  protocol. 


A  state  diagram  of  the  protocol  is  shown  in  Fig.  13.  The  original  transmitter  quits  the  communication  if 
the  destination  reports  bad  channel  conditions  and  there  are  no  available  cooperators.  Otherwise,  the 
original  transmitter  performs  transmission  of  adaptive  incremental  redundancy  possibly  helped  by  its 
neighbors.  If  a  nodes  with  a  smaller  number  of  hops  decodes  the  packet,  then  the  original  communication 
is  quit  and  the  new  node  takes  charge  of  packet  forwarding. 
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We  developed  a  simulator  which  describes  the  network  from  the  physical  to  the  routing  level. 
Performance  is  averaged  over  random  topologies,  where  nodes  are  uniformly  placed  in  a  fixed  square  area 
of  side  500m.  Packet  arrivals  in  the  buffer  of  each  node  are  generated  according  to  a  Poisson  process,  and 
the  destination  of  each  packet  is  randomly  chosen  among  all  the  other  nodes  in  the  network.  We  assume 
that  a  fixed  fraction  of  the  queue  is  left  for  packets  received  from  other  nodes  which  must  be  forwarded  to 
their  final  destinations.  A  timeout  is  also  set,  and  packets  which  have  not  yet  been  delivered  when  it 
expires  are  dropped.  Finally,  each  node,  as  previously  explained,  has  a  routing  table,  generated  via 
Shortest  Path  algorithm,  in  which  the  routes  towards  every  possible  destination  are  stored. 


The  physical  layer  is  based  on  a  Minimum  Mean  Square  Error  (MMSE)  version  of  the  layered  space-time 
multiuser  detection  (EASTMUD)  architecture  previously  considered  in  this  research  line. 


Node  number 


Figure  14:  End  to  end  throughput  for  the  considered  scheme. 


Fig  14  depicts  the  end— to— end  throughput  for  the  basic,  OR  and  CFIOR  protocols  as  a  function  of  the 
number  of  nodes  in  the  network,  each  equipped  with  only  one  antenna.  As  can  be  seen,  the  throughput 
grows  for  all  the  protocols  as  the  number  of  nodes  increases,  until  it  reaches  a  floor. 

In  fact,  as  the  number  of  nodes  is  increased,  the  network  is  initially  able  to  support  the  greater  traffic 
generated  by  the  transmissions,  whereas  beyond  a  certain  value  the  level  of  parallelism  that  the  network 
can  sustain  saturates.  This  is  due  to  the  effect  of  interference,  that  does  not  allow  other  communications  to 
start  (the  intended  receiver  reports  bad  channel  conditions)  and  results  in  an  increased  failure  rate  of  the 
delivery  attempts. 


The  OR  protocol  achieves  an  improved  throughput  with  respect  to  that  of  the  basic  protocol.  In  fact, 
although  a  slightly  larger  signaling  overhead  is  required,  the  capability  of  choosing  opportunistically  the 
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path  improves  the  probability  that  the  packet  is  eventually  delivered  to  its  destination.  Moreover,  the 
avoidance  of  bad  links  enables  more  efficient  communications,  and  thus  the  interference  load  decreases 
and  each  packet  is  delivered  in  a  shorter  time.  The  CHOR  protocol  outperforms  both  the  OR  protocol  and 
the  basic  communication  scheme.  The  possibility  for  the  neighbors  with  a  good  channel  to  the  next  hop 
relay  to  help  the  HARQ  process  by  sending  HARQ  packets  on  behalf  of  the  source  speeds  up  the  delivery 
process  and  reduces  the  probability  that  nodes  enter  backoff.  However,  cooperative  HARQ  may  result  in 
bursty  transmissions  if  the  number  of  HARQ  relays  is  large,  an  event  that  may  happen  in  densely 
populated  networks.  The  burst  of  cooperative  transmissions  may  generate  short-term  interference 
congestion  conditions,  that  may  be  harmful  for  the  other  ongoing  communications.  The  effectiveness  at 
the  intended  receiver  depends  on  the  tradeoff  between  the  amount  of  collected  redundancy  and  the  quality 
of  the  channel  experienced  by  each  HARQ  packet. 

As  a  final  remark,  we  underscore  the  importance  to  extend  the  cooperative  rationale  to  a  wider  network 
perspective,  and  to  consider  the  effect  of  the  surrounding  network  on  each  single  cooperative 
communication,  that  are  distinctive  features  of  our  work 
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SUMMARY  OF  ACCOMPLISHMENTS 


Our  primary  focus  is  on  exploiting  the  PHY-layer  gains  due  to  MIMO  at  the  higher  layers.  In 
line  with  this,  we  have  two  researeh  projeets.  These  projeets  target  the  exploitation  of  MIMO 
spatial  diversity  and  MIMO  spatial  multiplexing,  respeetively.  In  the  following  we  summarize 
our  aehievements  in  these  projeets. 


1.  Abstractions  of  MIMO  Diversity  Gain:  Our  purpose  in  this  project  is  to  assess  the  fidelity  of 
previously  considered  abstractions  (or  models)  of  MIMO  diversity,  by  comparing  the 
performance  predicted  by  these  models  with  the  practical  performance  observed  on  our  testbed. 
In  the  previous  report  we  have  included  the  description  of  the  models  we  have  investigated,  and 
some  of  our  observations.  In  this  report,  we  include  our  additional  experiments  and  results  with 
regards  to:  (i)  assessing  the  fidelity  of  network  simulations  (ii)  evaluating  how  the  routing  layer 
decisions  interact  with  each  model  when  diversity  gain  is  examined  on  multi-hop  flows.  For  (i), 
we  include  our  observations  on  single-link  flows,  and  for  (ii)  we  include  our  observations  when 
different  routing  metrics  are  used. 

2,  Medium  access  control  for  Multi-User  MIMO  Communications:  Our  purpose  in  this  project  is  to 
utilize  the  interference  cancellation  capability  using  MIMO  systems  and  to  increase  the  spatial 
reuse  in  a  multi-hop  network  by  allowing  multiple  concurrent  communications  in  a  region. 
Towards  this  goal  we  have  designed  centralized  and  distributed  algorithms  that  divide  the 
communication  links  into  separate  sub-topologies,  which  are  then  activated  in  a  TDMA  fashion. 
In  a  given  sub-topology,  the  concurrent  communications  are  successful  with  high  probability.  We 
also  show  that  previously  employed  models  for  multi-user  MIMO  communications  yield 
significantly  low  probabilities  of  successful  decoding  in  multi-hop  settings. 
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1.  CAPTURING  THE  PERFORMANCE  GAINS  DUE  TO  MIMO 
DIVERSITY 


In  this  project,  we  focus  on  the  modeling  of  MIMO  diversity  gain.  We  have  explained  the  details 
of  this  project  in  the  previous  progress  report.  Briefly,  we  examine  the  cross  layer  interactions 
between  the  PHY,  link  and  routing  layers  in  networks  with  MIMO  links  operating  in  the 
diversity  mode.  Towards  this,  we  first  create  models  that  represent  the  MIMO  diversity  gain  in 
terms  of  a  mapping  between  the  SNR  and  the  achievable  BER  (or  PER,  i.e.,  packet  error  rate).  In 
parallel,  we  also  take  measurements  on  an  indoor  802.1  In  testbed.  Second,  we  examine  the 
goodness  of  each  model  by  comparing  their  performance  with  the  performance  driven  by 
measurements,  on  multihop  networks  via  simulations.  Our  work  brings  to  light  the  following 
cross-layer  dependencies  that  affect  the  gains  due  to  MIMO  diversity:  (1)  PHY  layer  gains  due  to 
MIMO  diversity  do  not  always  carry  over  to  the  higher  layers,  (2)  the  use  of  other  PHY  layer 
features  such  as  PEC  codes  significantly  influence  the  gains  due  to  MIMO  diversity,  and  (3)  the 
choice  of  the  routing  metric  can  impact  the  gains  possible  with  MIMO. 


In  the  previous  progress  report,  we  have  shown  results  with  regards  to  the  observations  (1)  and 
(2)  above.  During  this  period,  we  have  added  more  results  that  validate  the  fidelity  of  our 
simulations  for  assessing  a  model  in  terms  of  representing  the  MIMO  diversity  gain.  These 
results  are  presented  in  part  A  below.  We  have  also  evaluated  these  models  using  different 
routing  metrics  in  multi-hop  networks.  Our  observations  are  presented  in  part  B. 


A.  Results  on  Single-Hop  Flows: 

In  order  to  assess  our  multi-hop  network  simulations  in  terms  of  how  good  they  are  for 
comparing  the  end-to-end  performance  of  each  model  to  that  with  the  measurement-driven 
model  (MDM),  we  compare  the  throughput  results  in  our  simulations  directly  with  those 
obtained  via  measurements.  Towards  this,  we  first  measure  the  link  throughput  with  SISO  and 
MIMO  on  different  links  in  our  testbed  at  all  three  rates.  Then,  we  simulate  several  single-hop 
topologies  and  obtain  the  throughput  with  SISO  and  MIMO  for  a  range  of  SNR  values.  We 
present  the  results  in  Pig.  1  for  two  sample  SNR  values.  We  observe  that  the  simulated 
throughput  with  MDM  differs  by  at  most  11%  from  what  we  measure  on  our  testbed.  This  result 
suggests  that  MDM  is  an  effective  indicator  of  throughput;  thus,  MDM  is  considered  to  be  a 
good  benchmark  to  quantify  the  throughput  in  multi -hop  settings.  We  also  observe  that  the 
throughput  with  M-3P  (SP)  differs  by  at  most  21%  from  the  throughput  with  MDM,  and  24% 
from  the  throughput  on  the  testbed  (at  3  bps/Hz  data  rate).  On  the  other  hand,  the  use  of 
simplistic  models  can  lead  to  a  drastic  overestimation  in  the  performance  improvement  with 
MIMO.  Por  instance,  with  M-2,  the  absolute  throughput  is  much  lower  (77%  lower  with  MIMO 
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and  3,000  times  lower  with  SISO)  than  what  is  observed  on  the  testbed.  Since  M-2  does  not 
account  for  FEC  codes,  the  packet  transmissions  with  this  model  experience  higher  losses 
leading  to  poor  throughput.  This  effect  is  more  pronounced  with  SISO  than  MIMO  (in  fact  the 
throughput  is  close  to  zero,  see  Figures  3  and  4)  and  thus,  the  diversity  gain  is  extremely  high. 
With  the  fast  fading  models,  M-3(FF)  and  M-3P(FF),  again  the  improvements  with  MIMO  over 
SISO  are  much  higher  than  what  is  observed  on  the  testbed.  This  is  because  the  errors  are  less 
bursty.  While  the  FEC  codes  are  more  effective  than  they  are  in  practice,  the  errors  tend  to  be 
spread  over  many  more  packets,  leading  to  poorer  throughputs  with  these  models.  With  M-3P 
(SF),  the  impact  of  the  environment  in  the  indoor  setting  (bursty  errors)  is  better  captured;  thus, 
the  behavioral  results  are  similar  to  what  is  observed  on  the  testbed. 

Note:  however  that,  the  throughput  behavior  with  the  models  are  not  exactly  identical  to  that 
with  measurements,  because  the  wireless  channel  could  be  different  in  reality  from  what  is 
implemented  in  the  simulator. 
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Figure  3:  SISO 
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Figure  4:  MIMO 
Throughput  at  22dB 


B,  Results  on  Multi-Hop  Flows  Using  Different  Routing  Metrics 

Next,  we  perform  simulations  on  multi-hop  networks  using  each  abstraction  at  the  physical  layer; 
we  do  so  on  multi-hop  flows  where  the  routing  protocol  (DSR  in  our  work)  chooses  the  shortest 
or  the  minimum-ETX  route  for  packet  delivery.  The  first  policy  refers  to  the  DSR  version  that 
uses  hop  count  as  the  routing  metric.  In  this  policy,  DSR  chooses  the  minimum  hop  route  from 
among  those  that  are  discovered.  When  this  route  fails,  other  cached  routes  may  be  used.  With 
the  second  policy,  we  compute  the  ETX  metric  to  be  the  inverse  of  the  link’s  reliability.  In 
particular,  ETX  is  quantified  in  terms  of  the  packet  delivery  rate  (PDR)  on  the  link  for  the 
corresponding  rate.  Intermediate  nodes  are  not  allowed  to  report  cached  routes  to  the  source  as  in 
the  default  ETX  implementation. 


42 


We  examine  the  ratio  of  the  achievable  end-to-end  network  throughput  using  MIMO  links  to  that 
using  SISO  links,  with  ETX-based  routing.  We  observe  that  with  ETX,  this  ratio  decreases  as 
compared  to  the  minimum  hop  routing  considered  earlier.  We  define  D  to  be  the  reduction  in  this 
ratio  compared  to  minimum  hop  routing  (Eig.  7).  We  consider  the  same  50-node  topologies  for 
evaluating  both  cases.  As  seen  in  the  figure,  the  gains  due  to  MIMO  decrease  with  the  use  of 
ETX  by  as  much  as  20%.  With  ETX,  the  use  of  the  more  stable  (albeit  longer)  routes  with  more 
reliable  links  (higher  average  SNR)  improves  the  SISO  performance  as  compared  to  minimum- 
hop  routes.  With  the  higher  average  SNRs  on  the  SISO  links,  the  EEC  codes  are  more  effective 
and  MIMO  does  not  provide  significant  additional  benefits.  Similarly,  we  observe  that  the 
conformance  of  the  M-3P  (SE)  to  MDM  in  terms  of  the  delay  performance  increases;  the  two 
curves  (corresponding  to  M-3P(SE)  and  MDM)  are  closer  to  each  other  in  Pig.  6  than  in  Pig.  5. 
This  is  because  ETX  links  operate  at  relatively  high  SNRs  where  the  PER  with  M-3P  (SP)  is 
closer  to  that  with  measurements. 
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Figure  7:  Performance  improvements  due 
to  MIMO  relative  to  SISO  decrease  with 
ETX-hased  routing 


This  work  is  currently  under  submission  to  ACM  CONEXT  2009. 


2.  MUETI-USER  MIMO  COMMUNICATIONS 

Space  division  multiplexing  with  MIMO  systems  allows  nodes  to  transmit  or  receive  multiple 
streams  at  a  time.  In  this  project,  our  goal  is  to  efficiently  utilize  spatial  multiplexing  in  multi¬ 
hop  multi-user  MIMO  networks.  We  have  stated  our  goal  and  approach  in  detail  in  last  year’s 
report.  In  brief,  we  design  a  topology  control  mechanism  that  enables  successful  multi-stream 
receptions  in  a  multi-user  MIMO  network. 

A,  Re-visiting  the  Prohlem 
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First,  we  briefly  re-state  our  problem.  In  multi-hop  MIMO  networks,  wherein  the  nodes  transmit 
using  one  of  their  antennas  (this  choiee  being  faeilitated  by  selective  diversity)  and  utilize 
Successive  Interference  Cancellation  (SIC)  for  receiving  their  target  signals  in  presence  of 
interference,  the  MU-SIC  (Successive  Interference  Cancellation  in  Multi-User  MIMO  networks) 
problem  seeks  to  find  a  maximal  number  of  communications  that  can  be  simultaneously  active, 
while  each  active  receiver  can  extract  its  target  signal  successfully. 

There  are  two  challenges  towards  realizing  the  above-mentioned  system.  First,  a  receiver  can 
afford  no  more  than  A-1  strongly  interfering  concurrent  transmissions  with  its  intended 
transmission  (where  A  is  the  size  of  the  antenna  array).  Furthermore,  the  success  of  the  detection 
of  each  strong  interfereing  signal  is  determined  based  on  its  SINR,  which  itself  depends  on  the 
strength  of  all  signals  in  the  received  compound  signal.  To  address  these  challenges,  or 
specifically,  to  ensure  that  SIC  can  successfully  be  used  at  the  receivers  to  remove  unwanted 
interference,  a  higher  layer  mechanism  must  control  which  links  can  be  active  at  the  same  time. 
Finally,  this  mechanism  must  allow  every  transmitter  to  acquire  medium  access  as  frequently  as 
possible,  in  order  to  limit  medium  access  delays. 

We  show  that  the  MU-SIC  problem  is  in  the  NP-hard  class. 

B,  Centralized  Solution  Approach  to  the  MU-SIC  Prohlem:  C-MUSIC 

We  design  C-MUSIC,  a  centralized  heuristic  algorithm  that  groups  the  communication  links  in  a 
given  multi-hop  network  into  a  minimal  number  of  groups,  such  that  in  each  group  SIC  can  be 
carried  out  successfully  with  an  arbitrarily  small  error  probability,  5. 

C-MUSIC  takes  as  input  the  interference  graph  G'  (V,  E')  of  the  network,  in  which  there  is  a 
vertex  in  V  for  every  directional  communication  in  the  network  and  there  is  a  directional  edge  in 
E'  between  each  node  pair  that  interferes  with  each  other.  (We  have  described  the  construction  of 
the  interference  graph  in  detail  in  prior  progress  reports;  therefore,  we  skip  the  details  here.)  It 
maintains  a  dynamically  changing  set  of  nodes  R'  which,  initially  contains  the  nodes  in  G'  (i.e., 
V).  It  executes  in  steps;  in  the  first  step  it  finds  a  maximal  subset  Rf  of  R',  a  maximal  subset  of 
G',  whose  elements  can  be  active  concurrently.  Rf  is  subtracted  from  R'  and  in  the  subsequent 
iterations,  a  maximal  subset  of  R'-  Rf  is  found.  C-MUSIC  terminates  when  R'  is  empty.  This 
construction  is  depicted  in  pseudocode  in  Eigure  8  below. 
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Figure  8:  Demonstrating  how  every  link  is  visited  and  evaluated  before  the  decision 
whether  it  must  be  added  to  a  group. 


We  show  that  C-MUSIC  has  a  performance  efficiency  of  Q,  where  Q  is  the  maximum 
opportunity  cost  in  the  network.  (The  opportunity  cost  Oe  for  link  e  is  the  maximum  number  of 
links  that  interfere  with  e,  which  can  be  active  simultaneously  if  e  is  deactivated.  From  among 
the  opportunity  costs  of  all  the  links  in  the  network,  Q  is  the  maximum  opportunity  cost.) 


C.  Distributed  Solution  Approach  for  MU-SIC  Problem:  D-MUSIC 

We  also  propose  a  distributed  variant  of  the  C-MUSIC  algorithm.  With  D-MUSIC,  nodes 
exchange  local  messages  and  determine  the  links  that  can  be  activated  concurrently  in  their 
neighborhood.  Thus,  while  the  sub-topologies  are  determined  on  a  global  scale,  the  process  of 
determination  is  performed  in  a  distributed  manner,  in  parallel,  and  locally  in  the  different  parts 
of  the  network. 

Each  node  associates  itself  with  one  of  the  Transmitter,  Receiver,  and  Overhearer  functionalities 
in  each  sub-topology.  Each  node  aims  to  become  a  transmitter  for  each  of  its  neighbors,  in  at 
least  one  sub-topology.  The  link  construction  is  initiated  at  the  Transmitter  nodes  (as  shown  in 
Eigure  9  below).  Upon  hearing  a  link  request  from  a  neighbor  u,  the  target  node  v  computes 
whether  the  candidate  transmitter  (u)  can  be  decoded  given  the  accumulated  level  of  interference 
in  the  network;  the  accumulated  interference  is  caused  by  the  communications  that  have  already 
become  part  of  this  sub-topology.  Towards  this  calculation,  v  considers  existing  transmitters  in 
terms  of  whether  they  are  received  with  a  stronger  or  weaker  power  compared  to  u.  If  link  (u,v) 
can  be  formed,  v  becomes  a  Receiver  and  acknowledges  the  link  request;  otherwise  it  remains  an 
overhearer  and  sends  a  negative-acknowledgement  to  u.  These  steps  are  shown  in  Eigures  10  and 
1 1  below. 
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D,  Performance  of  C-MUSIC  and  D-MUSIC 

We  evaluate  the  performanees  of  C-MUSIC  and  D-MUSIC  in  various  settings  using  OPNET 
network  simulator.  We  also  compare  the  performance  of  our  approaches  with  an  approach 
where,  A-1  other  links  are  allowed  to  be  simultaneously  active  in  a  neighborhood  of  a  receiver; 
as  discussed  earlier,  many  MAC  protocols  are  based  on  allowing  these  many  concurrently  active 
links  [2],  We  call  this  the  degree  of  freedom  based  approach  or  simply  DoF-based  topology 
control. 


In  Figure  12,  we  show  that  the  frequency  with  which  sub-topologies  are  activated  with  D- 
MUSIC  is  high  at  moderate  densities.  In  Figure  13,  show  that  topology  control  with  D-MUSIC 
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in  dense  multi-user  MIMO  networks  effeetively  bounds  the  probability  of  deeoding  error.  In 
Figure  14,  we  eompare  C-MUSIC  and  D-MUSIC  with  the  DoF-based  topology  eontrol.  We  find 
that,  although  the  DoF-based  topology  eontrol  finds  fewer  sub-topologies  (sinee  it  only  eonsiders 
the  number  of  interferers),  the  throughput  is  mueh  lower  compared  to  what  is  achieved  with  both 
our  algorithms.  This  is  because  the  decoding  error  probability  in  each  sub-topology  is  as  high  as 
0.7,  whereas  it  is  0.04  with  MUSIC.  Finally,  in  Figure  15  we  show  that  in  grid  topologies,  C- 
MUSIC  is  coerced  into  constructing  a  significantly  larger  number  of  sub -topologies  (2.2  times 
more  than  with  random  topologies  for  A=8),  each  consisting  of  much  fewer  links.  This  is 
because  in  a  grid  topology  each  node  is  equi-distant  from  its  neighbors,  and  the  signals  of 
neighbors  arrive  at  a  receiver  with  similar  powers.  Unlike  C-MUSIC  however,  we  observe  that 
DoF  constructs  a  similar  number  of  sub-topologies  with  both  the  random  and  the  grid 
deployments.  This  is  because  DoF  takes  into  account  the  node  degree  (which  is  7  on  average  in 
all  simulated  topologies)  and  not  the  SINR. 

This  work  is  currently  under  submission  to  IEEE  Infocom,  2010. 
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As  a  part  of  the  MURI  project.  Pis  Javidi,  Milstein  and  Cruz  have  used  sophisticated  methods 
from  stochastic  optimization  and  approximation  to  investigate  and  propose  optimal  cross  layer 
(PHY,  MAC,  and  Routing)  and  distributed  mechanisms  for  MIMO  ad  hoc  networks  [12],  [13], 
[14].  The  outcome  of  these  analytic  studies  resulted  in  the  design  of  a  family  of  scalable 
schedulers  which,  asymptotically,  use  minimum  total  power  to  transfer  the  information  bits  for 
all  the  end-to-end  connections  at  the  requested  rates.  To  this  end,  we  decomposed  the  multi-hop 
MIMO  network  into  multiple  “decoupled”  MIMO  broadcast  subsystems.  In  each  time  slot,  each 
subsystem  independently  decides  the  transmit  power,  the  antenna  weights,  the  transmission  rates 
and  the  forwarding  rules  of  information  bits  based  on  the  channel  state  information  and  previous 
decisions  of  other  neighboring  subsystems.  Given  a  decomposition  of  the  network  into  MIMO 
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“decoupled”  sub-blocks,  the  proposed  algorithms  are  shown  to  converge  asymptotically  to  the 
optimal  solution.  Even  though  asymptotic  in  nature,  these  algorithms  have  significant 
implications  for  the  design  of  MIMO  networks  as  they  1)  established  a  promising  benchmark, 
predicting  significant  improvements  for  many  network  scenarios,  and  2)  provide  insight  in  the 
design  of  cross-layer  optimal  resource  allocation  at  the  MIMO -MAC  and  routing  layers. 

Lott  and  Teneketzis  [8]-[9]  were  among  the  first  researchers  who  introduced  opportunism  in  the 
context  of  wireless  multi-hop  routing.  Independently,  Larssen  proposed  a  fully  implemented 
protocol  suite  for  Selection  Diversity  Forwarding  (SDF)  in  [10].  In  both  of  these  papers,  the 
source  identifies  a  set  of  potential  forwarders  and  multicasts  the  message  to  them.  The  successful 
recipients  respond  with  ACKs  in  the  order  in  which  they  were  listed  in  the  original  message 
header,  so  that  a  collision  between  ACKs  is  avoided.  The  source  identifies  the  best  amongst 
these  receivers,  according  to  some  predefined  criteria,  and  sends  a  forwarding  order  to  it.  Hence, 
the  actual  routing  decision  is  made  after  transmission  of  data.  While  [1]  suggests  time-invariant 
critera  such  as  hop-count,  transmission  rate,  and  transmission  energy  for  making  the  final 
forwarding  decision,  [10]  suggests  various  time-varying  criteria,  like  forward  progress,  cost 
progress,  queue  backlogs  etc,  for  making  the  final  forwarding  decision.  Extremely  Opportunistic 
Routing  (ExOR)  [11]  is  based  on  the  same  idea,  with  differences  in  the  selection  metrics. 

With  a  slight  twist,  [19]  considers  the  issue  of  diversity  routing  from  the  perspective  of  traffic 
stability,  i.e.,  guaranteeing  a  certain  vector  of  flow  rates  between  given  source-destination  pairs. 
We  would  also  like  to  point  out  that  our  work  in  [14]  philosophically  coincides  with  formulation 
in  [19],  with  the  main  point  of  difference  in  our  approaches  and  methodologies:  Our  work  [14] 
relies  on  constrained  convex  optimization  and  stochastic  approximation  frameworks  to  minimize 
(asymptotically)  the  expected  sum  power  subject  to  a  constraints  on  the  minimum  long-term 
average  rate  at  the  flow  level,  while  [19]  uses  a  Lyapunov-type  penalty  function  as  the  cost  of 
violating  the  QoS  requirements  to  establish  stability.  The  Lyapunov  drift  analysis  is,  then, 
applied  to  prove  the  stability,  while  working  within  a  stochastic  approximation  framework,  our 
formulation  allows  for  a  more  general  optimization  setting.  The  price  of  this  more  powerful 
methodology  is  the  asymptotic  nature  of  the  result.  However  our  more  recent  work  [20]  uses  the 
insight  obtained  from  [14]  to  combine  the  benefits  of  [8]  and  [19]  vis-a-vis  a  congestion  diversity 
measure.  Such  successful  merging  of  cross-layer  diversity  routing  schemes  underlines  the 
significance  and  promise  of  the  work. 

The  decision  per  hop  approach  raises  a  number  of  protocol  issues  [20]-[21]: 

Queuing  Delay  and  Protocol  Dynamics:  In  order  to  make  a  local  decision  on  the  best  next  hop 
node,  the  queue  backlogs  as  well  as  channel  quality  at  the  receiving  nodes  needs  to  be  fed  back 
to  the  transmitting  node.  In  other  words,  the  next  “best”  hop  a  dynamically  changing  option,  not 
only  based  on  channel  characteristics  and  variations  but  also  based  on  traffic  matrix.  This  can 
potentially  cause  oscillatory  behavior,  whose  long-term  impacts  might  be  negligible  but  cause 
undesirable  short  term  behavior.  In  other  words,  while  these  oscillatory  effects  tend  to  smooth 
out  for  long  flows,  they  significantly  disturb  the  behavior  of  short  flows. 
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Fig.  1.  A  multi-hop  wireless  network 

Congestion  Problems  and  Delay-Throughput  Trade-off:  In  many  cases,  there  is  an  inherent 
trade-off  between  packets  finding  the  “shortest”  route  to  the  destination  and  the  least  congested 
one.  We  have  been  able  to  show  examples,  where  choosing  “shortest”  routes  can  substantially 
decrease  the  delay  [15].  We  have  shown  that  including  queuing  delay  as  well  as  channel  quality 
(both  instantaneous  state  as  well  as  Baysian  models)  in  determining  the  “distance”  to  a  node  is 
necessary  to  guarantee  an  overall  desirable  behavior.  In  this,  our  primary  results  [15]  combine 
results  on  backpressure  routing  [17]-[19]  with  appropriate  concepts  from  opportunistic  routing 
to  improve  the  delay  performance  in  SISO  networks.  Our  work  attempts  at  taking  advantage  of 
diversity  in  congestion  levels  across  the  network  without  creating  unnecessary  balancing. 
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Fig.  2.  Delay  versus  traffic  rate  (intensity)  in  the  example  given  in  Fig.  1. 

Channel  Quality  Indicator  (CQI)  and  MIMO  scheduling:  Our  work  in  cross-layer  MIMO 
networking  relies  heavily  on  the  notion  of  channel  quality  to  coordinate  MAC  and  routing  layer 
functionalities.  In  particular,  we  have  relied  heavily  on  the  simplicity  of  our  MIMO-BC  sub¬ 
blocks  to  identify  and  adapt  instantaneous  rate  and  power  schedules.  The  choice  of  sub-gradient 
in  this  setting  was  key  to  obtaining  the  result. 

Channel  Quality  Distribution  (CQD):  In  order  to  make  a  local  decision  on  the  best  next  hop 
node,  the  channel  quality  information  is  needed.  This  potentially  includes  short  term  indicator 
(CQI)  as  well  as  the  channel  radio  characteristics.  While  the  opportunism  reduces  the  impact  of 
channel  estimation  error  in  CQI,  the  performance  of  opportunistic  algorithms  shows  a  large 
degree  of  sensitivity  to  CQD  [7].  An  important  extension,  our  work  [22]  addresses  the  question 
of  how  to  account  for  temporal  variation  in  channel  quality  as  well  as  effective  schemes  to 
construct  useful  probabilistic  models  to  capture  the  impact  of  diversity.  In  other  words,  given  the 
stationarity  of  CQI,  one  can  devise  many  algorithms  to  account  for  receiver  diversity  in  the  long 
term. 
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Fig.  3.  Delay  versus  the  size  of  the  “hole”  A-A,  henee,  network  size. 


High-SNR  Cross-Layer  Optimization  for  Delay-Limited  Data  Transmission  in 
Outage-Limited  Wireless  Systems 

In  our  previous  work  [7],  we  studied  the  optimal  operating  point  in  point-to-point  (quasi¬ 
static)  MIMO  channel  and  the  asymptotic  expression  of  the  total  probability  of  bit  error,  where 
errors  occur  either  due  to  delay  or  due  to  erroneous  decoding.  The  problem  setting  focused  on 
the  case  where  there  is  no  channel  state  information  at  the  transmitter  (no  CSIT)  and  no 
feedback,  and  on  the  static  case  of  fixed  operating  parameters.  In  our  recent  work  [2],  we 
generalize  this  study  to  include  other  outage-limited  communication  settings,  such  as  cooperative 
relay  [4,5]  and  fast- fading  channels  [3].  In  these  settings,  we  are  interested  in  the  asymptotic 
high-SNR  error  performance  when  the  delay  bound  requirement,  D,  is  finite  and  small.  Given 
that  the  asymptotic  expression  of  the  total  probability  of  bit  error  is  valid  without  requiring 
asymptotically  large  D,  it  is  then  meaningful  to  ask  about  the  optimal  coding  block  duration,  a 
question  which  is  not  answered  in  our  previous  studies  with  asymptotic  D.  To  have  a  valid 
expression  of  the  delay  violation  probability  at  finite  and  small  D,  we  assume  a  class  of  smoothly 
scaling  (with  SNR)  bit-arrival  processes.  This  class  of  processes  covers  many  interesting 
processes  used  for  traffic  modeling. 


Our  analysis  provides  closed-form  expressions  for  the  error  performance,  as  a  function  of 
the  channel  and  source  statistics.  These  expressions  identify  the  scaling  regime  of  the  source  and 
channel  statistics  in  which  either  delay  or  decoding  errors  are  the  dominant  cause  of  errors,  and 
the  scaling  regime  in  which  a  prudent  choice  of  the  coding  duration  and  rate  manages  to  balance 
and  minimize  these  errors.  That  is,  in  this  latter  regime,  such  optimal  choice  manages  to  balance 
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the  effect  of  channel  atypicality  and  burstiness  atypicality.  We  apply  the  results  in  the  different 
communication  settings  discussed  above. 


Buffer 

Wireless 

Fading 

Channel 

Source 

Encoder 

Decoder 

Destination 

Serve  RT  bits 
every  T  timeslots 


Bit  error  Bit  error 

(decoding)  (delay  violation) 


Fig.  4.  System  model 


Fig.  5.  Exponent  of  the  total  error  probability  for  SISO  fast  fading  chaimel  and  compound 

Poisson  bit-arrival  process 


Many-Sources  Large  Deviations  for  Max-Weight  Scheduling 

This  work  is  an  extension  of  our  previous  study  for  MIMO  multiple  access  channel  [6] 
and  our  recent  work  described  earlier.  In  [6],  we  studied  the  joint  optimization  of  the  MAC  layer 
and  the  physical  layer.  We  formulated  and  analytically  derived  bounds  on  the  optimal  operating 
point  and  the  asymptotic  (high-SNR  and  large  delay  bound  D)  error  performance  of  MIMO- 
MAC  channel  for  bursty  sources  with  delay  constraints.  The  adopted  system  model  brings 
together  the  four  types  of  gains:  diversity,  spatial  multiplexing,  space-division  multiple-access, 
and  statistical-multiplexing  gains. 

To  extend  the  work  in  [6]  to  the  case  when  the  delay  bound  D  is  not  asymptotically  large, 
in  [1]  we  look  at  the  queuing  performance  of  the  dynamic  queue -based  (Max-Weight)  scheduling 
policy  of  a  single  fixed-capacity  server,  when  the  arrival  processes  are  an  aggregation  of  multiple 
i.i.d.  flows.  We  study  a  particular  scaling  of  the  traffic,  known  as  many-sources  scaling,  where 
the  number  of  flows  scales  with  the  channel  capacity.  The  interested  queuing  performance  is  the 
asymptotic  buffer  overflow  probability.  Assuming  a  many-sources  sample  path  large-deviation 
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principle  (LDP)  for  the  arrival  processes,  we  establish  an  LDP  for  the  queue  length  process  by 
employing  Garcia's  extended  contraction  principle  that  is  applicable  to  quasi-continuous 
mappings.  In  the  future,  we  hope  to  use  this  result  [1]  and  the  result  in  [2]  to  extend  the  work  in 
[6]  to  multiple  access  channels. 


1 


Fig.  6.  Multiple-queue,  single-server  model 


Fig.  7.  Example  of  the  actual  and  the  lower  and  upper  bounds  of  the  exponents  of  the 

1  2 

buffer  overflow  probability  where  the  threshold  is  (b  ,  b  =1),  for  compound  Poisson  bit- arrival 
process  (at  various  arrival  rates)  and  for  two  users. 
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2.5  Fast  Power  Allocation  for  MIMO  Relays 
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Highlights  of  Research 


MIMO  relays  are  useful  building  bloeks  for  tactical  wireless  mobile  ad  hoc  networks.  Power 
allocations  in  space,  time  and  frequency  are  essential  for  MIMO  relays  to  perform  efficiently. 
Fast  computation  of  power  allocation  in  response  to  fast  changing  environment  such  as  channel 
state  information  and  network  topology  is  critically  important  for  network  performance.  Our 
achievements  in  the  past  year  include;  a)  the  proof  of  the  optimality  of  diagonalization  of 
multihop  non-regenerative  MIMO  relays  for  any  Schur  convex  or  Schur  concave  cost  functions, 
b)  the  establishment  of  a  generalized  water  filling  theorem  for  MIMO  link  power  allocation  with 
multiple  power  constraints,  and  c)  the  development  of  a  class  of  power  allocation  algorithms  for 
uplink  and  downlink  of  a  non-regenerative  MIMO  relay  system  with  multiple  users. 
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a)  Optimality  of  Diagonalization 


We  consider  a  multiple-hop  MIMO  relay  system  where  each  relay  performs  waveform 
amplification  and  transformation  without  digital  decoding  and  re-encoding.  Such  a  relay  is  also 
called  non-regenerative  relay.  Comparing  to  regenerative  relays,  non-regenerative  relays  are 
much  simpler  and  have  much  less  end-to-end  delay.  However,  a  disadvantage  of  non- 
regenerative  relays  is  the  effect  of  noise  accumulation.  To  minimize  the  effect  of  noise 
accumulation,  the  design  of  the  relay  matrix  for  each  relay  is  critical. 


Figure  1:  A  multi -hop  MIMO  relay  system 


As  shown  in  Figure  1,  the  output  waveform  vector  of  the  source  node  is  Xj  =  FjS  where  s  is  the 
vector  of  independent  source  symbols  and  Fj  is  the  precoding  matrix  at  the  source.  The 
waveform  vector  received  by  the  first  relay  node  is  yj  =  HjXj  -l-  Vj  where  Hj  is  the  channel 
matrix  of  the  first  link  and  Vj  is  the  noise.  The  relay  matrix  of  the  first  relay  is  F2 ,  i.e., 
X2  =  F2y2 .  Similar  expressions  hold  for  all  other  relays.  At  the  destination,  the  recovered  symbol 
vector  is  s  =  where  is  the  waveform  vector  received  by  the  destination  node  and 

W"  is  the  minimum  mean  square  error  (MMSE)  equalization  matrix.  Here,  we  assume  that 
interferences  between  nodes  that  are  two  or  more  hops  apart  are  negligible. 


One  important  question  is  how  we  should  choose  the  matrices  Fj ,  F2 ,  . . . ,  F^^  where  L  is  any 
positive  integer.  The  answer  to  this  question  in  general  depends  on  the  criterion  of  interest.  Such 
criteria  may  include  the  maximization  of  the  end-to-end  mutual  information  (i.e.,  capacity),  the 
minimization  of  the  averaged  mean  squared  errors  in  s ,  or  the  minimization  of  the  maximum 
mean  square  errors  in  s .  It  turns  out  that  all  these  criteria  belong  to  a  much  larger  class  of 
functions  called  Schur  convex  functions  or  Schur  concave  functions  in  majorization  theory  [1]. 
After  several  attempts  over  several  years,  we  finally  succeeded  to  prove  in  [2]  that  the  optimal 
structures  of  all  the  relay  matrices  under  the  Schur  convex  or  Schur  concave  criterion  subject  to 
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power  constraints  at  all  relays  are  such  that  the  effective  channels  between  s  and  s  are 
diagonalized  and  the  singular  values  of  all  channel  matrices  are  aligned  with  each  other 
according  to  their  values.  This  fundamental  result  is  a  multi-hop  generalization  of  many 
previously  published  results  for  two-hop  MIMO  relays,  e.g.,  see  references  in  [3].  It  is  also  a 
result  that  other  researchers  have  not  been  able  to  establish  in  their  studies  of  multi-hop  MIMO 
relays.  The  optimality  of  diagonalization  implies  that  the  singular  vectors  of  each  of  the 
unknowns  Fj,  F2,  ...,  F^  should  be  constructed  from  the  singular  vectors  of  the  corresponding 
channel  matrix,  and  the  only  remaining  unknowns  are  the  singular  values  of  each  of  Fj,  F2 ,  ..., 
F^ .  This  simplifies  the  computation  tremendously.  The  singular  values  of  Fj ,  F2 ,  . . . ,  F^^  can  be 
found  via  cyclic  water  filling  search. 


b)  Generalized  Water  Filling 


Water  filling  is  a  very  useful  concept  and  algorithm  for  allocating  power  for  a  MIMO  (and/or 
multi-carrier)  link  with  a  conventional  power  constraint.  But  for  tactical  MIMO  networks, 
multiple  power  constraints  may  be  required  for  a  MIMO  link  for  purposes  such  as  controlled 
interferences  to  other  nodes  in  the  network.  This  generalized  problem  can  be  formulated  as 
follows: 


maxQ>o  log|l  +  HQH"|  (1) 
s.t.,  T*  for  i  =  l,2,---,m 

where  H  is  the  channel  matrix  of  the  link  under  consideration,  Q  the  source  covariance  matrix, 
B,.  the  constraint  matrices  depending  on  applications.  For  the  conventional  problem  where 

m  =  \,  the  conventional  water  filling  algorithm  is  the  most  efficient  method.  But  when  m  >  1 , 
the  common  method  to  solve  (1)  is  the  general-purpose  convex  problem  solver  known  as  CVX  in 
Matlab.  In  [4],  we  established  the  following  theorem: 


Theorem  (Generalized  Water  Filling):  The  solution  to  (1)  is  given  by: 

Q  =  K“^V(I-2:“^)+V^K“^  (1) 

where  K  =  (^^^///B^^Bj)^^^  (assumed  to  be  non-singular),  V  and  L  are  determined  from  the 
SVD  ,  (0^  replaces  all  negative  diagonal  elements  by  zeros  and  leaves  all  non- 
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negative  diagonal  elements  unehanged,  and  \i  =  are  the  solution  to  the  following  dual 

problem; 


m 

max-log|I  +  HQH^  | )  -  P,)  (2) 

i=i 

s.t.  Q  =  K“^V(I-2:“^)+VK“^ 


The  dual  problem  is  convex  and  can  be  solved  efficiently  using  the  Newton’s  method. 
Furthermore,  when  the  parameter  space  of  the  dual  problem  is  much  smaller  than  that  of  the 
primal  problem,  which  is  a  typical  case,  the  dual  problem  can  be  solved  with  much  less  time  than 
the  primal  problem.  The  resulting  algorithm  is  much  faster  than  CVX. 


c)  Power  Scheduling  for  a  MIMO  Relay  System  with  Multiple  Users 


Also  shown  in  [4],  we  have  developed  several  power  scheduling  algorithms  for  uplink  as  well  as 
downlink  of  a  MIMO  relay  system  with  multiple  multi-antenna  users.  The  relay  node  performs 
waveform  amplification  and  transformation  without  digital  decoding  and  re-encoding,  which  is 
hence  non-regenerative. 


y. 


y 


K 


Figure  2:  Downlink  of  multiuser  MIMO  relay 


For  downlink,  our  algorithms  search  for  the  optimal  source  covariance  matrix  and  the  optimal 
relay  matrix  to  either  maximize  the  users’  sum  rate  subject  to  a  power  constraint  or  minimize  the 
total  power  consumption  subject  to  users’  rate  constraints. 
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Figure  3:  Uplink  of  multiuser  MIMO  relay 


For  uplink,  our  algorithms  optimize  the  source  covariance  matrices  of  multiple  users  and  the 
relay  matrix.  Several  cost  functions  are  considered  depending  on  the  choices  such  as  dirty  paper 
coding  (DPC)  or  zero  forcing  dirty  paper  coding  (ZFDPC)  at  the  sources.  The  source  covariance 
matrices  and  the  relay  matrix  can  be  searched  by  a  joint  gradient  search  or  a  cyclic  search.  In  all 
cases  except  for  an  algorithm  based  on  high-SNR  approximation  and  geometric  programming, 
the  algorithms  obtain  locally  optimal  results  of  the  original  problems.  To  increase  the  likelihood 
of  finding  the  global  optimal,  we  considered  multiple  random  initializations.  This  strategy  is 
naturally  suitable  for  parallel  computations  with  multi-core  microprocessors.  One  important 
observation  from  [4]  is  that  approximating  a  non-convex  problem  by  a  convex  problem  is  not 
necessarily  the  best  strategy.  Multiple  gradient  searches  with  multiple  random  initializations  can 
achieve  much  better  results  than  a  convex  approximation  based  approach. 


References: 


[1]  A.  W.  Marshall  and  I.  Olkin,  Inequalities:  Theory  of  Majorization  and  Its  Applications . 
Academic  Press,  1979. 

[2]  Y.  Rong  and  Y.  Hua,  “Optimality  of  diagonalization  of  multi-hop  MIMO  relays,”  IEEE 
Transactions  on  Wireless  Communications,  submitted  March  2009,  revised  June  2009. 

[3]  Y.  Rong,  X.  Tang  and  Y.  Hua,  “A  unified  framework  for  optimizing  linear  non- 
regenerative  multicarrier  MIMO  relay  communication  systems,”  IEEE  Transactions  on 
Signal  Processing,  accepted  June  2009. 

[4]  Y.  Yu  and  Y.  Hua  “Power  scheduling  for  a  MIMO  relay  system  with  multiple-antenna 
users”,  IEEE  Transactions  on  Signal  Processing,  submitted  April  2009,  revised  July 
2009. 


59 


2.6  Interference  Cancellations  in  Multiuser  and  Wireless  Relay 
Networks 


PI:  Hamid  Jafarkhani 
Students: 

Erdem  Koyuncu:  Partially  supported  by  MURl  (PhD  student) 

Feng  Li:  Partially  supported  by  MURl  (PhD  student) 

We  have  had  contributions  in  three  research  projects: 

(1)  Interference  cancellation  in  multiuser  networks 

(2)  MIMO  Multi-Hop  Channels 

(3)  Canceling  interference  in  wireless  relay  networks 

In  the  following  we  briefly  summarize  our  progress  in  eaeh  projeet: 

(1)  Interference  cancellation  in  multiuser  networks 


The  issue  of  interference  in  wireless  networks  has  recently  attracted  a  lot  of  attention.  We  had  proposed 
interference  cancellation  methods  for  any  number  of  users  and  any  number  of  antennas  at  the  users  and 
the  receiver.  We  continued  our  contribution  by: 


A,  Performance  Analysis  of  Multiple  Antenna  Multi-User  Detection 


We  derive  the  diversity  order  of  our  previously  proposed  multiple  antenna  multi-user  cancellation  and 
detection  schemes.  The  common  property  of  our  detection  methods  is  the  usage  of  Alamouti  and  quasi- 
orthogonal  space-time  block  codes.  For  detecting  J  users  each  having  N  transmit  antennas,  these  schemes 
require  only  J  antennas  at  the  receiver.  Our  analysis  shows  that  when  having  M  receive  antennas,  the 
array-processing  schemes  provide  the  diversity  order  of  N(M  -  J  -l-  1).  In  addition,  our  results  prove  that 
regardless  of  the  number  of  users  or  receive  antennas,  when  using  maximum-likelihood  decoding  we  get 
the  full  transmit  and  receive  diversities,  i.e.  NM,  similar  to  the  no -interference  scenario.  [33] 
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B,  Multiple-Antenna  Interference  Cancellation  and  Detection  for  Two  Users  Using  Precoders 


Also,  we  consider  interference  cancellation  for  a  system  with  two  users  when  users  know  each  other 
channels.  The  goal  is  to  utilize  multiple  antennas  to  cancel  the  interference  without  sacrificing  the 
diversity  or  the  complexity  of  the  system.  Before,  in  the  literature,  it  was  shown  how  a  receiver  with  2 
receive  antennas  can  completely  cancel  the  interference  of  two  users  and  provide  a  diversity  of  2  for  users 
with  2  transmit  antennas.  We  propose  a  system  to  achieve  the  maximum  possible  diversity  of  4  with  low 
complexity.  Our  main  idea  is  to  design  precoders,  using  the  channel  information,  to  make  it  possible  for 
different  users  to  transmit  over  orthogonal  spaces.  Then,  using  the  orthogonality  of  the  transmitted 
signals,  the  receiver  can  separate  them  and  decode  the  signals  independently.  We  analytically  prove  that 
the  system  provides  full  diversity  to  both  users.  In  addition,  we  provide  simulation  results  that  confirm 
our  analytical  proof  [35],  [36] 


(2)  MIMO  Multi-Hop  Channels 


MIMO  multi-hop  channel  plays  an  important  role  in  wireless  ad  hoc  networks.  In  this  work,  we 
investigate  the  resource  allocation  optimization  problem  in  order  to  achieve  throughput  maximization 
with  given  resource  constraints.  We  propose  algorithms  with  low  complexity  to  achieve  maximum 
capacity.  The  main  idea  is  to  determine  the  rank  of  the  optimal  transmit  covariance  matrix  and  the  optimal 
power  allocation  of  each  node  separately  using  our  low-complexity  algorithm.  Further,  we  reduce  the 
complexity  of  our  algorithm  by  adding  another  pre-processing  algorithm.  Using  our  algorithms,  we  find 
that  while  dynamical  allocation  of  time  and  power  could  increase  channel  capacity,  equal  time  and  power 
allocation  among  different  nodes  may  not  cause  much  capacity  loss.  [32] 


(3)  Canceling  interference  in  wireless  relay  networks 


In  another  effort,  we  consider  the  problem  of  interference  in  wireless  relay  networks.  There  have  been 
many  cooperative  schemes  proposed  for  relay  networks.  Most  of  these  schemes  are  designed  for  a  single 
user.  We  consider  a  relay  network  that  allows  multiple  users  transmit  simultaneously.  Our  main 
contributions  are: 


A,  Interference  Cancellation  in  Distributed  Space-Time  Coded  Wireless  Relay  Networks 


We  consider  the  interference  cancellation  (IC)  problem  in  multi-user  wireless  relay  networks.  First,  it  is 
shown  that  using  distributed  space-time  coding  (DSTC),  the  multiple  antenna  1C  scheme  that  we  had 
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proposed  for  systems  with  direct  transmissions  can  be  applied  to  relay  networks.  The  ML  decoding  after 
full  1C  can  be  performed  symbol  by  symbol.  Then,  by  allowing  1C  at  relays,  a  new  degree  of  freedom  in 
relay  network  design  is  discovered.  With  this  new  idea,  the  required  number  of  antennas  at  the  receiver 
for  full  1C  can  be  reduced  and  a  balance  between  diversity  and  delay  can  be  obtained.  [34] 


B.  Beamforming  in  Wireless  Relay-Interference  Networks  with  Quantized  Feedback 


We  consider  quantized  beamforming  in  wireless  amplify-and-forward  (AF)  relay  networks.  We  use  the 
Generalized  Lloyd  Algorithm  (GLA)  to  design  the  quantizer  of  the  feedback  information  and  specifically 
to  optimize  the  bit  error  rate  (BER)  performance  of  the  system.  Achievable  bounds  for  different 
performance  measures  are  derived.  First,  we  analytically  show  that  a  simple  feedback  scheme  based  on 
relay  selection  can  achieve  full  diversity.  Unlike  the  previous  diversity  analysis  on  the  relay  selection 
scheme,  our  analysis  is  not  aided  by  any  approximations  or  modified  forwarding  schemes.  Then,  for  high- 
rate  feedback,  we  find  an  upper  bound  on  the  average  signal-to-noise  ratio  (SNR)  loss  and  show  that  it 
decays  at  least  exponentially  with  the  number  of  feedback  bits,  B.  Using  this  result,  we  also  demonstrate 
that  the  capacity  loss  also  decays  at  least  exponentially  with  B.  In  addition,  we  provide  approximate  upper 
and  lower  bounds  on  the  BER,  which  can  be  calculated  numerically.  Simulations  are  also  provided,  which 
confirm  our  analytical  results.  We  observe  that,  for  R  relays,  our  designs  achieve  full  diversity  when  B  > 
logR  and  a  few  extra  feedback  bits  are  sufficient  for  a  satisfactory  performance  in  terms  of  the  array  gain. 
Simulations  also  show  that  our  approximate  BER  is  a  reliable  estimation  on  the  actual  BER  for  even 
moderate  values  of  B.  [37] 
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2.7  Feedback  MIMO  Systems 
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SUMMARY  OF  RESEARCH 


Our  work  is  concerned  with  feedback  in  MIMO  based  Ad  Hoc  networks;  development  of 
effective  quantization  methods,  analysis  of  feedbaek  systems  with  finite  rate  feedbaek,  analysis 
of  feedback  systems  with  imperfect  channel  state  information,  and  developing  protocols  that 
support  multiple  simultaneous  transmissions. 


a)  Average  SEP/BEP  analysis  with  Channel  Estimation  Errors  and  Feedback  Delay 


In  this  work,  we  analyze  average  symbol  error  probability  and  bit  error  probability  performances 
of  M-PSK  and  M-ary  rectangular  QAM  constellations  with  Gray  code  mapping. 


This  was  the  work  was  mainly  earried  out  in  the  first  two  years  of  MURI  projeet  and  the  journal 
paper  related  to  this  work  appeared  in  IEEE  transactions  on  communications  Jan.  2009  issue  [1]. 
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b)  Optimum  Codebook  Design  for  Minimizing  the  Average  SEP  Loss  and  Analysis  of  Loss 
Due  to  Channel  Quantization 


In  this  work,  we  focused  on  Multiple  Input  Single  Output  (MISO)  systems  where  channel  state 
information  (CSI)  is  conveyed  from  the  receiver  to  the  transmitter  through  a  finite -rate  feedback 
link.  The  importance  of  the  choice  of  performance  metric  and  the  effect  of  mismatch  in  the 
channel  statistics  assumptions  are  the  main  focus  of  this  work. 

Some  of  the  ideas  related  to  this  work  were  discussed  in  last  year’s  report,  and  the  journal  paper 
with  complete  details,  appeared  in  IEEE  transactions  on  signal  processing  May  2009  issue  [2], 


c)  Modeling  of  Feedback  Imperfections:  Distinction  between  Feedback  Delay  and 
Estimation  Error 

This  work  was  initiated  last  year  and  constitutes  a  major  portion  of  our  work  this  year.  The 
motivation  for  this  work  is  to  address  some  of  the  shortcomings  of  existing  work.  In  particular,  to 
develop  models  of  feedback  systems  that  more  closely  reflect  real  communication  systems  and  to 
provide  analytical  results  for  appropriate  and  meaningful  metrics.  We  made  significant 
improvements  on  the  modeling  of  imperfections  and  identified  the  proper  metric  namely  Packer 
Error  Probability  (PEP).  On  the  analytical  side,  we  have  clearly  pushed  the  boundary  and 
developed  new  tools.  The  expressions  are  complex  but  simple  enough  to  compute  and  so  are  still 
useful  in  that  they  can  obviate  the  need  for  computer  simulations  that  are  time  consuming.  A 
journal  paper  related  to  this  work  was  submitted  to  IEEE  transactions  on  signal  processing  in 
May  2009  and  revised  in  Jul.  2009  [3].  ].  A  conference  paper  related  to  this  work  also  appeared 
in  Globecom’08  [4].  Another  journal  paper  related  to  developing  an  analytically  tractable 
approximation  for  the  Gaussian  Q-function  also  appeared  in  IEEE  communications  letters  Sep. 
2008  issue  [5].  We  now  briefly  discuss  the  main  ideas  related  to  this  work. 

Multiplicative  fading  is  a  major  source  of  performance  degradation  in  multipath  wireless 
environment.  Channel  coding  and  interleaving  can  offer  some  protection  from  the  negative 
effects  of  fading.  However,  in  some  wireless  systems  data  has  to  be  organized  into  small 
packets,  which  are  confined  to  fixed  time  slots,  with  or  without  interleaving.  One  popular 
example  of  such  a  system  is  the  slotted  multiple  access  scheme.  It  is  important  for  the  system 
designers  to  know  the  impact  of  fading  on  the  performance.  An  important  metric  for  studying  the 
performance  of  a  non-interleaved  wireless  packet  data  transmission  is  the  average  packet  error 
probability  (PEP).  Packet  error  probability  is  also  increasingly  becoming  an  important  quality-of- 
service  parameter  for  the  wireless  networking  community  since  it  determines  how  frequently  the 
information  packet  has  to  be  re-transmitted. 
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Extensive  analytical  results  quantifying  the  impact  of  fading  on  average  symbol  and  error 
probability  (SEP/BEP)  are  available  for  various  modulation  schemes.  However,  in  slow  fading 
situations,  there  is  no  mapping  between  the  average  SEP/BEP  and  the  average  PEP. 
Consequently  knowing  average  SEP/BEP  does  not  help  in  understanding  the  average  PEP. 
Analysis  of  average  PEP  is  a  more  complicated  problem  compared  to  the  analysis  of  average 
SEP/BEP.  Analytical  quantification  of  packet  error  probability  has  received  considerable 
attention  in  the  literature.  Closed-form  expressions  for  PEP  have  been  derived  for  the  non¬ 
coherent  ESK  modulation.  The  non-coherent  ESK's  SEP,  conditioned  on  the  channel,  is  an 
exponential  function  and  taking  expectation  of  the  higher  powers  of  conditional  SEP  w.r.t.  the 
fading  random  variable  is  analytically  tractable.  However,  closed-form  expressions  are  not 
available  for  coherent  BPSK  and  other  constellations.  Conditional  PEP  (conditioned  on  a 
function  of  the  wireless  channel)  for  a  scheme  such  as  coherent  BPSK  results  in  integer  powers 
of  the  Gaussian-Q  function.  This  makes  the  analysis  challenging  because  in  order  to  derive  the 
average  PEP  expression,  one  has  to  integrate  the  integer  powers  of  the  Gaussian-Q  function  w.r.t. 
the  random  variable  that  captures  the  fading  environment,  an  analytically  difficult  exercise.  We 
also  note  that,  to  the  best  of  our  knowledge,  the  effect  of  channel  estimation  errors  on  PEP  has 
not  been  considered  in  the  literature. 


In  this  work  we  consider  the  problem  of  deriving  analytical  expressions  for  PEP  of  a  multiple 
input  single  output  system  with  various  forms  of  practical  imperfections.  We  later  show  that  this 
problem  captures  various  commonly  interested  performance  analysis  of  wireless  systems  as 
special  cases.  The  first  form  of  feedback  imperfection  considered  is  channel  estimation  error.  It 
is  now  a  common  practice  to  model  the  actual  channel  and  its  estimate  as  a  jointly  Gaussian 
random  process,  with  an  error  term  that  is  orthogonal  to  the  channel  estimate.  The  error  term 
associated  with  a  particular  channel  estimate  is  unknown  to  the  receiver  and  hence  it  becomes 
part  of  noise  when  the  performance  analysis  is  carried  out.  If  the  channel  under  consideration  is 
varying  at  symbol  level,  or  if  the  performance  criteria  is  average  symboEbit  error  probability, 
then  the  variance  of  the  error  term  will  be  simply  added  (along  with  the  symbol  dependency)  to 
the  variance  of  the  receiver  noise  resulting  in  an  effective  noise  term  with  variance  equaling  the 
sum  of  variance  of  receiver  noise  and  the  variance  of  the  estimation  error  term.  In  this  part  of  the 
work,  we  also  follow  the  standard  model  of  joint  Gaussianity  between  the  channel  and  its 
estimate,  but  adapt  it  to  the  packet  fading  model.  An  important  difference  is  that  in  a  packet 
fading  model  the  error  term  is  constant  for  the  entire  packet  while  each  symbol  experiences  a 
different  noise  sample  requiring  new  analytical  tools. 


The  second  form  of  feedback  imperfection  we  address  is  the  delay  between  constructing  the 
beamforming  vector  at  the  receiver  and  using  it  at  the  transmitter.  Another  well  accepted 
formulation  is  to  treat  the  impact  of  feedback  delay  in  a  manner  similar  to  estimation  errors,  i.e.. 


65 


actual  channel  and  its  delayed  version  are  assumed  to  be  jointly  Gaussian  with  an  unknown  (to 
the  reeeiver)  error  term  that  is  orthogonal  to  the  delayed  version.  Since  the  delay  related  error 
term  is  unknown  to  the  receiver,  similar  to  estimation  related  error  term,  during  performance 
analysis  it  becomes  part  of  noise  thus  removing  any  coneeptual  distinction  between  the  mismatch 
in  beamforming  due  to  feedbaek  delay  and  estimation  errors.  Though  mueh  of  the  past  work  on 
feedback  delay  effectively  make  the  delay  related  error  term  part  of  receiver  noise,  alternate 
options  were  considered  (primarily  in  the  context  of  adaptive  modulation).  However,  it  is 
important  to  note  that  much  of  the  work  treated  estimation  errors  and  feedback  delay  in  a  similar 
manner,  i.e.,  either  both  the  error  terms  are  assumed  to  be  known  or  unknown  to  the  reeeiver.  In 
this  work,  based  on  feedback  system  considerations  we  feel  it  is  appropriate  to  treat  the  errors 
due  to  feedbaek  delay  to  be  known  at  the  receiver,  while  the  errors  due  to  estimation  errors  are 
un-known  at  the  receiver.  This  modeling  approaeh  is  adopted  in  this  work  and  it  shows  that  the 
impaet  of  feedback  delay  on  beamforming  MISO  system  performance  can  be  less  severe  and  is 
also  eonceptually  quite  different  from  ehannel  estimation  errors. 


The  third  form  of  feedback  imperfection  considered  in  this  work  is  finite -rate  quantization  of  the 
channel.  To  summarize,  the  contributions  presented  in  this  work  are  threefold:  an  accurate 
eharacterization  of  estimation  errors  in  a  paeket  fading  context,  a  new  modeling  of  feedback 
delay  whieh  shows  improved  performance  for  a  beamforming  MISO  system  and  conceptually 
distinguishes  it  from  estimation  errors,  and  derivation  of  an  analytieal  expression  quantifying  the 
impact  of  channel  estimation  errors,  feedbaek  delay  and  channel  quantization  on  the  average 
paeket  error  probability.  All  these  eontributions  further  the  understanding  of  feedback 
communication  systems.  As  a  side  benefit,  the  analytieal  tools  developed  promise  to  be  of 
general  interest  with  broad  applicability. 


d)  Wireless  Channel  Modeling  and  Prediction 


In  this  work,  we  turn  our  attention  to  the  issue  of  minimizing  the  negative  impact  of  feedback 
delay  on  the  performanee  of  the  system.  An  obvious  solution  to  the  problem  of  feedback  delay  is 
to  predict  the  channel  and  then  quantize  the  predicted  channel.  In  this  chapter  we  study  the  role 
of  ergodieity  in  wireless  ehannel  prediction.  With  an  eye  on  developing  a  better  wireless  channel 
simulator  we  first  begin  with  a  general,  but  non-linear  and  non-traetable,  form  of  wireless 
ehannel  model.  Under  eertain  assumptions  we  then  consider  a  simplified  and  well  aecepted  linear 
sinusoidal  ehannel  model.  Following  the  sinusoidal  ehannel  model,  eonditions  under  whieh  the 
ergodic  assumption  is  valid  are  presented.  This  sheds  insight  into  when  statistical  channel  models 
that  employ  ensemble  averaging  are  appropriate. 
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Due  to  the  laek  of  ergodieity  in  a  typieal  real  world  wireless  ehannel,  Least  Squares  predietion, 
an  approaeh  based  on  time  averages  is  motivated  as  opposed  to  linear  minimum  mean  squared 
error  ehannel  predietion,  an  approaeh  based  on  ensemble  averaging.  We  then  study  methods  sueh 
as  Forward-Backward  and  rank  reduction  for  high  quality  channel  prediction.  Simulation  results 
are  presented  to  complement  the  analytical  expressions.  Simulation  results  also  illustrate  the 
improvement  in  channel  prediction  quality  [6]. 


e)  MIMO  Ad  Hoc  Networks 

Last  year  we  proposed  a  novel  asynchronous  Media  Access  Control  (MAC)  protocol, 
Opportunistic  MAC  (OMAC),  for  MIMO  ad-hoc  networks  [10].  The  proposed  solution  is  based 
on  closed  loop  minimal  feedback  antenna  selection  diversity  scheme  and  optimums  receive 
combining.  The  use  of  antenna  selection  diversity  contributes  to  a  reduction  in  the  feedback 
information  and  in  the  effective  interference  produced.  To  utilize  the  spatial  degrees  of  freedom 
offered  by  MIMO,  we  propose  the  use  of  a  novel  rank  based  metric  to  obtain  interference 
information  as  well  to  to  enable  multiple  simultaneous  transmissions  and  to  make  MAC 
decisions.  The  rank  of  the  interference  matrix,  Rj  is  used  as  a  metric.  Through  analysis  and 
simulation,  we  found  that  the  proposed  protocol  significantly  outperforms  802.1 1  MIMO  and 
obtained  high  spatial  degree  of  freedom  utilization. 


The  performance  evaluation  in  our  work  had  assumed  that  the  rank  estimation  mechanism  was 
error-free  (i.e.,  it  would  always  correctly  estimate  the  number  of  concurrent  communications).  In 
our  current  work  [7],  we  develop  a  robust  method  to  estimate  this  quantity,  based  on  a  rank 
estimation  algorithm  resilient  to  noise  and  interference.  Our  key  contribution  is  the  design  and 
performance  evaluation  of  a  novel  carrier  sense  scheme  tailored  to  MIMO  that  exploits  the 
spatial  structure  of  MIMO  waveforms.  The  proposed  carrier  sense  system  used  here  is  by  no 
means  specific  to  OMAC,  and  could  also  be  adopted  by  other  MIMO  ad  hoc  protocols.  The 

other  MIMO  MAC  protocols  based  on  carrier  sense  may  use  the  conventional  approach  (i.e., 
compare  the  average  received  power  across  the  antennas  against  a  certain  threshold)  that  does 
not  exploit  this  structure  or  may  propose  ad  hoc  techniques  that  work  only  when  transmissions 
are  slot  synchronous,  while  our  scheme  works  also  for  asynchronous  communications.  We  show 
that  this  carrier  sense  mechanism  can  help  realize  the  communication  parallelism  inherent  in 
MIMO  at  a  limited  cost  in  computational  complexity. 
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2.8  Hierarchical  modulation  for  jointly  combined  space-time  coding 
and  spatial  multiplexing  in  MIMO  system 


PI:  Laurence  Milstein 
GSR:  Seok-Ho  Chang 
1.  Introduction 

We  are  interested  in  the  cross-layer  design  of  a  tactical  communications  system  to  support  real-time  video 
transmissions.  The  basic  element  of  the  design  is  the  use  of  MIMO  technology  to  achieve  unequal  error 
protection  (UEP)  among  the  different  components  of  a  scalable  video  waveform.  As  will  be  discussed 
below,  scalable  video  signals  are  designed  so  that  different  parts  of  the  waveform  have  different  levels  of 
importance.  Thus,  UEP  is  a  natural  way  of  protecting  the  more  important  components  more  heavily  than 
the  less  important  components.  The  goal  of  the  research  is  to  jointly  use  the  key  properties  of  the  physical 
layer  and  application  layer  to  enhance  system  performance  through  the  use  of  MIMO  designs. 


Scalable  video  has  both  a  base  layer  and  an  enhancement  layer.  Successful  reception  of  only  the  base 
layer  provides  relatively  poor  reproduction  of  the  video,  but  the  quality  improves  monotonically  as  more 
of  the  enhancement  layer  is  successfully  received.  However,  if  the  base  layer  is  not  correctly  received,  the 
enhancement  layer  is  useless.  Thus,  the  base  layer  should  be  protected  more  strongly  than  the 
enhancement  layer. 


Since  the  base  layer  of  scalable  video  should  be  transmitted  more  reliably  than  the  enhancement  layer, 
space-time  coding  can  be  employed  for  the  base  layer,  while  spatial  multiplexing  can  be  used  for  the 
enhancement  layer.  In  this  research,  we  combine  the  UEP  that  arises  with  the  use  of  a  hierarchical 
constellation,  with  the  UEP  that  can  be  achieved  through  the  use  of  MIMO. 


11.  Proposed  scheme 


We  regard  the  hierarchical  constellation  as  the  superposition  of  two  subconstellations,  a  basic 
subconstellation  and  a  secondary  subconstellation,  where  the  Euclidian  distance  between  the  signal  points 
of  the  basic  subconstellation  is  greater  than  that  of  the  secondary  subconstellation.  As  an  example,  for  2 
by  2  MIMO  systems,  the  constellation  symbols  of  our  scheme  are  given  by 
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(1) 


5^[2n]  5^[2«  +  l] 
SB[2n]  Ssiln  +  l] 


Xi[2n]  +  X2A[2n]  -x*[2n  +  l]  +  X2^[2n  +  1] 

Xj  [2n  + 1]  +  X2b  [2«]  Xj  [2n]  +  X25  [2n  + 1] 


where  eaeh  row  corresponds  to  a  transmit  antenna,  each  column  corresponds  to  a  time  symbol,  and  other 
notation  is  defined  as  follows: 


S ^\n\ ,  Sg\n\ :  hierarchical  constellation  symbols  which  are  transmitted  on  antennas  A  and  B, 
respectively. 

x,[«] :  basic  subconstellation  symbol  which  is  transmitted  using  an  Alamouti  code. 

^2A  ’  ^2B  ■  secondary  subconstellation  symbols  which  are  transmitted  using  spatial  multiplexing  on 

antennas  A  and  B,  respectively. 


Note  that  the  basic  and  secondary  subconstellations  can  have  an  arbitrary  alphabet  size.  As  an  example. 
Fig.  1  (a)  uses  QPSK  for  the  basic  subconstellation,  where  the  minimum  Euclidian  distance  between 

signal  points  is  given  by  dj  (note  that  the  average  power,  12).  Fig.  1  (b)  shows  the  secondary 

subconstellation  using  16  QAM  whose  minimum  Euclidian  distance  is  given  by  J2  (note  that  the  average 

power,  P2  =  5^2  12).  Eastly,  Fig.  2  shows  hierarchical  4/64  QAM  which  is  the  superposition  of  the  basic 
and  the  secondary  subconstellations  shown  in  Figs.  1  (a)  and  (b). 


Note  that  from  (1),  we  are  effectively  using  fourth  order  diversity  for  the  two  base  layer  bits  of  each  64 
QAM  symbol,  and  using  spatial  multiplexing  for  the  four  enhancement  layer  bits  of  each  64  QAM 
symbol. 
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(a)  QPSK  (Almouti  coding) 


(b)  16  QAM  (Spatial  multiplexing) 


Fig,  1,  Basic  and  secondary  subconstellations. 
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Fig.  2.  Hierarchical  4/64  QAM. 


111.  Simulation  results 


We  evaluate  the  performance  of  the  proposed  hierarchical  4/64  QAM  scheme  using  the  mapping  of  (1). 
Note  that  for  each  symbol  duration,  the  proposed  scheme  transmits  2  bits  and  8  bits  for  the  base  layer  and 
the  enhancement  layer,  respectively.  For  comparison  purposes,  we  also  evaluate  the  performance  of  a 
hierarchical  2/32  QAM  scheme  which  employs  just  spatial  multiplexing,  and  a  hierarchical  4/1024  QAM 
scheme  which  employs  just  space-time  coding.  Note  that  both  schemes  also  transmit  2  bits  and  8  bits  for 
the  base  and  enhancement  layers,  respectively,  for  each  symbol  duration.  In  other  words,  the  above  three 
schemes  all  have  the  same  data  rate. 


We  use  the  progressive  source  coder  SPIHT  [2]  as  an  example,  and  provide  the  results  for  the  standard  8 
bits  per  pixel  (bpp)  512x512  Lena  image  [3]  with  a  transmission  rate  of  0.375  bpp.  To  compare  the 
image  quality,  we  use  peak-signal-to-noise  ratio  (PSNR),  defined  as 


PSNR  =  lOlog 


255^ 

e[d] 


(dB) 


(2) 


where  255  is  due  to  the  8  bpp  image,  and  E\p\  is  the  expected  distortion  of  the  image.  For  this  analysis, 

we  used  mean-squared  error  as  the  distortion  measure,  D  .  All  the  three  schemes  are  decoded  using 
maximum  likelihood  decoding  [1].  Fig.  3  shows  the  PSNR  performance.  It  is  seen  that  the  proposed 
hierarchical  modulation  which  combines  space-time  coding  and  spatial  multiplexing  outperforms  the  two 
conventional  schemes. 
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SNR  per  symbol  (dB) 

Fig.  3.  PSNR  performance  of  the  proposed  scheme. 


111.  Conclusion  and  future  research 


In  this  report,  we  proposed  a  specific  scheme  for  hierarchical  modulation  which  combines  two  different 
MIMO  techniques  to  minimize  the  expected  distortion  for  either  progressive  images  or  scalable  video  by 
employing  uniform  QAM  constellations  for  the  basic  and  the  secondary  subconstellations.  However, 
subconstellations  themselves  can  be  hierarchical  QAM  constellations.  As  an  example,  for  hierarchical  256 
QAM,  the  basic  subconstellation  can  be  4/16  QAM  and  the  secondary  subconstellation  also  can  be  4/16 
QAM.  Further,  when  we  employ  hierarchical  constellations  for  and  X2g\_n\  of  (1),  the  secondary 

subconstellation  symbols  which  are  transmitted  using  spatial  multiplexing  on  antennas  A  and  B  can  have 
distinct  minimum  Euclidian  distances.  In  [2],  specific  methods  of  multiplexing  hierarchical  constellations 
which  have  distinct  minimum  Euclidian  distances  was  proposed  to  achieve  multiple  levels  of  UEP.  We 
are  now  considering  using  the  techniques  of  [2]  in  conjunction  with  spatial  multiplexing.  As  an  example. 
Fig.  4  shows  hierarchical  4/16  QAM  constellations  having  different  minimum  Euclidian  distances  for  the 
secondary  subconstellations.  Note  that  Alamouti  coding  is  still  used  for  the  base  layer  of  the  scalable 
video. 
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(a)  Secondary  subconstellation  on  antenna  A  (b)  Secondary  subconstellation  on  antenna  B 
Fig.  4.  Hierarchical  4/16  QAM  constellations  having  different  minimum 
Euclidian  distances  for  secondary  subconstellations. 
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2.  Introduction 


In  one  aspeet  of  our  work,  we  eonsider  a  network  composed  of  an  infinite  number  of  routes  on 
an  infinite  plane.  The  first  approaeh  we  take  is  to  assume  that  eaeh  route  simply  eonsists  of  only 
one  link,  i.e.,  a  transmitter  (TX)  and  reeeiver  (RX)  at  distanee  R.  This  single-hop  model  ean  be 
eonsidered  as  a  snapshot  of  an  aetual  multi-hop  network  and  allows  us  to  evaluate  the 
performanee  in  terms  of  single-hop  metrics,  that  refleet  the  performanee  benefit  at  the  multi -hop 
level.  Sueh  metries  are  the  network  throughput,  defined  as  the  product  (spatial  density  of  TXs)  x 
(link  throughput),  whieh  addresses  the  need  to  paek  as  many  transmissions  as  possible  in  spaee; 
and  the  information  efficiency,  defined  as  the  produet  (transmission  distanee)  x  (link  throughput), 
whieh  eaptures  the  trade-off  present  in  a  multi-hop  network  where,  transmitting  farther  means  a 
paeket  needs  fewer  hops  to  reach  its  final  destination.  However,  for  a  fixed  TX  power,  the  RX 
signal-to-interferenee  (SIR)  ratio  is  lower.  The  single-hop  approach  allows  the  analytieal 
evaluation  of  the  aforementioned  metries,  as  a  funetion  of  various  physieal  layer  parameters, 
sueh  as  the  multiple-aeeess  (MA)  scheme,  e.g.  frequeney  hopping  (FH),  the  eoding  seheme,  e.g., 
eonvolutional  eoding,  and,  in  case  the  nodes  have  more  than  one  antenna,  the  multiple-input 
multiple-output  (MIMO)  teehnique,  e.g.,  spaee-time  eoding  or  spatial-multiplexing 


75 


Although  it  provides  useful  insight  on  how  physical  layer  choices  affect  the  link  performance, 
the  above  approach  has  its  shortcomings.  It  implicitly  assumes  that  the  source  of  communication 
and  its  respective  final  destination  lie  at  an  infinite  distance  from  each  other  and  focuses  on 
optimizing  the  performance  of  a  single  link,  with  the  hope  of  providing  a  benefit  at  the  end-to- 
end  level.  Since  the  distance  of  the  final  destination  or  the  number  of  hops  are  not  specified, 
there  are  no  guarantees  in  terms  of  end-to-end  delay  and  throughput.  The  second  part  of  our 
work  assumes  a  simplified  -  but  still  realistic  -  physical  layer  and  focuses  on  evaluating  the  end- 
to-end  delay  and  throughput  over  a  typical  network  route,  when  the  source  and  final  destination 
are  at  a  specified  distance  R,  and  a  number  of  relays  are  placed  on  the  line  between  them,  in 
order  to  forward  the  packets  that  originate  at  the  source.  A  simple  link-layer  protocol  is 
considered,  where,  if  a  packet  is  not  received  correctly  by  a  node,  it  is  stored  at  the  head  of  the 
queue  of  the  previous  node  in  the  route  and  retransmitted  at  the  next  available  opportunity. 
Central  to  our  analysis  is  the  notion  of  stability  of  the  relay  queues,  i.e.,  making  sure  that  their 
lengths  do  not  become  unbounded  over  time,  and,  as  a  result,  the  end-to-end  delay  becomes 
infinite.  We  address  issues  such  as:  determining  the  optimal  number  of  relays  and  their 
placements,  so  that  the  delay  is  minimized  and/or  the  throughput  is  maximized;  the  impact  of 
imperfect  relay  placement  on  the  delay;  the  maximum  allowable  MAC  probability  for 
backlogged  sources  and  the  maximum  allowable  packet  arrival  probability  for  non-backlogged 
sources  such  that  all  queues  in  the  network  are  stable. 

With  regard  to  the  network  topology,  we  assume  that  the  locations  of  the  sources  (TXs  in  the 
first  model)  are  drawn  independently  according  to  a  spatially  homogeneous  Poisson  process  of 
density  X  and  the  orientation  of  the  destinations  (RXs  in  the  first  model)  is  random.  The 
topology,  thus  the  interference  experienced  over  a  link,  can  actually  change  due  to  mobility  or, 
effectively,  due  to  random  access,  if  the  network  is  static.  In  both  cases,  a  link  is  characterized  by 
its  packet  success  probability,  which  is  an  average  measure  of  performance,  over  different 
network  topologies  and  channel  realizations.  We  assume  that  the  channel  between  any  two  nodes 
at  distance  r  includes  Rayleigh  fading  and  path-loss  according  to  the  law,  f'°,  where  b>2  is  the 
path-loss  exponent.  All  nodes  have  the  same  transmit  power  -  normalized  to  one  -  and  additive 
Gaussian  noise  is  disregarded,  i.e.,  we  are  considering  an  interference-limited  environment. 

3.  Summary  of  Achievements. 

In  the  journal  paper  [1],  we  explore  the  benefits  of  channel  diversity  in  wireless  ad  hoc  networks. 
Our  model  is  that  of  a  Poisson  point  process  of  transmitters,  each  with  a  receiver  at  a  given 
distance.  A  packet  is  divided  into  blocks  which  are  transmitted  over  different  subbands  that  are 
determined  by  random  frequency  hopping  (FH).  At  the  receiver,  a  maximum-likelihood  decoder 
is  employed  to  estimate  the  transmitted  packet/codeword. 

We  demonstrate  that  large  gains  in  terms  of  network  capacity  are  possible,  by  combining  FH 
during  packet  transmission  and  error  correction  coding  of  modest  complexity.  If  L  is  the 
Hamming  distance  of  the  convolutional  code  employed  at  the  TX,  X  is  the  density  of  TXs  and  M 
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is  the  number  of  subbands,  we  show  that,  asXIM  0,  the  eodeword  error  probability  seales  as 
O  ( (^  /  M)^).  This  implies  that  the  transmission  eapaeity  seales  as 

s  for  s  ^  0,  where  s  is  the  eonstraint  plaeed  on  the  eodeword  error  probability.  The  pre- 
eonstant  depends  on  the  geometry  of  the  symbol  eonstellation  and  is  also  approximately 
proportional  ioN'^  where  N  is  the  number  of  antennas  at  the  RX  and 

a  =2/b,  where  b  >  2,  is  the  propagation  exponent.  We  also  derive  upper  and  lower  bounds  on  the 
ergodie  eapaeity  C  of  the  typieal  TX-RX  link. 

Praetieal  physieal  layer  issues  are  diseussed  sueh  as  ehannel  estimation,  power  eontrol  (PC)  and 
ehannel  eorrelation.  We  demonstrate  via  simulation  that,  with  an  aeeeptable  rate  loss,  aeeurate 
ehannel  state  information  (CSI)  ean  be  obtained  for  deeoding.  With  respeet  to  PC,  it  is  shown 
that  ehannel  inversion  ean  aetually  improve  the  performanee,  sinee  the  error  eorreetion  eode 
protects  the  RX  from  the  deep  fades  of  its  nearby  interferers.  Finally,  the  impact  of  the  channel 
correlation  is  evaluated  as  the  number  of  subbands  and/or  the  number  of  dwells  is  decreased  and 
it  is  shown  that  the  gains  compared  to  slow  FH  are  still  significant. 

In  journal  paper  [2],  we  investigate  the  end-to-end  delay  performance  of  multi -hop  random 
networks.  We  consider  networks  where  each  route  consists  of  a  source,  a  number  of  relays  and  a 
destination  at  a  finite  distance,  and  the  locations  of  the  sources  and  relays  are  determined 
according  to  independent  Poisson  point  processes.  Nodes  are  equipped  with  queues  to 
accommodate  the  randomness  in  the  delivery  of  a  packet,  which  is  determined  by  the  level  of  the 
SIR  over  each  hop.  Given  a  TDMA/ ALOHA  MAC  protocol,  we  first  study  an  idealized  network 
model  where  all  routes  have  the  same  number  of  hops,  the  same  distance  per  hop  and  their  own 
dedicated  relays.  Assuming  that  a  stationary  regime  exists  for  this  network,  we  analytically 
evaluate  the  mean  end-to-end  delay  and  the  throughput  for  backlogged  and  non-backlogged 
sources.  A  key  point  in  our  analysis  is  taking  into  account  that  the  interference  level  over  each 
hop,  i.e.,  the  density  of  transmitting  nodes,  depends  on  the  packet  success  probabilities  and  vice 
versa.  The  derivation  of  closed-form  expressions  permits  the  optimization  of  the  mean  end-to- 
end  delay  with  respect  to  the  number  and  placement  of  the  relays.  The  benefits  of  this  relay 
selection  -  or  routing  -  strategy  are  then  verified  in  the  original  network  via  simulation. 

If  the  sources  are  backlogged,  we  find  that  the  delay  is  minimized  if  the  first  hop  is  much  larger 
than  the  remaining  hops,  e.g.,  in  a  delay-optimized  three-hop  route,  the  first  hop  covers  half  the 
total  distance.  On  the  other  hand,  in  the  non-backlogged  case,  the  delay  is  minimized  for 
equidistant  hops.  We  demonstrate  that  the  respective  optimal  numbers  of  hops  scale  sublinearly 
and  linearly  with  the  source-destination  distance.  We  also  discuss  stability  issues  in  a  multi-hop 
random  network  and  derive  sufficient  conditions  on  the  medium  access  probability  (MAP),  and 
the  traffic  load  in  the  case  of  non-backlogged  sources,  for  the  rate  stability  of  the  network. 

Simulations  show  that,  when  the  relays  of  each  route  are  selected  out  of  a  random  population,  the 
analytical  results  are  accurate,  provided  that  the  density  of  the  relay  process  is  sufficiently  large. 
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2.  Introduction 

In  this  work  we  have  investigated  interferenee  mitigation  teehniques  in  point-to-multipoint 
(broadeast)  MIMO  eommunieation  systems,  in  whieh  the  interferenee  mitigation  is  performed  at 
the  transmitter.  To  mitigate  the  interferenee  at  the  transmitter,  the  transmitter  must  know  the 
ehannel  eharaeteristies,  typieally  the  ehannel  impulse  response.  This  ehannel  state  information 
(CSI)  may  be  obtained  from  ehannel  measurements  performed  at  eaeh  of  the  reeeivers  by  means 
of  reeeived  pilot  signals  sent  by  the  transmitter.  Then,  the  CSI  must  be  sent  to  the  transmitter.  In 
sueh  a  scenario,  the  ehannel  time  variations  must  be  relatively  slow  so  that  a  reliable  estimate  of 
the  ehannel  eharaeteristies  is  available  at  the  transmitter.  In  some  systems,  the  uplink  and 
downlink  channels  are  identical,  e.g.,  the  same  frequency  band  is  employed  for  both  the  uplink 
and  the  downlink,  but  separate  time  slots  are  used  for  transmission.  This  transmission  mode  is 
ealled  time-division  duplex  (TDD).  In  TDD  systems,  the  pilot  signals  for  ehannel  measurement 
may  be  sent  by  eaeh  of  the  users  in  the  uplink.  In  our  treatment,  we  have  eonsidered  two  eases. 
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In  the  first  case,  we  assume  that  the  channel  characteristics  are  known  perfectly  at  the 
transmitter.  In  the  second  case,  we  assume  that  the  channel  measurement  is  corrupted  by  additive 
noise,  and  we  evaluate  the  effect  of  the  noisy  estimate  on  the  error  rate  performance. 

3.  Summary  of  Achievements 

In  journal  paper  [1],  we  investigate  the  effects  on  the  performance,  degradation  due  to  channel 
estimation  errors  in  a  point-to-multipoint  (downlink)  MIMO  communication  system  in  which  the 
interference  mitigation  is  performed  at  the  transmitter.  In  particular,  we  consider  a  MIMO 
system  which  employs  a  Tomlinson-Harashima  precoder  to  pre-process  the  transmitted  signal  so 
as  to  suppress  the  interference  at  each  of  the  decentralized  receivers.  The  QR  decomposition  of 
the  channel  matrix  is  employed  to  obtain  an  equivalent  channel  in  which  successive  interference 
cancellation  at  the  transmitter  can  be  used  to  remove  the  multiuser  interference.  However,  in  the 
presence  of  estimation  errors  in  the  channel  matrix,  the  multiuser  interference  cannot  be  perfectly 
cancelled.  We  derive  the  effect  of  channel  estimation  errors  on  the  error  rate  performance  of  the 
receivers  and  assess  the  performance  degradation  on  the  system  as  a  function  of  the  quality  of 
channel  estimate.  Computer  simulations  are  also  performed  to  corroborate  the  analytical  results. 

In  Chapter  9  of  the  book  on  New  Directions  in  Wireless  Communications  Research  [2],  we 
present  a  tutorial  survey  of  interference  suppression  techniques  in  point-to-multipoint  MIMO 
systems  in  which  the  interference  suppression  is  performed  at  the  transmitter.  More  specifically, 
we  consider  the  downlink  of  a  MIMO  system  which  transmits  simultaneously  to  multiple  users 
that  are  geographically  distributed.  The  transmitter  employs  Nt  antennas  to  transmit  to  K 
receivers,  where  Nt  >  K.  Each  user  is  assumed  to  have  a  receiver  with  one  or  more  receive 
antennas.  The  distinguishing  feature  of  the  MIMO  downlink  transmission  scenario  is  that  the 
geographically  separated  receivers  do  not  employ  any  coordination  in  processing  the  received 
signals.  Consequently,  it  is  desirable  to  suppress  the  multiuser  interference  by  pre-processing  the 
signals  to  the  transmitter. 

The  suppression  of  the  multiuser  interference  by  means  of  transmitter  processing  is  usually 
called  signal  precoding.  Signal  precoding  at  the  transmitter  may  take  one  of  several  forms, 
depending  on  the  criterion  or  the  method  used  to  perform  the  precoding.  The  precoding 
techniques  described  in  this  chapter  include  linear  processing  methods  based  on  the  zero-forcing 
and  minimum  mean-square-error  criteria,  and  nonlinear  processing  methods,  namely, 
Tomlinson-Harashima  precoding,  vector  precoding,  and  a  precoding  method  based  on  lattice 
precoding.  An  assessment  of  the  effectiveness  of  these  linear  and  nonlinear  precoding  techniques 
in  suppressing  multiuser  interference  is  also  provided. 
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Channel-aware  seheduling  has  emerged  as  a  promising  teehnique  for  improving  speetrum 
effieieney  in  wireless  networks.  In  [1],  an  ad-hoe  network  was  evaluated  that  eontends  for  a 
ehannel  using  random  aeeess  in  eonjunetion  with  ehannel-aware  distributed  opportunistie 
seheduling  (DOS).  With  DOS,  onee  a  sueeessful  ehannel  probing  has  been  made  (through  a 
sueeessful  ehannel  eontention),  the  sueeessful  link  may  deeide  to  eontinue  the  data  transmission 
if  the  observed  ehannel  eondition  is  “good”;  otherwise,  it  may  skip  the  transmission,  and  let  all 
the  links  re-eontend  for  the  ehannel.  For  time  varying  ehannel  eonditions,  different  links 
experienee  different  ehannel  eonditions  in  different  time  slots.  It  is  likely  that  after  further 
probing,  the  ehannel  ean  be  taken  by  a  link  with  a  better  ehannel  eondition,  resulting  in  higher 
network  throughput.  Henee,  the  multiuser  diversity  aeross  links  and  the  time  diversity  aeross 
slots  can  be  exploited  in  a  joint  manner. 


In  a  practical  system,  many  wireless  applications,  such  as  multimedia  traffic,  have  stringent 
delay  requirements.  In  this  work  we  have  focused  on  developing  channel-aware  distributed 
opportunistic  scheduling  (DOS)  that  takes  delay  into  consideration.  Without  delay  constraints,  it 
is  possible  that  the  system  may  spend  an  arbitrarily  long  period  of  time  on  channel  probing, 
looking  for  better  channel  conditions.  This  may  significantly  degrade  the  QoS  performance  of 
delay-sensitive  applications,  for  which  the  delay  performance  is  of  critical  importance.  A  main 
objective  of  this  study  is  to  obtain  a  rigorous  understanding  of  DOS  under  delay  constraints. 


In  [I],  we  study  DOS  under  the  average  delay  constraint  from  two  different  perspectives, 
namely,  network-wide  average  delay  constraint  and  user-specific  average  delay  constraint. 
Specifically,  average  delay  constraint  here  refers  to  an  ensemble  constraint  after  averaging  over 
many  packet  transmissions. 


The  average  delay  constraint  was  first  considered  from  a  network-centric  point  of  view,  where 
links  cooperate  to  maximize  the  overall  network  throughput  subject  to  the  constraint  that  the 
network-wide  average  probing  and  transmission  time  is  no  greater  than  a  given  time 
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constant,  a  .  Observe  that  the  optimal  seheduling  under  the  network-wide  average  delay 
constraint  is  equivalent  to  a  constrained  optimal  stopping  problem.  In  this  work,  the  "stochastic" 
Lagrangian  approach  is  utilized  to  convert  the  constrained  problem  into  an  unconstrained  one. 
Under  the  assumption  of  continuous  channel  statistics,  it  is  shown  that  the  continuity  of  channel 
fading  distribution  ensures  that  the  channel  exhibits  sufficient  randomization  to  close  the  gap 
between  the  constrained  primal  problem  and  the  uneonstrained  dual  problem.  We  then 
eharaeterize  the  eorresponding  threshold-based  optimal  scheduling  algorithm  and  its  throughput. 
It  was  diseovered  that  a  sharp  transition  exists  for  the  optimal  threshold  tied  to  a  critical  time 
eonstant,  a  ,  in  the  sense  that  if  a  is  less  than«* ,  the  optimal  threshold  is  upper-bounded  by  a 
funetion  of  a  ;  otherwise,  the  imposed  delay  eonstraint  has  no  impaet  on  the  optimal  seheduling, 
and  the  optimal  threshold  is  the  same  as  that  in  the  unconstrained  case.  An  iterative  algorithm 
was  developed  to  find  the  eritieal  time  eonstant,  a  ,  and  show  its  eonvergenee. 


Fig.  3:  Maximal  throughput  for  different  values  of  a  at  different  average  SNR,  p 


Figure  3  illustrates  the  maximal  throughput  under  different  time  constraint  values,  a  at  different 
average  SNR,  /?  =  5,10,15,20,25.  In  this  figure,  the  throughput  of  the  threshold  policy,  is 
plotted  in  solid  lines  and  the  throughput  of  the  optimal  threshold  policy  for  the  unconstrained 
case,  X  is  plotted  by  dashed  lines.  The  circle  denotes  the  maximal  throughput  corresponding  to 
the  optimal  scheduling  policy  under  the  average  delay  constraint,  x*  .  We  first  observe  from 
Fig. 3  that  in  general,  <  x*  <  x* .  Another  important  observation  from  Fig.  3  is  that  x^ 
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intersects  with  x*  at  its  peak  point,  which  is  exactly  the  critical  time  constant,  a* .  Beyond  this 
point,  the  optimal  scheduling  policy  does  not  need  to  take  the  imposed  time  constraint  into 
account,  i.e.  the  optimal  scheduling  policy  with  time  constraint  will  be  the  same  as  that  without 
time  constraint  when  the  time  constraint  a  is  beyond  this  critical  point.  It  can  also  been  seen 
from  Fig.  3  that  the  maximal  throughput  x*  is  an  increasing  function  of  average  SNR,  p  for  a 

given  «  .  This  is  because  when  the  average  SNR  increases,  the  transmission  rate  becomes  higher 
and  yields  a  higher  throughput. 


Next,  we  explore  distributed  scheduling  under  the  average  delay  constraint  from  a  user-centric 
perspective,  where  each  link  seeks  to  maximize  its  own  throughput  subject  to  the  constraint  that 
the  expected  user  delay  should  be  no  greater  than  its  own  delay  constraints^,  m  =  1,2,...M  .  We 

treat  the  problem  threshold  selection  under  the  individual  delay  constraint  as  a  non-cooperative 
game.  We  show  that  the  Nash  equilibrium  exists  for  this  constrained  non-cooperative  game.  Our 
results  reveal  that  there  exists  a  vector  of  critical  time  constant  ,  as  the  counterpart  to  that  for 

the  network-centric  case,  such  that  only  when  all  the  delay  constraint  >  s* ,  then  the  delay 
constraints  have  no  impacts  on  the  optimal  threshold  in  the  Nash  equilibrium.  We  further  provide 
an  iterative  algorithm  to  find  the  Nash  equilibrium  and  prove  the  convergence. 


Table  3 

OPTIMAL  THRESHOLD  OF  E.ACH  INDIVIDUAL  LINK  UNDER  INDIVIDUAL  .AVER.AGE  TIME  CONSTRAINT 


Link 

Q* 

(pm  =0.125) 

r'* 

V?7l 

^  ’ 

V777 

Link  1 

66.12  (p  =  2) 

0.1504 

0.1504  (£Vm  =  67.0) 

0.077(am  =  64.0) 

Link  2 

64.81  (p  =  4) 

0.2188 

0.2188  (Qm  =  67.5) 

0.175(q>„  =  64.5) 

Link  3 

64.14  (p  =  6) 

0.2652 

0.2652  (a„  =  68.0) 

0.263(q>„  =  65.0) 

Link  4 

63.70  (p=  8) 

0.3006 

0.3006  (Qm  =  68.5) 

0.298(Qm  =  65.5) 

Link  5 

63.39  (p  =  10) 

0.3293 

0.3293  (am  =  69.0) 

0.326(Qm  =  66.0) 

Table  3  is  presented  to  show  the  optimal  threshold  of  each  individual  link  under  individual 
average  time  constraint  in  Nash  Equilibrium  state.  The  numerical  result  in  table  3  shows  that 
when  the  imposed  constraint  of  all  the  links  are  bigger  than  then  the  critical  time  constraint 

of  a* ,  the  optimal  threshold  of  all  links  with  time  constraint  (shown  in  4-th  column)  will  be  the 
same  as  the  scenario  without  time  constraint  (shown  in  3rd  column).  At  the  other  hand,  when  the 
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imposed  time  eonstraint  is  less  than  the  (shown  in  5-th  column),  the  optimal  threshold  of  all 
the  links  will  be  less  than  the  unconstrained  threshold,  . 


In  [1],  we  only  considered  the  average  delay  constraint.  It  would  also  be  of  great  interest  to 
consider  individual  packet  lifetime  constraint.  Different  from  the  average  delay  constraint  that 
specifies  an  ensemble  constraint  after  taking  average  over  many  packet  transmissions,  the  packet 
lifetime  constraint  is  imposed  on  scheduling  of  individual  packets  and  is  essentially  a  constraint 
on  sample-path  realizations.  We  are  currently  working  along  this  line,  and  had  some  initial 
results. 
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2.12  Bargaining  in  the  MISO  Interference  Channel 

PI:  A,  Lee  Swindlehurst 
GSR:  Matt  Nokleby 
Summary 

A  common  model  used  for  describing  wireless  ad  hoc  networks  at  the  physical  layer  involves  assuming  a 
collection  of  links  that  are  simultaneously  active  and  thus  interfering  with  one  another  during  a  given  time 
slot.  The  problem  becomes  one  of  allocating  the  resources  of  the  network  (power,  beamformer  weights, 
etc)  so  that  a  particular  performance  metric  (e.g.,  overall  throughput,  fairness,  etc.)  is  optimized.  In  this 
work,  we  have  examined  the  MISO  interference  channel  using  cooperative  bargaining  theory.  Bargaining 
approaches  such  as  the  Nash  bargaining  solution  and  the  Kalai-Smorodinsky  solution  have  previously 
been  used  in  wireless  networks  to  strike  a  balance  between  max -sum  rate  efficiency  and  max -min  equity 
(fairness)  in  users’  rates.  However,  cooperative  bargaining  for  the  MISO  interference  channel  has  only 
been  previously  studied  extensively  for  the  two-user  case.  We  have  developed  an  algorithm  that  finds  the 
optimal  Kalai-Smorodinsky  beamformers  for  an  arbitrary  number  of  users  in  the  general  MISO 
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interference  channel  setting.  We  also  consider  joint  scheduling  and  beamformer  selection,  using  gradient 
ascent  to  find  a  stationary  point  of  the  Kalai-Smorodinsky  objective  function.  When  interference  is  strong, 
the  flexibility  allowed  by  scheduling  compensates  for  the  performance  loss  due  to  local  optimization. 
Finally,  we  have  explored  the  benefits  of  power  control,  showing  that  power  control  provides  non-trivial 
throughput  gains  when  the  number  of  transmitter/receiver  pairs  is  greater  than  the  number  of  transmit 
antennas. 

Background 

Many  of  the  recent  advances  in  multi-antenna  communications  systems  involve  multi-user  scenarios, 
which  present  a  complicated  problem  partially  because  performance  criteria  are  difficult  to  characterize. 
There  is  not,  for  example,  a  single  data  rate  or  bit-error  probability  to  optimize.  Instead,  we  can  only 
maximize  composite  performance  measures  such  as  the  network  sum  rate,  max -min  fairness,  or  quality- 
of-service  requirements.  Ultimately,  the  choice  of  objective  function  is  often  somewhat  arbitrary. 

To  meet  this  challenge,  researchers  have  begun  to  apply  game  theory  [1],  a  mathematical  idealization  of 
human  decision-making,  to  problems  in  multi-user  communications.  Game  theory  provides  a  systematic 
framework  for  the  study  of  decision-makers  with  potentially  conflicting  interests,  as  well  as  solutions  for 
such  conflicts.  Accordingly,  a  game -theoretic  analysis  can  provide  a  tractable,  structured  approach  to 
resource  allocation.  Researchers  have  employed  game-theoretic  ideas  to  design  “fair”  power-control 
schemes,  develop  decentralized  network  algorithms,  and  otherwise  solve  resource-allocation  problems  in 
communications  networks  [2-8]. 

In  this  work,  we  have  focused  on  the  multiple-input  single-output  (MISO)  interference  channel.  In  the 
MISO  interference  channel,  several  communication  links,  each  involving  a  multi-antenna  transmitter  and 
a  single-antenna  receiver,  are  simultaneously  active.  This  scenario  models,  for  example,  ad  hoc  networks 
with  multiple  antenna  nodes  where  the  receivers  employ  fixed  beamformers.  Game -theoretic  solutions  for 
the  MISO  interference  channel  based  on  bargaining  have  been  considered  in  [9,10],  but  only  for  the  two- 
user  case.  Our  particular  focus  has  been  to  maximize  network  performance  according  to  the  Kalai- 
Smorodinsky  solution  [1 1],  a  cooperative  bargaining  approach  closely  related  to  the  well-known  Nash 
bargaining  solution  [12].  For  our  problem,  the  fundamental  idea  of  the  Kalai-Smorodinsky  (K-S) 
approach  is  to  maximize  users’  rates  while  ensuring  that  users  experience  the  same  fraction  of  the  rate 
they  would  achieve  without  interference.  In  practice,  the  K-S  solution  defines  a  compromise  between 
efficiency  (defined  herein  in  terms  of  maximizing  the  sum  rate)  and  equity  (maximizing  the  minimum 
rate). 

Our  primary  contribution  is  an  algorithm  that  efficiently  finds  the  K-S  solution  for  an  arbitrary  number  of 
users,  rather  than  just  the  two-user  case.  We  transform  the  rate -maximization  problem  to  a  series  of 
convex  programming  problems,  allowing  us  to  find  the  beamformers  that  achieve  the  rates  defined  by  the 
K-S  solution.  A  drawback  of  the  K-S  solution  is  that,  when  interference  becomes  strong  for  a  single  user, 
all  users’  bargained  rates  tend  toward  zero.  To  avoid  this,  we  have  also  investigated  the  use  of  joint 
scheduling  and  beamformer  selection  under  K-S  bargaining,  which  introduces  a  temporal  degree  of 
freedom  for  avoiding  interference.  Scheduling  also  convexifies  the  feasible  rate  region,  which  is  an 
important  consideration  in  cooperative  bargaining.  However,  the  need  to  jointly  address  scheduling  and 
beamformer  selection  complicates  the  optimization,  preventing  us  from  easily  finding  the  K-S  solution. 
We  therefore  devise  a  gradient-based  algorithm  to  find  a  stationary  point  of  the  K-S  objective  function. 


84 


While  we  sacrifice  global  optimality  to  include  scheduling,  the  performance  advantage  of  employing 
time-division  multiplexing  significantly  outweighs  the  potential  loss  of  optimality  when  the  interference 
is  strong. 

Simulation  Results 

To  illustrate  the  performance  of  our  approach,  we  use  a  simulation  with  4  transmit  antennas  and  K  users. 
In  each  simulation,  we  randomly  place  K  transmitter/receiver  pairs  on  the  unit  square.  The  channel 
coefficients  are  independently  drawn  from  the  zero-mean,  unit-variance,  complex  Gaussian  distribution. 
The  channel  gains  are  computed  according  to  a  standard  path-loss  model  based  on  distance.  In  Figures  1 
and  2  we  examine  algorithm  performance  in  terms  of  efficiency  and  equity  for  K  =  {2,  4,  6,  8,  10}. 

We  compare  our  proposed  K-S  algorithms  with  rates  obtained  by  the  max-min,  max-sum,  and  TDMA 
algorithms.  The  TDMA  rates,  represent  a  baseline  for  the  scheduled  K-S  solutions.  Figure  5  shows  the 
average  mutual  information  per  user,  averaged  over  100  realizations  for  each  value  of  K.  Not  surprisingly, 
the  average  rate  is  highest  under  sum  rate  maximization.  Both  K-S  approaches  degrade  as  we  increase  the 
number  of  users,  but  eventually  the  scheduling  approach  gives  a  better  average  rate  in  spite  of  the  fact  that 
it  gives  only  a  stationary  point.  In  Figure  6  we  examine  the  minimum  mutual  information  across  all  links, 
averaged  over  the  same  100  realizations.  Max-min  (again  unsurprisingly)  gives  the  highest  minimum  rate, 
followed  by  the  K-S  approaches.  Max -sum  gives  the  worst  minimum  rate,  which  drops  nearly  to  zero 
beyond  K  =  2.  The  K-S  solution  allows  us  to  maintain  the  sum  rate  while  still  protecting  the  weakest 
links. 
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Figure  1 :  Average  mutual  information  per  user. 


Figure  2:  Average  mutual  information  of  the  worst-case  user. 
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2.13  Power  Allocation  and  Bit  Loading  for  Spatial  Multiplexing  with 
Imperfect  CSI 

PI:  A,  Lee  Swindlehurst 
GSR:  Michael  Larsen 
Summary 

In  narrowband  multiple-input  multiple-output  communication  systems  when  the  channel  state  information 
(CSI)  is  known  perfectly  at  the  transmitter  and  the  receiver,  and  with  no  restrictions  on  codebook  design, 
the  well-known  water-filling  solution  maximizes  information  throughput  given  a  fixed  transmission 
power  through  use  of  the  singular  value  decomposition  (SVD).  In  practice,  when  practical  codebooks  are 
used  and  perfect  CSI  is  unavailable  due  to  mobility  and  channel  estimation  errors,  adaptations  to  power 
levels  and  bit-loading  schemes  may  need  to  be  made  in  order  to  maintain  reasonable  performance.  To 
analyze  this  scenario,  we  have  derived  simple  expressions  that  detail  the  impact  of  imperfect  CSI  on  the 
signal  and  interference  plus  noise  powers  for  subchannels  obtained  via  the  SVD-obtained.  These 
expressions  are  used  to  derive  thresholds  on  the  level  of  CSI  imperfections  and  noise  power  tolerable  in 
SVD-based  multiplexing  systems  using  M-ary  quadrature  amplitude  modulation.  In  particular,  using  our 
analysis,  one  can  determine  in  advance  whether  or  not  spatial  multiplexing  will  provide  any  rate 
performance  gain  compared  with  simple  beamforming  solutions.  In  addition,  approximately  optimal 
subchannel  power  levels  are  found  for  the  imperfect  CSI  scenario. 


87 


Background 


Given  perfect  channel  state  information  (CSl)  at  the  transmitter  and  the  receiver  of  a  multiple-input 
multiple-output  (MIMO)  communications  link,  along  with  no  restrictions  on  codebook  design,  the  well- 
known  water-filling  solution  maximizes  information  throughput  given  a  fixed  transmission  power.  The 
water-filling  solution  relies  on  the  use  of  the  singular  value  decomposition  (SVD)  to  separate  the  MIMO 
channel  into  independent  single-input  single-output  subchannels,  which  separation,  in  principle,  enables 
interference-free  data  multiplexing  [1-2].  Based  on  the  signal-to-noise  ratio  (SNR)  of  the  resulting 
subchannels,  bit-loading  and  coding  schemes  are  then  devised  to  approach  the  available  channel  capacity. 


In  practice,  two  difficulties  arise  that  must  be  considered  when  designing  an  SVD-based  communication 
scheme.  First,  it  is  not  possible  to  obtain  perfect  CSl  since  realizable  channel  estimates  are  formed  from 
noisy  measurements.  Mobility  in  the  communication  system  adds  additional  complications  for  CSl 
estimation,  since  the  MIMO  channel  changes  rapidly  when  the  transmit  or  receive  arrays  are  in  motion. 
When  noisy  or  outdated  CSl  is  used  in  conjunction  with  an  SVD-based  multiplexing  method,  the  MIMO 
subchannels  become  coupled  resulting  in  potentially  severe  subchannel  power  loss  and  mutual 
interference.  As  a  consequence,  subchannel  bit-loading  levels  selected  based  on  the  assumption  of  perfect 
CSl  may  no  longer  meet  given  error  performance  constraints.  For  these  constraints  to  be  met,  back-off 
strategies  are  often  employed  to  reduce  the  subchannel  bit  rates  and  create  an  error  margin  that  accounts 
for  the  CSl  mismatch. 

However,  without  knowledge  of  the  levels  of  the  signal  loss  and  interference  induced  by  the  subchannel 
coupling,  it  is  not  clear  how  to  systematically  carry  out  such  reductions. 


A  second  difficulty  arises  from  the  finite  nature  of  the  symbol  constellations  used  for  bit  allocation.  In 
order  for  the  optimality  of  the  water-filling  power  levels  to  hold,  subchannel  bit  loading  requires  an 
infinite -length  codebook  with  continuous  modulation  order.  The  use  of  discrete  codebooks  and  finite 
constellations  such  as  quadrature  amplitude  modulation  (QAM)  results  in  additional  performance  loss 
when  used  with  the  power  levels  chosen  by  the  water-filling  method  which  were  selected  assuming 
idealized  symbol  constellations  [3]-[5].  When  practical  codebooks  are  to  be  used,  subchannel  power  level 
selection  incorporating  information  about  the  codebook  may  help  reduce  such  losses  due  to  mismatches 
between  theory  and  practice.  Capacity-optimal  transmission  strategies  when  only  imperfect  CSl  is 
available  have  been  studied  in  [6]-[9]  and  the  references  therein.  Minimum  symbol-error-rate  (SER) 
optimization  and  adaptive  modulation  have  also  been  studied  in  [10],  [11].  These  papers  generally 
advocate  multiplexing  or  beamforming  using  other  information,  such  as  channel  covariance  information, 
to  guide  their  design,  and,  unlike  SVD-based  schemes,  often  require  complex  decoding  schemes  at  the 
receiver. 
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Several  other  studies  seeking  to  retain  the  elegant  structure  of  MIMO  SVD-based  solutions  have  focused 
on  the  impact  that  imperfect  CSl  and  finite  modulation  have  on  the  SVD-based  schemes.  The  effect  of 
imperfect  CSl  on  the  symbol  error  rate  (SER)  of  MIMO  SVD-based  multiplexing  methods  has  been 
examined  in  [12]  for  BPSK  transmission,  and  its  effect  on  the  bit  error  rate  (BER)  using  more  general 
constellations  including  M-ary  QAM  (M-QAM)  has  been  analyzed  in  [5],  [13]-[16].  With  an 
understanding  of  the  effect  of  imperfect  CSl  and  finite  modulation  schemes  on  SVD-based  multiplexing 
methods,  the  primary  design  task  becomes  the  selection  of 

subchannel  power  and  bit-loading  levels.  For  instance,  having  derived  an  expression  for  the  BER,  [14] 
proposes  an  ad-hoc  method  to  select  power  and  bit-loading  levels  in  each  subchannel  to  meet  a  chosen 
BER  constraint.  For  the  special  case  of  uniform  power  allocation  over  all  of  the  MIMO  subchannels,  [5] 
develops  an  expression  for  an  adjusted  per-subchannel  signal-to-interference-and-noise  (SINR)  ratio  that 
may  be  used  to  select  subchannel  bit-loading  levels  [17].  In  [4],  optimal  bit  loading  using  finite  symbol 
constellations  assuming  perfect  CSl  is  discussed  and  an  ad-hoc  method  is  proposed  for  approximating 
that  optimum. 


We  have  developed  a  closed-form  expression  for  the  subchannel  SINR  when  using  imperfect  CSl  in 
SVD-based  signaling  schemes.  This  result  generalizes  the  SINR  derivation  of  [5]  allowing  for  non- 
uniform  power  allocation  and  correcting  an  overestimation  of  the  interference  power  due  to  CSl 
inaccuracies.  In  addition,  this  new  expression  allows  for  different  levels  of  CSl  error  at  the  transmitter  and 
receiver,  whereas  previous  methods  typically  assume  the  CSl  is  identical  at  both  ends  of  the 
communications  link  or  that  perfect  CSl  is  available  at  the  receiver.  Once  obtained,  our  SINR  expressions 
allow  us  to  (1)  analyze  the  impact  of  imperfect  CSl  on  the  MIMO  SVD  subchannels;  (2)  find  thresholds 
for  the  amount  of  channel  uncertainty  and  measurement  noise  tolerable  when  multiplexing  over  the  SVD- 
based  subchannels;  and  (3)  perform  robust  bit-loading  design  through  the  systematic  selection  of  the 
subchannel  power  levels.  For  the  special  case  of  identical  CSl  at  both  ends  of  the  link,  the  SINR 
expression  allows  us  to  derive  approximately  optimal  subchannel  power  levels  for  M-QAM  signaling 
under  a  subchannel  SER  constraint,  a  result  directly  applicable  to  robust  bit-loading  design. 
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Figure  3  -  A  comparison  of  the  average  bits  per  transmission  achieved  using  the  uniformed  and  informed 

bit-loading  schemes  when  accounting 

for  subchannel  outages  with  confidence  level  0.1. 


Simulation  Results 

In  Figure  3,  we  show  a  sample  simulation  result  to  illustrate  the  validity  of  our  analysis.  We  assume  a 
four-antenna  transmitter  and  a  four-antenna  receiver,  and  we  assume  that  both  the  transmitter  and  receiver 
are  using  an  outdated  channel  estimate  whose  coefficients  are  equal  to  the  true  coefficients  plus  a 
Gaussian  perturbation  term  of  a  given  variance.  Note  that  our  analysis  holds  for  more  general  scenarios, 
where  the  CSI  at  the  transmitter  and  receiver  is  not  the  same,  but  this  example  will  illustrate  the  point. 

We  compare  the  average  throughput  (capacity  x  (1  -  outage  rate))  for  four  different  schemes:  (1) 

WFUBL,  which  is  the  standard  water-filling  solution  that  ignores  the  presence  of  the  CSI,  (2)  WFIBL, 
which  is  the  approach  detailed  in  [14]  that  takes  the  CSI  error  into  account,  (3)  BFIBL,  which  is  the 
performance  of  standard  beamforming  with  transmission  power  modified  to  account  for  the  CSI  error, 
and  (4)  UPIBL,  which  does  optimal  power  allocation  and  bit  loading  taking  CSI  into  account.  The 
average  bits  per  transmission  achieved  by  these  schemes  are  plotted  in  Figure  3  versus  the  variance  of  the 
CSI  error.  The  vertical  blue  line  is  the  threshold  derived  by  our  analysis  showing  the  level  of  CSI  error 
above  which  it  is  no  longer  advantageous  to  perform  spatial  multiplexing.  The  threshold  accurately 
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predicts  the  breakpoint,  which  occurs  for  perturbations  that  are  about  13dB  below  the  average  gain  of  the 
channel  coefficients.  Note  that  our  proposed  optimal  power  allocation/bit  allocation  scheme  (UPIBL) 
provides  the  best  performance  across  all  levels  of  CSl  error. 
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2.14  Instantaneous  and  Average  Sum-Rate  Optimizations  in  Multi¬ 
user  MIMO  Subnetworks 


James  R.  Zeidler,  Principal  Investigator 
Mike  Jensen,  Principal  Investigator 
Adam  Anderson,  Graduate  Student 


The  sum-rate  throughput  was  analyzed  with  linear  processing  (beamforming)  in  multiple-input 
multiple-output,  single-  and  multi-user  wireless  links  for  subnetworks  in  the  context  of  larger 
mobile  ad-hoc  networks  (MANET). 


The  complexity  that  MANET  scheduling  algorithms  will  encounter  for  various  types  of  PHY 
layer  links  was  determined.  This  included  point-to-point  (P2P)  links  including  the  single  user 
channel  (SUC),  broadcast  channel  (BC),  multiple  access  channel  (MAC),  interference  channel 
(IC),  and  also  for  more  complex  topologies  generally  referred  to  as  hybrid  channels  (HC)  that 
includes  the  MIMO  X  channel.  The  analysis  first  derives  the  number  of  all  possible  half-duplex 
link  configurations  for  a  given  MANET  link  topology,  and  then  derives  the  sum-rate-optimal 
beamforming  for  each  of  the  topologies  using  explicit  channel  state  information  for  all  users. 
Figure  1  shows  the  number  of  possible  schedules  as  a  function  of  the  number  of  nodes  in  the 
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network.  This  plot  shows  that  the  eomplexity  of  seheduling  grows  signifieantly  as  the  network 
size  inereases. 


Maintaining  aeeurate  CSI  for  all  users  beeomes  a  signifieant  eomplexity  issue  for  networks  with 
time  varying  links.  When  the  input  information  neeessary  to  derive  these  optimal  beamformers 
beeomes  outdated  or  erroneous,  signifieant  loss  at  the  PHY  layer  is  seen  whieh  would  have  a 
easeaded  loss  effeet  on  the  higher  network  layers.  Figure  2  shows  the  loss  in  sum-rate  for  CSI 
based  beamformers  when  the  transmitting  nodes 


Fig.  1 :  Number  of  possible  sehedules  for  various  links  in  an  N  -  node  ad-hoe  network 


(a)  Equal  Per-Node  Power  (ENP) 


Fig.  2:  Instantaneous  sum-rate  throughputs  with  erroneous  CSI  at  the  transmitter 
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have  outdated  CSI  that  is  quantified  by  the  ehannel  varianee  over  time,  .  Measured  ehannel 

eharaeteristies  from  the  BYU  data  base  are  used  in  the  simulations  to  ensure  aeeurate 
representation  of  the  ehannel  spatial  eorrelation. 


An  alternate  beamformer  (RCDI-HC)  that  uses  ehannel  distribution  information  in  the  form  of 
spatial  correlation  was  formulated  to  provide  increased  stability  for  multi-user  communication  in 
these.  Since  the  channel  distribution  will  normally  remain  stable  for  a  longer  period  than  the 
instantaneous  CSI,  the  required  frequency  of  channel  updates  is  reduced  accordingly.  This  may 
lead  to  more  robust  scheduling  performance  for  tactical  networks  with  routes  that  can  be 
maintained  over  longer  time  intervals.  It  should  be  noted  however  that  although  the  required 
channel  updates  are  less  frequent,  they  require  updates  of  the  channel  covariance  matrix  rather 
than  the  instantaneous  channel  state  and  are  thus  more  complex.  This  is  discussed  in  section  1, 
and  methods  to  reduce  the  complexity  of  the  covariance  matrix  updates  using  simple  models 
such  as  the  Kronecker  model  are  being  investigated  under  this  project. 


To  compare  the  performance  of  the  two  sum-rate  beamformers  for  different  link  topologies,  we 
first  summarize  performance  results  for  the  perfect  CSI  (RCI-HC)  beamformer  in  Table  1.  Next, 
we  examine  the  performance  of  the  RCDI-HC  sum-rate  beamformer  that  based  on  channel 
distribution  information  in  the  form  of  spatial  correlation  in  Table  2.  The  results  show  that  the 
relative  trends  in  sum-rate  for  the  different  link  topologies  are  similar  to  those  observed  in  CSI- 
based  beamformer  in  Table  1.  Furthermore,  the  average  sum  rate  for  RCDI-HC  beamformers  for 
each  link  is  roughly  70-80%  of  their  RCI-HC  counterparts.  One  of  the  most  important 
observations  regarding  the  data  in  Table  2  is  that  for  a  statistically  stationary  channel,  the 
resulting  sum-rates  are  not  a  function  of  channel  time  variations  and  will  not  suffer  performance 
loss  due  to  delay  in  feeding  back  information  to  the  transmitter.  This  means  that  while  RCDI-HC 
offers  lower  initial  sum-rate,  it  maximizes  the  duration  over  which  this  average  sum-rate  is  valid. 
Research  into  the  temporal  stability  of  the  CDI  based  beamformers  in  typical  channel  conditions 
is  in  progress. 


Table  1 
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Average  sum-rate  (bits/sec/Hz)  for  various  links  using  RCDI-HC 


Link 

Equal 
Per-Node 
Power  (ENP) 

Equal  Total 
Power  (ETP) 

Single-User  (SUC) 

7.22 

9.98 

Broadcast  (BC) 

8.47 

11.37 

Multiple-Access  (MAC) 

11.73 

11.71 

Interference  (IC) 

10.03 

1 1.00 

Hybrid  (HC) 

11.38 

12.56 

Table  2 

Average  sum-rate  (bits/sec/Hz)  for  various  links  using  CSI-Based  Precoding 


Link 

Equal 
Per-Node 
Power  (ENP) 

Equal  Total 
Power  ( ETP) 

Single-User  (SUC) 

10.00 

13.23 

Broadcast  (BC) 

12.04 

16.27 

Multiple-Access  (MAC) 

16.16 

16.24 

Interference  (IC) 

12.76 

13.76 

Hybrid  (HC) 

16.33 

18.07 

Maximum  Rale 

16.77 

18.12 

Nonlinear  Broadcast  (DPC) 

13.20 

17.66 

References 


[1]  A.  Anderson,  J.  R.  Zeidler  and  M.  Jensen,  “Instantaneous  and  Average  Sum-Rate 
Optimizations  in  Multi-user  MIMO  Subnetworks”,  submitted  to  IEEE  Transactions  on  Wireless 
Communications 


95 


2.15  Outage  Probability  Characterization  of  Multi-User  MIMO 
Networks 

James  R.  Zeidler,  Principal  Investigator 
Bhaskar  Rao,  Principal  Investigator 
Sagnik  Ghosh,  Graduate  Student 

1.  Introduction 


In  this  work,  we  are  interested  in  providing  stable  transmission  in  time -varying  multi-user 
MIMO  networks.  One  of  the  most  compelling  problems  in  rapidly  changing  networks  is  the  issue 
of  channel  feedback  to  all  of  the  nodes.  It  is  well  known  that  Channel  State  Information  at  the 
Transmitter  and  Receiver  (CSIT  and  CSIR,  respectively)  at  every  node  can  significantly  increase 
the  throughput  of  the  network.  However,  the  amount  of  feedback  required  for  this  is  considerable 
and  can  consume  significant  amounts  of  bandwidth  in  the  system.  Furthermore,  in  networks 
with  high  mobility,  by  the  time  all  of  this  information  has  been  communicated  to  every  node,  that 
information  is  already  outdated  due  to  the  movement  of  the  nodes.  In  this  task  we  have 
investigated  different  types  of  channel  information  that  are  more  robust  to  changes  in  the 
channel.  Specifically,  our  work  deals  with  the  channel  statistics,  which  vary  much  more  slowly 
than  the  channel  itself.  Furthermore,  it  is  possible  to  attain  statistical  information  a  priori  based 
on  the  location  of  the  nodes.  Our  goal  was  to  determine  how  Channel  Distribution  Information 
(CDI)  can  be  utilized  to  achieve  reliable  communication  for  mobile  ad  hoc  networks. 


In  previous  years,  we  have  studied  how  to  use  CDI  to  maximize  the  average  rate  in  MIMO 
Broadcast,  MIMO  Multiple  Access,  and  MIMO  hybrid  networks  [I].  This  year,  having  explored 
the  average  rate  maximization  problem,  we  moved  to  studying  the  problem  from  an  outage 
probability  perspective.  This  formulation  has  many  advantages.  With  CDI,  while  we  can  assure 
some  average  rate,  we  can  never  guarantee  that  a  communication  scheme  will  achieve  a  certain 
rate  100%  of  the  time  (which  would  require  perfect  CSI).  In  an  outage  probability  framework, 
we  are  willing  to  accept  this  reality,  but  we  quantify  the  outage  on  the  links  and  guarantee  a 
certain  rate  or  SNR  the  remainder  of  the  time.  This  framework  allows  us  to  have  better 
characterize  our  network  for  purposes  of  design  and  analysis. 


Our  work  thus  far  has  concerned  minimizing  the  total  power  consumed  by  the  network  while 
meeting  certain  communication  requirements.  This  problem  is  non-trivial,  as  we  cannot  simply 
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meet  our  requirements  by  blasting  power  out  of  all  our  transmitter  nodes.  This  will  eause  too 
much  interference  in  the  network  to  other  nodes,  so  we  must  intelligently  use  the  spatial 
characteristics  of  the  channel  to  minimize  this  interference.  For  networks  with  perfect  CSI  at  all 
the  nodes,  this  problem  has  been  studied  in  [2].  We  look  at  this  problem  if  we  know  the  CDI  at 
all  the  nodes,  which  brings  us  much  closer  to  implementing  such  a  system  in  a  practical  high- 
mobility  MIMO  network. 


2,  Summary  of  Achievements 


First  and  foremost,  in  any  discussion  about  outage  probability  (the  probability  that  we  meet  a 
target  SINR  on  each  link),  we  desire  an  appropriate  closed-form  expression  to  be  able  to  analyze 
and  optimize.  In  the  first  part  of  this  work,  we  consider  the  popular  Rayleigh  fading  model  on  the 
antennas,  taking  general  correlation  structures  into  account.  These  correlation  structures  manifest 
themselves  as  covariance  matrices  for  each  link,  and  we  assume  all  the  nodes  have  this 
information  (or,  equivalently,  that  a  central  arbiter  has  this  information  and  controls  all  the 
nodes).  We  further  assume  that  on  each  link,  we  utilize  linear  transmit  and  receive  beamforming. 
Under  these  assumptions,  we  derive  a  closed-form  expression  for  the  outage  probability  for  both 
the  single-user  and  multi-user  MIMO  problems. 


Once  we  have  our  desired  outage  expression,  we  look  at  minimizing  the  sum  power  in  the 
network  subject  to  constraints  on  the  outage.  This  means  we  want  each  link  to  meet  its  SINR 
threshold  x%  of  the  time,  where  we  are  able  to  set  x.  We  first  look  at  the  single-user  problem, 
where  only  the  channel  statistics  are  known  at  both  the  transmitter  and  receiver.  Such  a  scenario 
can  occur  when  the  nodes  are  moving  so  rapidly  that  even  the  receiver  has  a  difficult  time 
capturing  the  channel.  In  addition,  channel  estimation  at  the  receiver  usually  requires  a  pilot 
signal  as  well  as  significant  computational  complexity,  which  the  designer  may  wish  to  avoid.  In 
such  a  setup,  when  we  have  CDI,  we  are  able  to  solve  for  optimal  transmit  and  receive 
beamformers  and  give  the  power  required  to  meet  an  outage  constraint.  See  Fig.  1  below.  We 
can  see  from  the  plot  that  most  of  the  CDI  outage  curves  are  actually  below  the  CSI  curve 
(assuming  perfect  channel  knowledge  at  the  transmitter  and  receiver).  From  this,  we  can  see  an 
additional  reason  we  would  want  to  use  such  a  setup:  if  we  are  willing  to  accept  a  given  outage 
on  the  link,  we  can  actually  transmit  with  lower  power  using  CDI  than  if  we  knew  the  channel 
perfectly  (using  CSI). 
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Single-User  MIMO 


Figure  15 


We  then  moved  on  to  looking  at  the  multi-user  MIMO  network  problem  and  minimizing  the  sum 
power  in  the  entire  network.  Using  some  tight  bounds,  we  are  able  to  find  “very  good”  transmit 
and  reeeive  beamformers,  and  we  solve  for  the  optimal  power  allocation  for  each  node.  See  Fig. 

2  below.  We  can  see  that  our  GDI  outage  solution  gives  enormous  power  savings  over  the  CSI 
case,  even  moreso  than  for  single-user  MIMO,  for  low  SINR  thresholds.  The  tradeoff  here  is  that 
if  our  system  requires  a  higher  SINR  threshold,  our  solution  becomes  infeasible  (which  can  be 
seen  when  the  curves  go  to  infinity)  and  we  are  not  able  to  guarantee  the  outage  constraints.  In 
today’s  systems,  however,  the  SINR  thresholds  on  the  lower  end  of  the  plot  are  reasonable, 
showing  that  if  we  use  GDI  information,  with  some  outage  we  can  save  a  lot  of  power  in  our 
system.  This  surprising  result  in  both  the  single-user  and  multi-user  MIMO  setups  makes  the 
study  of  the  outage  probability  framework  interesting  since  we  can  now  make  a  quantifiable 
tradeoff  between  outage  and  power.  The  two  additional  benefits  of  this  framework  is  that  all  this 
information  must  be  calculated  only  once  while  the  channel  statistics  are  valid  (compared  to  the 
channel  itself  being  valid),  yielding  huge  savings  in  computational  complexity,  and  as  discussed 
above,  the  large  reduction  in  feedback. 
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Multi-User  MIMO,  L=3:  Knowing  the  true  channel  vs.  the  statistics 
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2.16  Wideband  MIMO  Channels  and  Cubature  Kalman  Filters 


PI;  Simon  Haykin 

Graduate  Students:  Nelson  Costa  and  I  Arasaratnam 

Under  the  Canadian  node  of  the  MURI  Project,  we  have  made  significant 
contributions  to  the  MURI  Project  in  two  areas: 

1 .  MIMO  Wideband  Channel  Models,  with  three  contributions  spearheaded 
by  Nelson  Costa  and  Simon  Haykin,  as  highlighted  in  the  material 
embodied  under  the  following: 

(i)  Invited  Lecture  : 

N.  Costa,  and  S.  Haykin,  “A  correlation -based  wideband  MIMO  channel 
model,”  Proc.  General  Assembly  of  the  International  Union  of  Radio 
Science  2008  (URSI 2008), 

7-16  Aug.  2008,  Chicago,  IL.  (invited) 

(ii)  Invited  Tutorial: 

N.  Costa,  "  Efficient  and  Accurate  Matrix  Channel  Models," 

The  Annual  MURI  Meeting,  UCSD,  October  2008. 
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(iii)  New  Book: 

N.  Costa  and  S.  Haykin,  Modeling  of  Multiple -Input,  Multiple -Output 
Channel  Models,  accepted  for  publication,  Wiley,  2009. 

2.  Cubature  Kalman  Filters,  representing  a  new  generation  of 
nonlinear  filters,  which  is  the  closest  approximation  to  the  optimal 
Bayesian  filter  in  a  second-order  sense  and  under  the  Gaussian 
assumption,  exemplified  by  the  following  paper: 

I.  Arasaratnam  and  S.  Haykin.  "  Cubature  Kalman  Filters,"  IEEE  Trans. 
Automatic  Control,  June  2009 


2.17  Many-to-many  Communication  in  MANETs 

PI:  J.J.  Garcia-Euna-Aceves 

Over  the  past  year,  we  continued  our  work  on  many-to-many  communication  in  MANETs.  The 
findings  in  this  project  cover  two  areas: 

•  The  study  of  fundamental  limits  for  the  dissemination  of  information  over  MANETs 
when  many-to-many  communication  is  allowed. 

•  The  modeling  of  approaches  that  enable  many-to-many  communication  in  MANETs 
by  exploiting  processing  and  storage  complexity  in  mobile  nodes. 

1 . 1  Eundamental  Eimits  of  Wireless  Networks 

We  completed  our  work  on  the  characterization  of  the  optimal  interference-free  capacity  of  a 
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wireless  network  subjeet  to  unieast  traffic,  and  showing  how  this  order  capacity  can  be  attained 
in  the  presence  of  interference  by  means  of  [1].  The  importance  of  this  work  is  that  no  prior  work 
had  focused  on  first  establishing  what  is  the  optimal  capacity  of  a  wireless  network  in  the 
absence  of  MAI,  and  then  determining  whether  that  capacity  is  attainable  when  MAI  is  present. 

In  stark  contrast  to  the  existing  literature,  our  analysis  presents  the  possibility  of  actually 
increasing  the  capacity  of  ad-hoc  networks  with  the  number  of  nodes,  even  while  the 
communication  range  tends  to  zero. 


We  modeled  a  random  network  with  n  nodes,  a  homogeneous  communication  range  of 

r(n),  and  unicast  traffic  for  k  source-destination  (S-D)  pairs.  In  the  absence  of  interference,  such 
a  network  corresponds  to  a  random  geometric  graph  (RGG)  with  an  edge  between  any  two  nodes 
separated  by  a  distance  less  than  r(n).  We  defined  a  combinatorial  interference  model  based  on 
RGGs,  and  use  it  to  express  all  the  protocol  models  used  in  the  past  and  a  model  that  we  later  use 
to  show  that  the  optimal  capacity  of  a  wireless  network  is  indeed  attainable.  We  generalized  prior 
results  by  Gupta  and  Kumar  and  our  own  results  for  wireless  networks  with  MPR,  and  showed 
that  the  optimal  capacity  of  wireless  networks  is  attainable  in  the  presence  of  MAI.  We  showed 
that  MPTR  achieves  the  optimal  capacity  of  0(n  r  (n)/k).  This  constitutes  a  gain  of  0(nr  (n)) 
over  any  previously  reported  feasible  order  capacity. 


We  also  completed  our  work  on  the  order  capacity  of  wireless  networks  using  network  coding 
(NC).  We  had  previously  shown  that  NC  used  with  MPT  and  MPR  renders  the  same  order 
capacity  as  simply  using  MPT  and  MPR.  Our  new  results  [2,  3,  4]  show  that  NC  provides  no 
order  gain  for  multicasting  and  broadcasting  compared  to  what  can  be  attained  with  traditional 
multicasting  and  broadcasting  based  on  store-and-forward  routing.  Widely  cited  experiments  [5, 
6]  have  been  reported  recently  in  which  NC  has  been  used  successfully  in  combination  with 
other  mechanisms  to  attain  large  throughput  gains  compared  to  approaches  based  on 
conventional  protocol  stacks.  These  empirical  results  have  led  many  to  believe  that  the 
combination  of  NC  with  wireless  broadcasting  can  lead  to  significant  improvements  in  the 
multicast  order  throughput  of  wireless  networks.  However,  the  exact  characterization  of  the 
multicast  order  capacity  of  NC  in  wireless  networks  has  remained  an  open  problem  since  its 
introduction  ten  years  ago,  with  only  limited  results  having  been  reported  to  date  on  the  subject. 
We  undertook  the  characterization  of  the  multicast  and  broadcast  throughput  order  of  wireless 
ad-hoc  networks  in  presence  of  network  coding.  We  considered  a  network  consisting  of  n  nodes 
distributed  randomly  in  the  network  space,  with  each  node  acting  as  a  multicast  source  of  a 
group  of  m  randomly  chosen  nodes  in  the  network.  We  have  shown  that,  under  the  protocol 
model,  the  per-session  multicast  capacity  of  random  wireless  ad  hoc  network  in  the  presence  of 
arbitrary  NC  has  a  tight  bound  of  0(1  /A/mnlog(n))  when  m  =  0(  n  /log(n)  )  and  0(l/n)  when  m  = 
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Q(  n  /  log(n)  ).  In  addition,  we  showed  that,  under  the  physieal  model,  the  per-session  multieast 
eapaeity  of  random  wireless  ad  hoe  network  with  arbitrary  NC  has  a  tight  bound  of  0(1/  Vmn  ) 
when  m  =  O  (n  /log(n)  ),  and  0(l/n)  when  m  =  Q(n/log(n)).  It  has  already  been  established  in 
the  literature  that  the  above  bounds  are  achievable  on  the  basis  of  traditional  store-and-forward 
routing  methods.  Hence,  our  results  demonstrate  conclusively  that  the  throughout  gain  due  to  NC 
for  mutlicasting  and  broadcasting  is  bounded  by  a  constant  factor! 


Despite  our  negative  result  on  the  multicast  order  throughput  for  NC,  the  constant-factor  gains 
that  can  be  attained  in  some  cases  with  NC  over  store-and-forward  routing  should  not  be  ignored, 
and  they  may  be  of  importance  in  practical  settings.  However,  our  work  has  also  shown  that 
embracing  concurrency  at  the  physical  layer  by  means  of  MPR  and  MPT  provides  the  same 
multicast  order  throughput  with  or  without  NC,  and  that  MPR  and  MPT  can  be  used  to  attain  the 
optimal  unicast  capacity  in  wireless  networks  [1].  Hence,  the  effectiveness  of  using  a  scheme 
based  on  NC  to  attain  higher  multicast  throughput  should  be  compared  against  approaches  based 
solely  on  physical-layer  concurrency.  In  addition,  the  signaling  overhead  incurred  by  the 
different  approaches  should  be  evaluated  carefully.  Instantiating  the  necessary  NC  state  at 
relaying  nodes  may  incur  signaling  overhead  that  outweighs  any  potential  throughput  gains  that 
NC  is  intended  to  provide  in  the  forwarding  of  data  packets. 


The  last  part  of  our  research  on  the  capacity  of  wireless  networks  consisted  on  focusing 

on  techniques  that  can  attain  the  benefits  of  distributed  MIMO  systems  with  limited  cost. 

We  developed  and  studied  [7]  an  interference  management  technique  for  the  downlink  of  a 
wireless  cellular  network  with  which  D  (D  <  K)  independent  data  streams  can  be  broadcasted  to 
D  out  of  M  mobile  stations  with  single  antenna  such  that  these  data  streams  do  not  interfere  with 
each  other.  We  demonstrated  that  D  can  be  any  number  up  to  the  maximum  value  of  K,  as  long 
as  M  is  large  enough.  Therefore,  interference  management  is  capable  of  achieving  the  maximum 
multiplexing  gain,  as  long  as  there  is  a  minimum  number  of  mobile  stations  in  the  network. 
Surprisingly,  by  taking  advantage  of  fading  channels  in  multiuser  environments,  the  feedback 
requirement  to  transmit  K  independent  data  streams  is  proportional  to  K,  and  the  encoding  and 
decoding  scheme  is  very  simple  and  similar  to  that  of  point-to-point  communications. 


The  original  multiuser  diversity  concept  was  based  on  searching  for  the  best  channels 

to  use,  while  our  new  approach  shows  that  searching  simultaneously  for  the  best  and  worse 
channels  can  lead  to  significant  capacity  gains.  This  technique  can  asymptotically  achieve  the 
capacity  of  dirty  paper  coding  (DPC  )  when  M  ^  oo.  In  general,  we  can  have  D  mobile  stations 
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implementing  our  interferenee  management  seheme,  where  D  depends  on  the  number  of  mobile 
stations  in  the  network.  If  D  <  K,  then  the  rest  of  K 

-  D  mobile  stations  require  to  perform  cooperative  decoding  in  order  to  transmit  K  independent 
data  streams.  Our  proposed  multiuser  diversity  scheme  provides  a  tradeoff  between  multiuser 
diversity  and  cooperation  among  mobile  stations.  This  proposed  distributed  MIMO  scheme  does 
not  require  mobile  stations  to  cooperate,  as  long  as  there  are  enough  mobile  stations  in  the 
network.  It  achieves  optimal  K  maximum  multiplexing  gain  in  the  downlink  of  cellular  systems 
as  long  as  K  «  M  .  If  there  are  not  enough  mobile  stations  in  the  network,  partial  cooperation 
among  them  is  required  to  achieve  the  maximum  multiplexing  gain. 


1 .2  Modeling  of  Approaches  for  Many-to-Many  Communication  in  MANETs 


The  vast  majority  of  analytical  models  for  medium  access  control  (MAC)  protocols  (e.g., 

see  [8,  9])  have  assumed  an  ideal  physical-layer  model,  in  which  nodes  within  a  given 
transmission  range  receive  packets  with  probability  1  if  there  are  no  other  concurrent 
transmissions,  and  packet  transmissions  fail  with  probability  1  if  there  are  any  concurrent 
transmissions.  This,  of  course,  is  not  realistic  in  most  practical  situations,  especially  in  mobile  ad 
hoc  networks  (MANETs),  where  a  packet  can  be  successfully  received  when  the  received  power 
is  larger  than  a  given  threshold,  the  received  power  levels  may  show  significant  variations 
around  a  mean  power,  and  even  in  the  presence  of  multiple  access  interference  (MAI)  caused  by 
concurrent  transmissions,  receivers  can  decode  packets  with  probability  less  than  1 .  The 
decoding  probability  together  with  MAC  protocol  behavior  determine  the  performance  of  MAC 
protocols  in  MANETs. 


In  reality,  the  received  signal  is  a  combination  of  many  replicas  arriving  over  multiple  paths 
between  the  transmitter  and  receiver.  The  signal  on  these  different  paths  can  interfere  with  each 
other  constructively  or  destructively;  this  multi -path  effect  causes  the  received  signal  and  power 
to  become  variables  of  space.  Especially,  if  the  transmitter  or  the  receiver  moves,  channel  fading 
causes  the  received  signal  and  power  to  become  variables  of  time  as  well.  In  addition  to  above 
propagation  impairments  caused  by  imperfect  channel  conditions,  noise  and  MAI  caused  by 
other  concurrent  transmissions  also  impact  the  probability  of  correct  packet  reception.  Many 
MAC  protocols  employ  carrier  sensing  to  mitigate  MAI  by  listening  channel  status  before  packet 
transmissions.  As  a  result,  a  node  transmission  probability  depends  on  the  channel  conditions 
and  is  different  in  carrier-sensing  MAC  protocols  than  in  MAC  protocols  without  carrier  sensing. 
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While  the  importance  of  the  physical-layer  effects  on  the  performance  of  MANETs  is  well 
recognized,  most  prior  analytical  models  of  MAC  protocols  operating  in  MANETs  avoid  their 
characterization  for  the  sake  of  simplicity.  For  example,  Carvalho  et  al.  [10]  use  linear 
approximations  for  the  relationship  among  probabilities  of  the  channel  being  busy,  a  node 
transmitting,  and  a  packet  being  received  successfully,  and  do  not  account  explicitly  for  the 
effect  that  network  density,  node  mobility  and  other  physical-layer  factors  have  on  the 
performance  of  MAC  protocols. 


A  few  works  have  attempted  to  analyze  those  realistic  physical  layer  factors  and  incorporate  their 
effect  in  the  modeling  of  MAC  protocols.  Pham  et  al.  [11,  12]  have  tried  to  take  imperfect 
channel  conditions  into  consideration;  however,  they  only  analyze  IEEE  802.1 1  DCF  and  assume 
only  one  specific  channel  condition,  namely  Rayleigh  channel  fading.  In  addition,  they  assume 
that  a  transmission  fails  if  more  than  one  packet  is  transmitted  concurrently,  while  in  reality  a 
successful  packet  reception  is  determined  by  the  signal  to  interference  plus  noise  ratio  (SINR). 
Zheng  et  al.  [13]  generalized  the  effect  of  imperfect  channel  conditions  by  making  the  strong 
assumption  that  bits  are  transmitted  with  a  fixed  error  probability,  which  provides  a  wider 
application  area  for  the  model;  however,  their  assumption  is  also  not  practical,  because  there  are 
many  network  parameters  in  MANETs  that  can  impact  the  packet  reception  probability,  such  as 
traffic  rate,  network  density,  and  mobility. 


We  have  worked  on  developing  a  model  that  takes  into  account  imperfect  channel  conditions  and 
multiple  access  interference,  and  proposed  a  generalized,  parameterized  framework  for 
representing  the  interaction  between  the  physical  (PHY)  and  MAC  layers.  The  focus  in  this  work 
is  not  in  providing  an  exact  representation  of  specific  MAC  protocols.  Rather,  we  aim  to  model 
generic  MAC  protocols  in  which  either  an  efficient  interference  prevention  scheme  is  employed 
or  not,  and  extract  common  properties  of  those  protocols.  We  have  used  two  well  known 
examples  of  those  two  types  of  MAC  protocols,  namely  IEEE  802.1 1  DCF  and  Aloha,  and  verify 
the  correctness  of  our  models  via  simulations.  Our  first  results  have  been  submitted  for 
publication,  and  our  plan  is  to  apply  the  model  to  the  study  of  different  MIMO  schemes.  We  aim 
at  the  analysis  via  simulation  and  our  analytical  framework  of  specific  protocols  based  on  MIMO 
techniques.  More  specifically,  we  are  studying  the  behavior  of  AMPTR  (adaptive  multi-packet 
transmission  and  reception),  which  is  a  MAC  protocol  aimed  at  exploiting  MPR  and  MPT. 
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Conclusions 


The  primary  challenge  addressed  in  this  project  has  been  to  determine  how  to  best  utilize 
the  potential  gains  of  the  multiple  antennas  at  the  transmitters  and  receivers  to  improve 
the  overall  performance  of  the  tactical  mobile  ad  hoc  network.  One  of  the  areas  where 
significant  progress  has  been  made  is  the  design  of  MIMO  antennas  for  conditions  where 
the  aperture  in  which  the  antennas  reside  is  restricted  such  that  the  current  distributions 
for  the  different  elements  overlap  within  the  aperture. 

Multiple  antennas  at  the  nodes  of  the  network  have  been  shown  to  increase  link  capacity 
by  allowing  multiple  parallel  radio  channels  in  the  same  spectral  band.  This  spatial 
multiplexing  gain  can  provide  higher  link  capacity  and  facilitate  processing  of  multiple 
packets,  but  these  gains  must  be  balanced  against  the  increased  interference  that  is 
generated  in  a  multi-user  network.  It  is  possible  to  limit  the  interference  to  other  users  in 
the  network  by  utilizing  the  antennas  to  provide  beamforming  gain.  Alternatively,  the 
antennas  can  be  utilized  to  provide  diversity  and  coding  gain  using  space-time  coding 
techniques.  Tactical  networks  may  also  exploit  the  spatial  selectivity  of  MIMO  equipped 
nodes  to  null  unwanted  interference  from  jammers. 

One  of  the  areas  where  significant  progress  has  been  made  is  in  the  area  of  joint 
beamforming  and  space-time  coding  to  improve  overall  network  performance.  It  has 
been  firmly  established  that  optimal  network  performance  is  not  achieved  simply  by 
maximizing  the  link  throughput  since  increasing  the  point-to-point  capacity  does  not 
exploit  the  spatial  multiplexing  gains  that  can  be  realized  by  MIMO  equipped  nodes. 
Analytical  results  are  obtained  to  determine  when  spatial  multiplexing  provides 
performance  gains  relative  to  simple  beamforming  solutions  and  optimal  subchannel 
power  levels  are  derived  as  a  function  of  the  accuracy  of  the  channel  state  information 
available. 

The  manner  in  which  power  is  allocated  to  these  MIMO  nodes  in  space,  time  and 
frequency  has  been  analyzed  to  determine  the  best  approach  to  use  in  response  to  changes 
in  the  channel  state  and  the  network  topology.  In  addition,  the  feedback  schemes 
required  to  implement  beamforming  in  these  networks  has  been  analyzed  for  quantized 
feedback  systems  as  a  function  of  the  number  of  bits  of  information  provided  per  update. 
The  impact  of  performance  limitations  due  to  both  the  feedback  delay  and  the  channel 
estimation  errors  has  been  evaluated.  The  impact  of  channel  fading  was  quantified  for  a 
variety  of  modulation  and  coding  schemes. 

The  problem  of  allocating  the  network  resources  in  a  multi-user  MIMO  network  is 
dependent  on  the  specific  performance  metric  of  most  interest.  The  allocation  of  power 
to  the  different  channels  is  defined  by  the  beamformer  weights  and  can  be  modified  to 
satisfy  different  network  metrics  such  as  the  overall  throughput,  fairness,  end-to-end 
delay,  etc.  In  this  project  we  have  examined  the  use  of  cooperative  bargaining  theory  to 
define  the  tradeoffs  between  the  max-sum  rate  efficiency  (throughput)  and  max -min 
equity  (fairness)  for  an  arbitrary  number  of  users  in  MISO  interference  channels.  Power 
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control  was  shown  to  provide  significant  throughput  gains  when  the  number  of 
transmitter/receiver  pairs  is  greater  than  the  number  of  transmit  antennas. 


MIMO  nodes  were  utilized  to  aehieve  unequal  error  protection  among  the  different 
components  of  a  sealable  video  waveform  in  order  to  provide  a  natural  way  to  protect  the 
most  important  components  of  the  video  waveform  at  the  expense  of  the  less  important 
eomponents.  This  was  aecomplished  by  the  use  of  a  base  layer  and  an  enhancement 
layer  where  spaee-time  eoding  is  used  to  proteet  the  base  layer  and  spatial  multiplexing  is 
used  for  the  enhancement  layer. 

The  incorporation  of  spatial  multiplexing  gains  provided  by  MIMO  equipped  nodes  into 
the  overall  network  design  results  in  an  inherently  crosslayer  solution  that  exploits  the 
physical  layer  information  in  the  MAC  layer  for  scheduling  and  routing  decisions.  This 
however  requires  a  elear  understanding  of  how  the  channel  state  information  (CSI)  is 
incorporated  into  the  MAC  layer  and  how  the  accuracy  and  temporal  stability  of  the  CSI 
affect  the  performance  of  the  network  protocols  used  to  make  seheduling  and  routing 
deeisions.  The  accuracy  of  the  CSI  obtained  is  impacted  by  the  temporal  variations  in  the 
ehannel  that  arise  in  realistic  environments,  espeeially  in  networks  with  significant 
mobility  of  the  nodes. 

In  addition,  there  is  an  inherent  difference  in  the  time  seales  over  whieh  CSI  and  routing 
updates  are  made.  In  this  project  the  use  of  channel  statisties  in  conjunction  with,  or  in 
place  of,  ehannel  state  information  has  been  evaluated  in  order  to  inerease  the  temporal 
stability  of  the  physical  layer  information  provided  to  the  upper  layers  of  the  network. 

We  have  obtained  recent  results  on  utilizing  ehannel  distribution  information  (CDI)  for 
beamforming  that  have  been  shown  to  significantly  reduce  the  frequency  that  channel 
information  must  be  feedbaek  to  the  transmitter.  When  CDI  is  used  in  place  of  CSI  on  a 
given  link,  the  network  saerifices  some  performance  in  maximum  throughput  for  that 
speeific  link  relative  to  that  attainable  with  perfect  CSI,  but  may  provide  improved 
performanee  when  there  are  significant  errors  in  the  CSI.  In  addition,  signifieant 
reductions  in  the  power  required  to  maintain  a  given  level  of  performance  are  aehieved. 
The  CDI  is  based  on  the  ehannel  statisties  and  ehanges  more  slowly  than  the  coherence 
time  of  the  channel.  Reducing  the  frequency  of  the  updates  required  has  been  shown  to 
signifieantly  reduce  the  overall  level  of  interferenee  created  to  other  users  in  the  network. 


During  the  past  year  we  have  developed  beamforming  algorithms  based  on  channel 
distribution  information  (CDI)  under  a  variety  of  ad  hoe  networking  seenarios  (e.g. 
interference  channel,  broadcast  channels,  multi-user,  and  hybrid  ehannels)  and  compared 
the  results  to  those  obtained  using  CSI.  The  inereased  temporal  stability  of  the  CDI 
beamformer  requires  fewer  updates  of  the  channel  statisties  than  the  channel  state,  but 
the  feedback  update  now  requires  information  defining  the  ehanges  in  the  ehannel 
covariance.  In  this  project  we  are  evaluating  various  matrix  compression  teehniques  to 
reduee  the  feedbaek  complexity  for  CDI  -based  signaling  approaches  so  that  a  realistic 
comparison  of  the  overall  system  eomplexity  of  these  teehniques  can  be  determined  in 
realistie  channel  eonditions. 
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Modeling  of  the  wideband  MIMO  ehannel  has  been  eondueted  under  this  projeet  and 
various  models  such  as  correlation  and  matrix  based  channel  models  have  been  analyzed. 
The  ability  to  predict  channel  variations  over  time  has  been  evaluated  using  linear  and 
nonlinear  prediction  techniques  and  also  by  using  cubature  Kalman  Filters.  A  MIMO 
system  that  incorporates  a  Tomlinson-Harashima  precoder  was  utilized  to  pre-process  the 
transmitted  information  so  as  to  suppress  the  interference  at  each  of  the  decentralized 
receivers.  A  QR-decomposition  of  the  channel  matrix  was  employed  to  obtain  an 
equivalent  channel  in  which  successive  interference  cancellation  at  the  transmitter  can  be 
used  to  remove  the  multi-user  interference.  Analytical  results  for  the  error-rate 
performance  were  obtained  as  a  function  of  the  channel  estimation  errors  and  were 
verified  by  computer  simulation. 

A  variety  of  routing  strategies  have  been  evaluated  during  the  past  year  to  incorporate 
physical  layer  information  into  protocols.  A  cooperative  paradigm  that  exploits  diversity, 
hybrid  automatic  retransmission  requests  and  MIMO  techniques  in  a  multi-hop  network 
scenario  has  been  developed.  In  addition,  decode-and-forward  physical  layer  network 
coding  architectures  have  been  developed. 

Testbeds  have  been  developed  to  evaluate  how  the  routing  layer  decisions  interact  with 
the  physical  layer  model  for  multi-hop  flows.  MAC  protocols  have  been  developed  that 
utilize  the  interference  cancellation  capability  of  multi-user  MIMO  systems  to  increase 
the  spatial  reuse  in  a  multi-hop  network  by  allowing  multiple  concurrent  communications 
in  a  region.  The  overall  multi-hop  network  has  been  decomposed  into  separate  sub¬ 
topologies  decoupled  from  each  other. 

Factors  such  as  the  queuing  delay  have  been  considered  in  addition  to  the  channel  quality 
and  channel  quality  distribution  to  define  delay-throughput  trade-offs  for  the  network. 
Opportunistic  routing  is  used  to  improve  delay  performance.  Closed  form  expressions 
for  the  error  performance  of  the  network  are  provided  as  a  function  of  the  channel  and 
source  statistics.  In  addition,  channel  aware  distributed  opportunistic  scheduling 
techniques  have  been  developed  for  delay  constrained  networks. 

Analytical  studies  have  been  made  to  determine  the  number  of  hops  that  should  be  used 
between  the  source  and  destination  nodes  of  an  ad  hoc  network.  Reducing  the  number  of 
hops  reduces  the  overall  network  overhead  but  this  requires  increased  transmission 
ranges  with  a  resulting  increase  in  the  level  of  interference  to  the  other  nodes.  Single  hop 
network  metrics  such  as  the  network  throughput  (the  product  of  the  spatial  density  of 
transmitters  and  the  link  throughput)  and  the  information  efficiency  (the  product  of  the 
transmission  distance  and  the  link  throughput)  were  analyzed  to  evaluate  the  impact  of 
space-time  coding  vs  spatial-multiplexing  for  selected  multiple  access  and  coding 
schemes.  The  importance  of  utilizing  the  channel  diversity  that  can  be  obtained  for  the 
MIMO  relays  in  conjunction  with  error  correction  coding  is  established. 

These  results  were  then  extended  to  multi-hop  networks  where  queues  are  added  to  each 
node  and  the  end-to-end  delay  and  throughput  performance  are  evaluated  for  idealized 
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backlogged  and  non-backlogged  networks.  Closed  form  expressions  for  these  idealized 
multi-hop  networks  are  obtained  that  allow  the  mean  end-to-end  delay  performance  to  be 
optimized  with  respect  to  the  number  and  placement  of  relays.  The  analytic  results 
obtained  were  verified  by  network  simulation.  The  rate  stability  of  these  networks  was 
also  evaluated  and  the  conditions  are  obtained  that  ensure  that  the  lengths  of  the  relay 
queues  do  not  become  unbounded  over  time,  thus  resulting  in  an  infinite  end-to-end 
delay. 

Finally,  to  address  the  inherent  problem  of  maintaining  ehannel  feedback  to  all  the  nodes 
of  a  time-varying  multi-user  network  we  have  evaluated  the  use  of  an  outage  probability 
(the  probability  that  we  meet  a  target  SINR  on  each  link)  as  a  metric  that  can  maintain 
stable  transmission  in  a  rapidly  time  varying  seenario.  CSI  allows  communications  for 
each  link  to  be  maintained  at  the  expense  of  increased  power  to  maintain  aeeurate  CSI  in 
poor  ehannel  eonditions.  By  accepting  a  specified  outage  probability  for  the  links  in  the 
network,  the  undesirable  interference  that  results  from  the  attempts  to  increase  the 
transmitted  power  to  maintain  eonnectivity  on  a  link  with  poor  channel  conditions  is 
avoided.  The  channel  distribution  information  (CDI)  eannot  ensure  successful 
eommunieation  at  all  times  for  each  link  in  a  time-varying  fading  channel,  but  it  can 
guarantee  a  certain  signal  to  interferenee  plus  noise  ratio  (SINR)  that  will  provide  a 
reliable  link  with  a  specified  outage  probability.  The  use  of  CSI  and  CDI  information  has 
been  analyzed  for  the  outage  probability  metric.  It  was  shown  that  the  use  of  CDI  results 
in  significant  reductions  in  computational  complexity,  feedback  and  transmitted  power. 
The  full  impaet  of  the  use  of  CDI  in  mobile  ad  hoe  networks  is  a  topic  of  current  research 
as  discussed  above. 
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Lee  Swindlehurst 
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Editor-in-Chief,  IEEE  Journal  of  Selected  Topics  in  Signal  Processing  Member 
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Journal  of  Ad  Hoc  Networks  (Nov.  2008);  “Elnderwater  communications  and  Wireless 
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2007  (awarded  at  ICC,  May  2008) 


Appendix  2:  Technology  Transfer 

Michael  Jensen 

We  continue  to  work  with  Rayspan  Corporation,  a  start-up  company  with  venture  capital 
funding,  to  develop  miniature  antenna  technologies  for  use  in  WLAN  systems.  We  have 
delivered  a  reconfigurable  antenna  topology  which  is  currently  being  evaluated  for  commercial 
application.  We  have  also  delivered  evaluation  software  to  assist  in  the  identification  of  suitable 
antenna  topologies  for  MIMO  and  diversity  communication. 

Srikanth  Krishnamurthy 

Professor  Krishnamurthy  is  currently  working  with  CISCO  to  apply  the  understanding  gained  via 
measurements  on  Space  Time  Block  Coding  and  Space  Division  Multiplexing  modes  of 
operation  with  MIMO  from  the  MURI  project  to  provide  improved  overall  throughput  by  auto¬ 
configuration  and  management  of  WLANs. 
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